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INTRODUCTION 


I  . 

The  objectives  of  this  program  are  to  extend  and  advance 
the  technique  of  nuclear  transmutation  doping  in  silicon  to 
a  technology  capable  of  successfully  compensating  the  residual 
boron  always  found  in  detector  grade  silicon,  to  develop  a 
basic  understanding  of  the  transmutation  process  in  silicon, 
to  investigate  the  effects  of  radiation  damage  which  always 
accompany  the  transmutation  process  and  to  devise  optimum 
techniques  for  the  removal  of  this  radiation  damage.  Although 
this  research  is  directed  toward  boron  compensation  in  extrinsi 
monolithic  focal  plane  array  material,  the  results  of  this 
technology  are  currently  being  utilized  in  a  number  of 
detector  applications  and  it  is  expected  that  additional 
device  applications  will  continue  to  be  discovered. 

Briefly,  the  transmutation  doping  process  involves 
irradiation  of  high  purity  silicon  in  a  nuclear  reactor  which 
provides  a  high  flux  of  thermal  neutrons.  Of  the  three  stable 
silicon  isotopes,  i  (  re  1  at  i  ve  abundance  of  .VI),  as  a  result 

of  thermal  neutron  capture  and  beta  decay,  is  the  only  isotope 
to  be  transmuted  into  an  impurity  dopant  ('^l’j.  The  nuclear 
details  of  this  process  have  been  reviewed  previously  and 
will  not  be  repeated  here.  ^  he  mention,  however,  that  the 

production  rate  of  phosphorus  has  been  determined  to  be 

-  4  "S  1 

1.676  x  10  P  atoms/cm'  per  thermal  ncutron/cin" . 

By  successful  measurement  of  the  neutron  fluence,  the 
phosphorus  concentration  added  by  the  irradiation  process 


can  bo  determined  to  a  very  high  degree  ol'  accuracv. 
Furthermore ,  since  the  range  of  thermal  neutrons  in  silicon 
is  verv-  long  compared  to  typical  ingot  dimensions  and  imicc 
the  target  Si  isotopes  are  randomly  and  uni  formic  division! 
the  resultant  dopant  distribution  is  far  superior  in  unilorei 
to  any  other  doping  process  iui»  in  existanee,  noth  macro 
scopically  and  microscopically  under  those  c  i  reams  t  nice  s  in 
which  the  transmutation  produced  dopant  coneon t ra i i on  is 
large  compared  to  the  residual  impuritv  conceit  t  i".i  i  i  on  s  in  the 
starting  material.  Idle  microscopic  uniformity  leads  io  an 
increase  in  reverse  bias  breakdown  voltage  and  farms  the 
basis  for  the  application  of  Nlii-Si  by  the  power  device  indu: 

for  the  extrinsic  lit  detector  arrays,  the  goal  is  -  a  use 
\T1)  technology  to  compensate  the  residual  boron  concent  r  it  i  on 

in  silicon,  wiiich  in  high  quality  float  cone  is  tvpicallv  of 

1  t 

the  order  of  It)  ”B/cm'  ,  without  compensating  tiie  deep  a  c  cep  1 0 
(either  (la  or  In)  which  are  res  pons  i  b  1  e  lor  the  photo,  uiduci: 
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Section  J]J  presents  a  description  of  the 
experimental  techniques  and  apparatus  employed 
throughout  this  work. 

Section  I\  presents  a  description  of  the 
experimental  results  and  a  discussion  ci 
possible  i nterpret at  ions. 

Section  V  summarizes  progress  and  accomplishment 


II.  TIlIiORMT  I  CAL  ANALYSIS 

The  theoretical  analysis  of  the  behavior  of  tlic  electrical 
parameters  of  silicon  as  a  function  of  transmutation  doping 
flucnce  presented  in  the  previous  final  report^  has  been 
extended  under  the  present  contract  and  will  be  presented 
in  this  section.  The  homogeneity  of  probe  resistivity 
as  a  function  of  flucnce  or  compensation  ratio  has  also 
been  greatly  refined  and  essentially  completed.  Ibis  model 
has  also  provided  a  means  of  calculating  the  maximum 
resistivity  which  can  be  obtained  by  NTU  compensation  before 
mixed  typing  occurs  given  the  initial  resistivity  fluctuation 
in  the  starting  material.  A  more  refined  calculation  of 
the  number  of  displacements  for  a  given  neutron  flucnce  is 
also  presented  which  takes  into  account,  in  a  more  realistic 
way,  the  energy  spectrum  of  the  neutrons  and  ionization 
energy  losses  of  the  secondaries. 

A.  Critical  flucnce  and _ I;  1  cot  r  i  cn  1  Parameter  f.xt  rema 

The  critical  flucnce  necessary  to  produce  NIP 
donors  sufficient  to  take  a  p-type  sample  to  intrinsic, 
i.c.,  exact  compensation,  has  been  defined  previously 
a  s 

n-  •  •  v'* 

where  K  is  the  phosphorus  production  rate  given  a- 
1  x  10  '  I’  atoms /cm'  per  n/cnC .  ' 

This  critical  flucnce,  which  is  a  sign'll  leant 
parameter  for  NIP  componsat  ion  problem:  ,  wa-  determined 


1 


previously  only  for  the  ease  where  the  starting 
material  initial  acceptor  concentration  was  sufficiently 
large  that  the  intrinsic  carrier  concentration  could  be 


neglected.  This  is,  of  course,  a  reasonable  assumption 
for  most  float  zone  material  available  today,  however, 
a  conceivable  experimental  technique  to  increase  the 
accuracy  of  NTD  compensation  using  a  series  of  successive 
irradiations  and  measurements  to  reach  the  desired 
compensation  ratio  would  ultimately  produce  material 
for  which  the  above  assumption  is  not  valid.  U'e  have, 
therefore,  extended  the  calculations  to  consider  this 
poss i b i 1  i  ty . 

When  the  initial  hole  concentration,  p  ,  before 

.VTH  doping  is  large  compared  to  the  intrinsic  carrier 

concentration,  n;,  from  Iiq .  (1)  and  p  =  \  ^  -  N(), 

=  P  /K  =  K  /K,  (  f 

c  1  o  p  o 

where  =  (eg  )  \  K  is  the  donor  production  rate 

and  i  is  the  initial  resist  ivi  tv.  This  is  the 
o 

approximate  relation  between  initial  resistivity  and 

critical  flucnce  used  previously. 

We  now  wish  to  express  this  flucnce  a  s  a  fund  i  <>  n 

of  n  in  an  exact  wav.  I  rom  the  usual  expression  I oi 
o  '  1 


resistivity,  =  (nenf)  +  pen  )  , 
electron  concent  rat  ion  using  tip  - 
expression  for  the  resistivity  in 
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K  p 

=  jr _ 

■> 

bn  ■  “  +  p' 
1  1 


where  b  =  u  / w  . 

n  p 

Solving  lor  the  hole  convent  r;it  ion ,  we  obtain 


I  (  h  / ,  )  +  !  (K  /r]“  -  -1  bn  -  |  ■  } 

IP  I'  1  i 


But  another  expression  for  p  can  be  obtained  using 
np  =  n .  “  and  the  neutrality  condition  p  ■*  X,,  =  n  •*  N , 


p  -  >,  (X.  -  X.  J  4  !  (X.  -  x  r  4  -in., ‘I  •  !  . 

I  1 '  ’  '  '  I 

liquations  (1!  and  ( 3  )  must  be  equal  and,  altei 
■  ubst  i  t  lit  i  ng  Iq.  ill  into  (31,  uc  obtain 

(  K  /  r  )  4  !(K  /r)2  -  4bn  tl  ^  |  ^  -  K^.  4  |  (  K<f  ^  <  d  n  .  ‘  ; 


Sett  in  it  K 


v  and  k  /,  =  x,  uc’  nisi:  to  solve  Id.  (I 


Id  i  Squaring  both  sides  el  lq .  (el,  uc-  I  i  ini  that 


ibn  .  “  -i  \  \ /  \  -  -  ibn  . 

I  v 


y  ( y  4  y  y  4  i  n  ;  ,  4  _n 


six  4  \/  x 


nil  .  ;  -i  _  n  . 

I  ) 


uherc  ue  have  used  Iq .  (id  in  the  !  a  ^  t  step  ahe\ 


•d  i  v  i  nr  I  o  i  \  ,  uc  I  i  nd 


i 


2(  b  +  n  n  •  ~ 

1 


This  equation  is  exact  and  analogous  to  the  approximate 
l.q .  (  2  )  . 

By  comparing  l:.q .  ("1  with  liq.  (.2),  it  becomes  clear 

that  the  second  term  in  l.q.  ( ~  1  is  the  intrinsic  carrier 
concentration  correction  to  the  approximate  1'q.  (2). 

The  percentage  error  in  critical  fluence  introduced  by 
neglecting  n.  can  therefore  be  estimated  from 

i 

"  error  =  100  x  2nd  term  of  l.q.  ( "’ 1  / 1|>  . 

he  have  plotted  this  error  in  Tig.  1  as  a  function 
of  the  initial  resistivity,  o  ,  of  the  undoped  starting 
material.  he  see  that  for  initial  resistivities  below 
do, 000  icon,  the  error  in  critical  fluence  is  below 
1'.’. .  figure  1  can  also  he  used  t  o  calculate  a 
correction  term  where  the  critical  fluence  is  estimated 
from  f q .  (21  by  the  relation 

1  =  1 K  / K ,  <  111  \) 

P  o 

where  the  correction  term,  \,  is  given  to 

good  approximat  ion  by  \;  rror  in  lig,.  I  100. 


i 


%  ERROR  IN  CRITICAL  FLUENCE 
BY  NEGLECTING 


Once  the  critical  flucnce  is  determined,  the 
extrema  and  singularities  in  the  electrical  parameters 
vs.  flucnce  (see  Pigs.  .V  h  of  Ref.  1)  can  he  determined 
relative  to  In  fact,  the  flucnce  difference  between 

<j>  .  and  the  fluences  corresponding  to  these  extrema  are 
intrinsic  constants  characteristic  of  the  semiconductor. 

The  maximum  in  the  resistivity  vs.  flucnce  has 
been  reported  previously.^  The  parameters  of  interest 
are  given  by  (liq.  (It)  of  Ref.  1) 


.  \Jb  =  n .  /p  / 
i  v  1  v  n 


b  KpX/bnj2  =  l-eni(wn/Wp),-’|'  1 


H  y  i 


(b  -  n/x/b 


where  we  have  defined  1 .  by  the  equation 

<i> .  =  n  •  /  K . 
i  i 

In  the  above  equations,  is  the  flucnce  necessary 
to  reach  the  maximum  resistivity  possible  which  occurs 
at  a  flucnce  which  is  slightly  lower  than  the  flucnce 
necessary  to  produce  exact  compensation.  When  the  NTP 
flucnce  is  equal  to  ,  the  annealed  resistivity  is  a 
maximum,  o  ,  given  by  I.cj .  (P)  and  corresponds  to  a  hole 

carrier  concentration  given  by  liq.  (8). 

An  inspection  of  figs.  I  and  S  of  Ref.  1  shows 
that  the  flail  coefficient  and  Mall  mob  i  1  i  t  v  vanish  at 


a  second  important  !  luence  which  we  shal  1  cal  1  i  . 

Previously,  the  Hall  coc  I' f  i  c  i  ent  s  ,  Rj  j ,  Hall  mobilities, 

U||,  and  Hall  carrier  concentrat  ions,  p^  and  n  j  j ,  were 
,  1 

C  ivc n  h v 


p  (  p“  -  n  j  " b  “  i  / e  ( p  “  +  n  .  'b  H  ,  (p-n) 


n  (  n  .  ‘  -  n“h“  1/c  (.n  •  +  n'b)  ,  m  p) 

l  ! 


uptP“  -  ni"b")/(p  +  n  j  "b  )  (.'Rj)--  0  ) 


.  n  (,n~b~  -  n  .  “  )  /b  ( n  b  +  n  . 


ir,.  m 


Pj,  =  (P~  +  n."b)-/p(p  -  n  .  " b ”  1 


~)  l  l 


n(|  ^  (n  b  +  n  .  “1  ~/n  (n"b  -  n  .  '  )  .  (1 

l‘ h e  equations  either  vanish  or  have  a  singularity  at 
a  hole  concentration  yiven  by 

n  =  n  .  b  .  (1 

*  ■"  i 

To  find  the  corresponding  I*  luence,  we  use  the  expression 

for  carrier  concentration  as  a  function  of  fluence  derive 

,  1 

pro  v ions  1  v 
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where  I 


fi  I  \  i  ir;  t  he  abe\  e 


0  c  P(K«frc) 


(19) 


and  substituting  Eq.  (16)  into  Eq .  (19),  we  find  the 

fluence  necessary  to  produce  this  singularity  is  given 
by  in  the  expression 


K 


(20) 


The  Hall  coefficient  reaches  a  maximum  on  either 
side  of  its  vanishing  point,  as  shown  in  Fig.  S  of 

Ref.  1.  These  maxima  can  be  found  by  differentiating 
p j j  and  n^  and  by  noting  that  p  and  n  do  not  contain 
singularities  at  these  two  points.  Therefore  dp/d$  f  0 
at  these  maxima  and  since 


dp. 


dp 


I  ■  «c. 


tlien  these  maxima  must  occur  for  dp^/dp  =  0  and 
diijj/dn  =  0.  Differentiating  Eqs.  (14)  and  (15)  then 
yield  the  solutions  for  the  actual  carrier  concentrations 
at  these  maxima 
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and 
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*- 
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ITT 
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lls  ini;  then  Hq.  (19),  and  its  n-type  analog.,  which  is 

obtained  from  liq.  (IS)  and  given  hv 
*> 

n .  "  -  n  “ 

i  =  --1— —  , 

Kn 


The  equations  of  this  subsection  and  the  equations 
of  Ref.  1  which  describe  the  behavior  of  p,  n,  ,  Rj|, 
and  U|j  as  a  function  of  NTH  fluence,  t,  form  the  has  i  s 
for  any  theory  relating  to  changes  in  room  temperature 


electrical  parameters  as  a  function  of  fluence.  In  the 
next  subsection,  the  behavior  of  resistivity  i  nhonogene i t y 
variations  with  changes  in  compensation  ratio  is  analysed 
for  NTH  compensation  using  the  basic  equations  presented 
in  this  subsection.  A  second  application  of  these 


basic  equations  is  also  presented. 


/ 


B. 


Resistivity  Homogeneity  Model 


Under  a  previous  contract,  the  variation  in 
resistivity  fluctuation,  Ap/p ,  was  determined 
experimentally  and  justified  by  an  approximate  cal¬ 
culation.^  The  theory  has  since  been  developed  completely 
and  will  he  summarized  in  this  section.  lor  a  more 
comprehensive  treatment  see  Ref.  4. 

he  begin  by  assuming  that  the  compensation  in  the 
highly  zone  refined  p-tvpe  detector  material  is 
negligible.  U:c  will  also  neglect  the  burn-up  of  boron 
by  neutron  transmutation.  he  further  assume  that 
any  doping  inhomogeneity  added  by  transmutation  is  very 
small  compared  to  the  boron  inhomogeneity  present  in 
the  starting  material.  These  assumptions  have  been 

justified  elsewhere  as  being  reasonable  for  detector 

.1 

grade  float  zone. 

Under  the  above  assumptions,  we  know  immediately 
the  behavior  of  the  impurity  inhomogeneity  as  a  function 
of  flu  e  nee,  <f> .  Th  is  is  illustrated  in  fig.  . 

Up  =  K 4>  is  the  impurity  concentration  of  phosphorus 
addl'd  by  \ID  and  At.  p  =  o.  :  s  t  he  imperil' 
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Corresponding  to  our  maximum  and  minimum  impurity 


concentrations  represented  schematically  in  1' i  g .  2 
are  minimum  and  maximum  resistivities.  These  resistivities 
can  he  measured  hy  spreading  resistance  or  four  point 
probing  alter  various  N’TD  fiuences  and  annealing  as  was 
presented  in  Ref.  1.  lor  a  particular  resistivitv 
trace  across  a  wafer,  it  is  easy  to  measure  [,  and 
p  in*  We  then  wish  to  calculate  how  these  two  quantities 
vary  with  fluence.  We  will  define  the  mean  of  the  maximum 
and  minimum  resistivities  and  the  maximum  fluctuation 
as 


p  =  (p 


max 


min 


)/: 


(25) 


and 


Ap  =  fp  -  p  ) 
u  max  mm 


(2M 


We  also  define  the  fraction 


Ap/p  = 


7  max 


m  j  n 


0  +  > 

max  1  min 


O' 


as  a  measure  of  the  relative  resistivity  fluctuation. 
We  now  need  to  invest  igate  the  behavior  el 
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where  <f>  is  the  NTH  thermal  neutron  fluence.  Since 
Q  4>  ami  K  and  R  are  essentially  constant,  we  wish 


to  calculate  the  behavior  of  r  .  .  ,  as  a  function 

max (min) 

of  Q. 

This  behavior  is  illustrated  schematically  in 

Fig.  5.  The  resistivity  increases  with  fluence  and 

compensation  ratio  until  a  maximum  is  reached.  This 

occurs  just  before  exact  compensation  (where  Q  =  1 

and  =  Njj)  .  The  carrier  concentration  decreases 

simultaneously  toward  n.  at  0=1.  After  exact 
'  i 

compensation,  further  NTH  doping  causes  a  decrease 
in  p  while  the  electron  concentration  increases.  The 
shaded  region  to  the  left  of  Q  =  1  represents  the 
original  resistivity  fluctuation  in  the  starting  material 
defined  by  and  p  while  t he  shaded  region  to 

the  right  of  Q  =  1  represents  this  same  fluctuation 


after  the  material  has  been  ovc rcompensated .  Both 

P  and  p  .  follow  the  same  resist  ivitv  vs.  fluence 
max  min 

curves  as  the  NTD  compensation  is  increased. 

To  calculate  these  effects,  we  define  the  mean 

boron  concentration  and  boron  fluctuation  as 

B  =  B  . 
jj  =  ni  a  x _ m  i  n 

a  n  d 
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l  p’.urc  a.  Schematic  representation  < 
the  ina v  i mini:  and  minimum  n 
t  ions  and  carrier  concent 


These  parameters  arc  constant  since  the  burn  up  of  boron 


is  negligible  compared  to  the  rate  ol  pho--  idioms 
production,  h.  these  parameters  can  he  calculated  I  rem 
''max  aiu'  ’min  * 01  T*’e  -starting  mate:  laJ  from  the 
fo 1 1 ow i ng  cquat  i  ons  ^  ' 
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lor  p ■  type  atul  the  negative  of  the  above  for  n-tvpe 
mater i al . 

from  the  above,  the  carrier  concentrations 


corresponding  to  r  (O')  and  p  .  (Q1  arc  clearlv 

1  '  max  '  min' 
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from  which  the  resistivities  can  he  determined  as 


IQ) 
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ns i ng  Pq .  (13  1  o f  Re f .  1  . 

We  can  now  easily  calculate  .p/e  from  Pq  .  (3S)  by 

allowing,  Q  to  vary  in  fq.  (31)  through  (38)  and 
remembering  our  definition  of  1  (now  as  a  Inaction 


of  O) 


Ai  ■/*  K"  r: 


,  (0  )  t  (  0  ) 

'  max  mm 

i  (Q)  ' 

'  max  '  min 


The  results  of  the  above  equal  ions  have  been 
calculated  numerically  and  are  represented  in  fig.  1 
for  i  p  1  vpe  sample  with  mean  resistivity  ot  J.'-tUt  i 


Several  curves  are  shown  corresponding  to  init  lal 
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resistivity  fluctuations,  100  Ap/p,  of  5,  20,  40,  and  60". 
Figure  5  shows  similar  results  except  that  the  mean 
resistivity  value  of  the  starting  material  was  taken 
to  be  10,000  fi-cm. 

The  following  points  should  be  noted  from  Figs. 

4  and  5.  Both  sets  of  curves  are  quite  similar 
indicating  that  the  results  are  reasonably  independent 
of  starting  material  resistivity.  For  p-type  material, 
the  fluctuation  always  increases  as  a  compensation 
ratio  of  unity  is  approached.  The  situation  would  be 
amplified,  however,  if  the  phosphorus  compensation 
were  doped  conventionally  since  this  would  introduce  an 
additional  impurity  fluctuation,  AP,  into  Fqs.  (34)  and 
(35).  These  curves  therefore  represent  a  best  possible 
situation.  For  a  P/B  ratio  of  approximately  two,  the 
final  NTD  resistivity  fluctuation  equals  the  starting 
material  fluctuation.  For  higher  Q  values,  the  fluctuations 
decrease.  Therefore,  for  best  uniformity,  NTD  compensated 
1R  material  should  have  a  Q  value  greater  than  two.  This 
result  is  in  agreement  with  conclusions  based  on 
detectivity  uniformity  measurements  vs.  Q  ratio  by  others. 

Figure  6  shows  the  decrease  in  100  Ap/p  for  the 
n-tvpc  side  of  intrinsic  (Q  =  1)  for  an  initial  resistivity 
of  18952  $2-cm.  Two  solid  curves  are  shown  representing 
initial  resistivity  by  fluctuations  of  1.23  and  7.12". 

These  values  bracket  the  initial  fluctions  of  the 
experimental  samples  reported  previously.' 


Reasonab 1 e 


'i<rure  6.  Comnar  i  son  of  resistivity  fluctuation  model 

with  experimental  data.  Circles  represent  NTD 
material  while  the  cross  represents  a  con¬ 
ventionally  compensated  wafer. 


agreement  between  experiment  and  theory  is  obtained 

therefore  justifying  the  assumptions  concerning  the 

starting  material  made  previously.  The  cross  shown 

in  Fig.  6  represents  a  conventionally  mc.l  t  -  compensated 

sample  resistivity  fluctuation.  It  is  clear  that  NT!) 

compensation  is  superior  in  agreement  with  the  model. 

C.  Maximum  Resistivity  Model 

* 

As  a  second  example  of  the  use  of  the  basic  N'TD 
doping  equations,  we  will  present  in  this  section  a 
model  which  calculates  the  maximum  mean  resistivity  which 
can  be  obtained  by  N'TD  compensation  before  the  initial 
resistivity  fluctuation  causes  mixed  typing,  i.e.,  part 
of  the  wafer  n-type  and  part  p-type.  These  results  arc 
important  for  other  detector  applications  which  rely 
upon  close  compensation  and  high  resistivity  for  proper 
operation  such  as  avalanche  visible  star  light  sensing 
detectors  or  x-ray  detectors. 

An  inspection  of  Fig.  3  shows  that  both  <  m.lx  and 
p  .  move  to  the  right  over  the  peak  at  pM  as  the 
fluencc,  and  hence  Q,  is  increased.  When  the  Q  value 
for  the  first  point,  p  ,  equals  unity,  then  p  =  i.  • 
and  any  further  compensation  will  produce  mixed  typing. 
This  mixed  typing  will  continue  until  sufficient 
irradiation  causes  the  Q  for  p  to  pass  Q  =  1  at 
which  point  the  mixed  typing  ceases  and  the  crystal  is 
entirely  n-type. 


IVe  can  easily  picture  from  Fig.  3  that  increasing 

the  initial  value  of  Ap/p  corresponding  to  the  shaded 

area  on  the  left  will  lower  the  average  resistivity 

obtainable  when  Pmnx  is  equal  to  pp.  We  have  calculated 

this  effect  and  show  the  results  in  Fig.  7.  If 

initially,  Ap/p  =  o,  then  the  maximum  value  of  p  is  pj. 

for  p-tvpe  and  p  ..  for  n-type  (see  Fig.  3).  •This  explains 

the  discontinuity  on  Fig.  7  at  Ap/p  =  o. 

The  procedure  to  calculate  the  rest  of  the  curves 

in  Fig.  7  is  slightly  different  for  n-type  and  p-type. 

For  p-typc,  p  and  p  .  are  easilv  obtained  from 
in  d  X  min 

p  and  Ap  using 


p 

max 


H  and  B  - „  can  now  be  calculated  from  p  r  ■  , 
max  min  max  (nun) 

using  Fqs.  (31)- (33).  Therefore,  R  and  aB  are  calculated 
as  before.  For  p-typc,  the  highest  value  of  p  occurs 
when  Pmax  =  p j | •  This  occurs  when  Q  for  pm  is  equal 
to  which  is  found  from  the  expression 

Qjl  =  1  -  f  (P,,2  -  n  j  “  i  /p1(B  ]  (39) 

where  pIf  =  fh  tu  from  Hq .  (8).  liquation  (39)  is  obtained 

by  solving  liq.  (S)  for  N  -  Np  and  tlien  substituting 
for  N ^  -  N'p  the  ({uantitv 

Na  -  Np  =  B  ( 1  -  Q)  . 


■»  *7 


It  is  clear  from  the  preceding  section  that 


min 


B  ( 1 


Q) 


and  that 

N\  -  Nd  =  B  ( 1  -  Q, j  +  AB/B). 

'max 

Note  that  we  arc  now  fixing  the  upper  resistivity  at  p.j 

and  allowing  the  fluctuation  aB  to  determine  the  lower 

resistivity  point  in  Pig.  3.  Prom  this  point  on,  the 

resistivities  are  calculated  from  -  N^  as  before 

and  once  p  .  and  p  arc  determined,  then  Ar/p  is 
mm  max  ’ 

easily  obtained.  The  procedure  is  identical  on  the 
n-type  side  except  that  now  the  higher  resistivity  is 
fixed  at  p.  at  Q  =  1  and  the  lower  resistivity  is 
determined  by  the  fluctuation  AB. 

An  inspection  of  the  results  of  Pig.  7  shows 
that  for  an  initial  resistivity  fluctuation  of  10", 
the  highest  p-type  average  resistivity  possible  is  just 
over  200,000  fl-cm  while  n-type  resistivities  above 
140,000  ft-cm  are  not  to  be  expected  without  mixed 
typing.  Claims  of  achieving  higher  resistivities  than 
these  by  Nil)  compensation’’  suggest  either  that  mixed 
typing  had  occurred  or  that  incomplete  annealing  of 
the  radiation  damage  contributed  to  the  resistivity. 

D .  Defect  Product  i  on _ Rate  Calculations 

In  this  section  t lie  calculation  of  the  rate  of 
displacements  by  fast  neutron  radiation  damage  is 


r 


extended  from  a  model  presented  previously.*  A  first 
principles  calculation  of  the  fast  neutron  lattice  damage 
involves  a  knowledge  of  five  factors: 

1.  The  fast  neutron  reactor  spectrum,  .,*>  ( I ^  )  . 

2.  The  spectrum  of  primary  knock-on  (PRO)  silicons 
produced  by  interaction  of  the  neutrons  and  lattice 
atoms . 

3.  The  fraction  of  I'KO  energy  delivered  to  the  lattice 
which  implies  a  calculation  of  ionization  energy  loss 
of  the  PKO’s. 

•1 .  The  average  energy  required  for  displacement  or 

vacancy- interst it ia  1  pair  production  (the  displace¬ 
ment  threshold). 

5.  The  disposition  of  vacancies  and  i nters i t i t a  1 s 

after  production  (defect  recombination  and  clustering). 
Previously,  we  had  estimated  the  neutron  energy 
spectrum  in  a  very  crude  way  in  terms  of  the  (Id  ratio, 
i.e.,  the  total  flux  below  and  above  the  (Id  absorption 
edge.  Also,  ionization  energy  losses  of  the  i'KO's  as 
they  come  to  rest  in  the  displacement  cascade  were 
ignored.  This  is  a  reasonable  assumption  for  heavy 
lattice  atoms  but  doubtful  in  the  ease  of  si  1  icon." 

These  factors  have  been  corrected  to  some  extent  in  t he 
present  model,  however,  no  attempt  to  calculate  item  .> 
will  be  made  since  this  is  extremely  complicated, 
hx per i men t a  1  evidence  will  he  presented  in  a  later 


3u 


section  which  suggests  that  about  90"  of  the  displacement 
defects  are  lost  at  room  temperature  through  annealing 
and  recomb  inat ion . 

We  shall  begin  by  briefly  reviewing  the  model 

presented  previously  since  many  of  the  calcul at iona 1 

steps  are  similar  and  will  serve  to  define  terms.  W'e 

6  - 7 

begin  with  a  well  known  result  from  conservation  ol 
energy  and  momentum  for  elastic  collisions  which  relates 
the  incident  neutron  energy,  li  to  the  silicon  atom 
recoil  energy,  !: , , 

I-.,  =  XH  sin"  -  (10) 

where  o  is  the  angle  through  which  the  neutron  is 
scattered  in  the  center- of-mass  coordinate  system. 

The  parameter  X  is  a  kinematic  constant  dependant  on 
the  masses  of  the  two  particles  and  is  given  by 

X  =  4Mj.M7/  CM  1  +  M3)2  (  1  1  ) 

where  Mj  is  the  neutron  mass  while  M-,  is  the  silicon 

atom  mass.  For  silicon,  x  =  4 (11 (28)/ Cl  +  28)“  =  O.losl. 

We  shall  attempt  to  keep  track  of  the  notation  used  in 

both.  Refs.  (>  and  ".  Therefore, 

h,  -  incident  neutron  energv  “  I. 

1  n 

li,  =  silicon  I’KO  recoil  energy  =  1  . 

Irom  l-.q.  (•Id)  it  is  clear  that  a  head-on  collision  !  *  ' 

will  v i e 1 d  the  maximum  possible  energy  transfer  to  the 


silicon  a  t  om 


i:2  =  Xlij  =  !;  n  (max )  (.12) 

(The  symbol  over  indicates  maximum  while  the 

inverted  symbol  "v"  will  indicate  minimum.)  In  any 

radiation  damage  problem,  there  will  be  a  minimum  recoil 

v 

energy  transferred  to  the  lattice  atom,  F,,,  below  which 
no  radiation  damage  will  occur  which  is  known  as  the 
displacement  threshold0  (item  4) 
v  v 

U2  =  xEl  =  Cd*  <•!•') 

v 

There  is  also  a  minimum  bombardment  neutron  energy,  F  , 

v 

which  can  produce  displacements  which  have  1'9  =  Lj  as 
their  maximum  possible  recoil  energy  as  given  by  l.q .  (4  3). 

To  investigate  displacements,  it  is  necessary  to 
know  how  the  neutron  interacts  with  the  lattice  atoms. 

This  is  expressed  in  the  differential  neutron  scattering 
cross  section,  do/dF-,.  For  a  moderated  reactor  neutron 
flux  spectrum,  very  few  neutrons  have  energies  above  the 
fission  peak  which  lias  a  maximum  at  about  1.5  McY.  he, 
therefore,  do  not  need  to  be  concerned  with  ladiation 
damage  caused  by  high  energy  nuclear  reactions  such  as 
(n,p),  (n,u)  and  (n,2n)  etc.  The  only  concerns  are 

therefore  thermal  neutron  capture  ( n , > )  and  recoil 
which  we  have  treated  previously'  and  clast  ic  scattering 
of  fast  neutrons  whose  damage  prnpert  ies  we  will  calculate 


here  . 


We  shall  assume  that  elastic  scattering  is  completely 


isotropic,  i.e.,  all  scattering  angles  are  equally 
6  ■*  7 

likely.  This  is  a  reasonable  assumption  for  the 

neutron  energies  of  interest  here,  however,  at  higher 

neutron  energies,  some  forward  scattering  is  known  to 

•  <-  6-7 
exist . 

Under  these  assumptions,  the  energy  dependence 
of  the  differential  cross  section  is  given  by^ 


°e/U:i 


o 


^2  —  ^'1 
f.2  >  Xl-j 


(4-n 


where  n  is  the  total  elastic  scattering  cross  section 
e 

-  2  4  7 

which  we  take  to  be  about  3  barns  (3  x  10  'em"). 

The  number  of  displacements  per  second,  dNL/dt, 
then  is  given,  for  a  monoenerget  i c  neutron  beam,  by(l 

dN^/dt  =  NjO  j<|>v  (43) 

2  2  7 

where  N^,  =  the  number  of  target  atoms  =  5  .\  10‘'“Si/cm'  , 

"> 

4>  is  the  monoenergic  neutron  flux  (n/cm'/sec)  ,  oj  is 

b  -  7 

the  displacement  cross-section  calculated  from  the 
interaction  cross-section  using  Hq.  (14) 


(  4(>  ) 
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and  v  Ls  known  as  the  damage  runction  '  and  represents 
the  mean  number  of  displacements  per  incident  neutron. 

In  general,  the  calculation  reduces  then  to  a 
calculation  of  v  which  depends  on  the  resultant  energy 
distribution  (item  2)  of  the  I'KO's  P(lij,l;?)  where 
P  C  Li  p  I:  1  d  F.  2  =  probability  that  a  P  KO  will  recoil  into 
dl:~  at  !•,  for  incident  neutron  energy  li.  . 

According  to  the  K inch  in- Pease  (k.-P.)  mode  I ^  , 

neglecting  any  energy  loss  of  the  PKO's  by  ionization, 

v  .  vr 

v(1't)  =  b?/21U  for  1.7  >•  li ,, 

therefore,  the  mean  damage  function  is  given  by1' 


v 


vUV!  PCM,  ,ii:)di:,  . 


1  4S 


for  our  simple  illustrative  case  of  a  monoenerg i c  neutron 
beam  and  hard-sphere  isotropic  scattering,  all  recoil 
energies  are  equally  likely/’  therefore,  p  is  a 
constant  and  since  it  is  a  probability  function 

v 

p  di-:,  =  p  f  r:2  -  i:  ,  j 

o  r 


1 


1 


1 


K,  -  i;7 


u:i  '  1;d 


Since,  for  any  reasonable  neutron  energy,  Ml 
we  have  P  -  =  I-,"1  or 


P(E,]  -  1/H, 


Therefore,  the  average  recoil  energy,  f:,,  is  given  hv 


E,  P(n0)dE, 


i  -  l]  ’ 

4-  Y.^/2 

i: ,  “  v 

-  !' 


Thus  . 

-  __  k  2  ^  j 

”  ~T~  =  ~2 


°r  Hp-  Hp,max/2  =  — —  in  the  notation  oi 


Re f .  1  . 


Also, 


vlE,)  P  dll. 


and,  using  Pqs.  (4  81  and  (491  , 


1  1 


since  li,/2  =  f^,  we  have  the  well  known  result 


f  /  ’  i  f  '  - 1 

27"1, 1  p'-’d 


Also  note  that 


it: 


/ 


di:. 


ITT  ' ' 2  =  °e 


v 


Therefore,  from  I:.q .  (45),  and  using  (42), 


d\ 


L) 


TF 


N, 


rv(£)  • 


(50  ) 


It  is  instructive  to  apply  this  monoenerget  i  c 
neutron  model  to  a  row-2  reactor  position  (see  next 
section)  to  determine  the  damage  rate  if  we  assume 
that  all  of  the  fast  neutrons  are  fission  neutrons 
with  an  energy  of  about  1  McV  and  determine  the  fast 
flux  from  the  known  thermal  flux  by  using  a  Cd  ratio  of 
about  10.  The  mean  thermal  flux  in  row- 2  is  about 
5  x  10  '  /cm  “/sec .  For  a  Cd- ratio  of  10,  the  fast 


,12 


flux  is,  therefore,  5  x  10  “n  j-/cm“/scc .  Using  hq .  (50) 
and  assuming  a  displacement  threshold  of  25  eV  and  a 

scattering  cross  section  ol  5  barns,  we  t  i nd  that 

'  F: 


IF 


(0.1  552)  (  1  x  1  0  ) 

1(25) 

1 552  d i spl accment s/neut  ron 


(511 


and 


5d 


\'TOe<j)V 


dNn 


=  (5  x  10“'"')  (.5  x  1  0  "  “  1 )  (3  x  1  0  1  2 )  (1332) 

=  5.09  x  10 14  disp/cm'Vsec  . 


We  shall  see  later  that  ionization  losses  of  the  I’KO's 
will  reduce  this  damage  rate  by  about  a  factor  of 
three.  The  above  model  was  presented  previously  for 
an  in-core  reactor  irradiation  and  represents  an  upper 
limit  to  the  radiation  damage  from  fast  neutrons 
which  can  be  expected. 

We  must  now  consider  the  moderated  neutron  energy 
spectrum  typical  of  the  graphite  reflector  region  of 
MURR.  To  a  good  approximation,  the  flux  of  a  moderated 
spectrum  follows  l/iij  up  to  a  fission  cut-off  energy 
F. and  is  zero  for  H ^  -  F(-.(l  Therefore, 

I  u/tpiwn  'i  i  l;t 


d<J)/dlij  = 


(5 


\’ 

where  L  =  i'n  (Xlij./li,)  and  the  differential  form  above 
represents  the  flux  spectrum  to  a  good  approximation. 
That  this  differential  form  is  a  good  representation  can 
be  seen  by  integration 


v 
F. , 


i: 


r 


f,A 


cl  i|) 


l'  mod 


fni:, 


=  t  i 

mod 

v 

Y.J  ■ 


is  the  total  number  of  neutrons  per  see  per 


where  d>  , 
mod 

cm"  which  can  produce  displacements  regardless  of  energy. 

The  calculation  of  the  damage  rate  is  similar 
to  the  above  except  that  the  distribution  of  1’KO's 
is  now  different  and  is  given  by^ 


This  equation  leads  to  Eq .  (30)  of  Ref.  1  as  follows. 

From  Ref.  6, 


v 


v(E2) 


dE 


2 


and  inserting  Eqs.  (44)  and  (47)  we  have 


dE. 


fp (max) 
2Ed 


(56) 


Then,  using  Eq.  (55)  and  the  definition  of  mean  recoil 
energy,  E^,  we  have 


t 

l 

E„ 


E2  P(E2)dE2 


E  -  F 
2  L2 


(V  -  e,“) 


L  -  1  2  A  E  p ( L  -  1) 


v 


But  E2  =  AEp  >>  E2  and  L  -  £n(E2/E2)  >>  1  since  lip  >  1  'V 
there  fore , 


AEp  ( A  E  p ) 


A  E , 


2  ( A E f ) L  2L 


or  in  terms  of  the  notation  used  previously  and  realizing 
that  E  (max)  =  F.p/x  ,  we  have 
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It 


Ep (max) /2 


£n  [E  (max) /Ej] 


as  in  Eq.  (30)  of  Ref.  1.  Combining  the  results  of 
Eqs.  (56)  and  (57),  we  find  that 
A  E  r 


4  Ed  Cn(AEf/Ed) 


(58) 


and 


dNr 


ur 


=  N„o  <k  ,v 
T  e  mod 


(59) 


Since  the  flux  spectrum  had  not  been  determined 
in  Ref.  1,  we  had  used  E^-  =  1.5  MeV  and  using  a  Cd-ratio 
of  10,  <*>mod  =  3  x  lO^n/cmVsec.  h'e  find  that  (we  had 
also  taken  E,  =  12  eV  but  will  use  25  eV  here) 


-  =  (.1332) (1 , 5  x  10°) 


4(25)  fn 
222.3  disp/neutron 


(.1332) (1 . 5  x  106)  ] 

25 


and 


dN 


dt 


-  =  (5  x  1  0 2 2 ) (3  x  10"24)  (3  x  1012)  (222.2) 


=  lx  1014  dis  p/cm /  sec 

which  is  about  a  factor  of  6  lower  than  the  monocnergct i c 
case  illustrated  previously. 

As  a  final  calculation,  we  shall  use  a  l/l:d  fit  to 
the  MURR  reactor  design  data  flux  group  calculations  for 
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row-2  of  the  reflector  and  will  also  take  ionization  losses 


of  the  PKO's  into  account  in  the  flavor  of  the  K-P's 
model.  The  best  1/P^  flux  spectrum  fit  ( liq .  (53)]  to 

the  design  data  has  been  determined  for  row- 2  to  be 
given  by  the  parameters 

K  =  ll>moc]/k  =  7.8  x  1011n/cm2/sec 

P.f  =  3.2  7  MeV. 


12  '■> 

The  above  parameters  then  give  $  .j  -  7.62  x  10  n/cmVsec 
and  L  =  f n (X E^/b^ ) *  9 . 7 7  giving  a  thermal  to  fast  ratio 
of  about  4  x  10^/7.62  x  10^“  =  5.25  in  reasonable 
agreement  with  a  Cd  ratio  of  8-10  (the  thermal  flux  is 
here  taken  as  the  peak  thermal  flux  in  row-2). 

To  account  for  electronic  energy  losses,  a 
characteristic  ionization  energy  threshold  for  the 
material,  P.,  is  estimated  and  then  the  K-P  model  is 

i 

modified  in  the  following  way,*'1 

v 

(  I  .  /  2  I  .  ,  Ii,  <  I  , 


v(t:,i  = 


V 


(  6  0  ) 


f.  '2b, 


i  - 


Since  the  stopping  power  is  some  function  of  the  number 
of  electrons  which  is  in  turn  a  function  of  the  lattice 
atomic  mass,  one  estimate  of  h.  is  obtained  from  the 
re  1  at  i  onsh  i  p(l 

]  .  I kev )  M  ,  ( a . m . u .  i  -  2  S  ke\  . 

A  second  estimate  for  sen;  iconduct  ors  and  :nsnlat"!> 
lias  been  g.  i  veil  by  So  i  f  r  as 


i  1 


r 


e.  = 

i 


NL  E 
_ £_ 


m  8 
e 


(1839)  (28)  4 

1  8 


=*  26  keV, 


in  good  agreement  with  the  above,  where  E ) '  is  the  mean 
energy  difference  between  the  middle  of  conduction  and 
valence  bands  which  we  take  to  be  about  4  eV.  The  effective 
threshold  for  ionization  could  be  as  much  as  a  factor  of 
three  higher,  however,  we  will  see  that  our  final  result 
is  relatively  insensitive  to  E . . 

We  now  proceed  as  before.  Equation  (55)  can  lie 
approximated  as 


P(E2)  -  4 


1 


1 


E 


E. 


( 6  3 ) 


since  L  >>  1,  and  E-,  =  AEj..  Then  using  E 
(6  5), 

.  n. 

E 


.  (48)  and 


v  = 


4 


'2  /  1 
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2Rd  \ R2 


dE. 


v 

R2  =  E 


'  1 

•  f  1 

,  \ 

dE  7 

E  ) 

E..  cl  ' 

i 

\  2 

‘2  / 

and  realizing  that  E .  >> 

we  find 

that 

J 


I 


1 

2E  ,L 


i.  -  ^ . 

E.  fn  +  -4- 
1  r;i  2E. 


Letting  B  =  E^/E?  =  E^/Xti^,  we  find  that 


[  p2  +  2 B  fn  Cl/3)  ] 
d 


4Ed  fn  (XEf/Ed) 


[&  +  2 B  fn  (1/6)] 


It  should  be  noted  that  this  result  is  identical  to 
Eq.  (58)  except  that  an  ionization  factor  [B"  +  2b  fn  1/6] 
has  been  added.  We  can  now  proceed  with  a  calculation  of 
the  damage  rate.  Using  the  parameters  which  describe 
our  flux  spectrum, 

Ef  =  3.27  MeV 

4>  ,  =  7.62  x  1 0  ^n/ cm ‘‘/sec  , 


E.  =  27  keV,  we  find  that 

l  ’ 

6  =  E./aEj.  =  6.  199  x  I0'2 

and  that  our  ionization  correction  term  is  u i ven  bv 


(  6 2  +  2  6  fn  1/6 )=  0.4067  , 


therefore,  from  Eq .  (64), 


»  =  0.406: 


(.1332) (3.27  x  10 


4(25)  fn 


(.  1  332)  (3.  2"  x  IP1’  )' 


0.4067  (446) 


181.1  d  i  sp/ iicu  t  ron 


1  3 


This  result  is  only  slightly  lower  than  our  previous 


model  where  w?c  had  taken  i;^-  =  1.5  MeV  and  ignored 

7 

ionization  losses.'  Calculating,  now,  the  damage 
rate,  we  find  that 

N't  ct  <p  ,  v 

I  e  mod 

(5  x  1  0  2  2 )  (5  x  1 0  ~  2  4  )  (7.62  x  1 0 1 2 ) ( 1 8  1  .  4  ) 
2.07  x  1014  d i sp/cm2/ sec . 

We  believe  this  last  damage  rate  to  he  the  most  accurate 
and  feel  that  the  number  can  not  be  improved  without  a 
detailed  measurement  of  the  neutron  spectrum  at  the 
irradiation  position  and  using  experimental  energy 
loss  data  of  silicon- ions  in  silicon.  These  results 
should  give  the  correct  damage  rate  to  within  a  factor 
of  2  or  3 . 

It  must  be  mentioned  that  the  damage  rate  in  the 
RSA1:  Bulk  Fool  variable  flux  facility  fnext  section)  is 
lower  for  the  following  reasons: 

1)  The  thermal  flux  f,  x  in  n/ent  /sec  or  about 
SO  times  lower  than  row- 2. 

2)  The  spectrum  is  somewhat  more  moderated  since 
the  Cd  ratio  is  about  AM. 

A  rough  estimate  of  the  displacement  rate  lor  RS\I  is 
therefore  1  x  lh'‘  to  2  x  1 0 4  “  d  i  sp/er.f' /sec  . 


i  t 


The  analytical  fit  to  the  reactor  design  data 

given  on  p.  -11  assumes  that  the  flux  per  unit  energy 

range  falls  as  1  / H  to  a  cutoff  energy  Hj-  and  is  then 

terminated.  The  four  flux  groups  and  their  energies 

used  to  determine  the  parameters  are: 

Troup  4  0-0.625  eV  4  x  10^n/cm“/sec 

Troup  3  0.625  ev-5.55  keV  0.027  1  0  '  "’n/cra"/  sec 

Troup  2  5.53  keV-821  keV  0.510  x  1 0 1 3n/cm2/sec 

Troup  I  821  keV-IO  MeV  0.1545  x  1 0  ] '\'n /ew2 / see 

Tro  iip  1  is  inconsequential  since  these  neutron 
energies  are  below  the  silicon  displacement  threshold. 

The  1,/li  spectrum  fit  yields  the  same  value  of  K  for 
groups  2  and  3.  The  value  of  1; ..  was  selected  for  group  1 
so  that  a  1/1;  spectrum  will  yield  the  total  flux  for  this 
group  when  integrated  over  the  group.  h'e  know  from 
spectrum  measurements  at  the  beryllium  reflector  that 
the  1/1:  s pcetrum  holds  for  groups  2  and  3,  however, 
the  flux  in  group  1  is  more  nearly  represented  by  the 
fission  spectrum  fit 

dt/dl.  =  k'  e  '  si  nh  yfTV. . 

V  in 

In  this  model,  the  parameter  k^  would  be  adjusted 
so  that  the  integral  of  this  function  over  group  1 
would  yield  the  total  flux.  Hn for t un a t e I v ,  other  cal¬ 
culations,  using  this  spectrum,  would  hn\o  t<,  K  p<  i 
formed  numerically.  ‘hnce  the  niirlu  t  .1 1  -p  I  a.  ervnt  • 
i  i ndependont  o r  cm  r".  \  above  !  (  ,  t  la  •  r  i.dod.ii  imi  ■ 
w  1 1  !  d  have  !  1  t  t  1  e  (  1  f  ■  t  <>n  1  ■  i j r  i  i  ui  1  t  . 


U’e  summarize  the  results  of  the  above  three  model'' 

in  Table  3.  Tt  is  clear  from  Table  3  that  the  order 

of  magnitude  of  the  displacement  rate  is  reasonably 

insensitive  to  the  nature  of  the  model  and  that  our 

earlier  estimates  bracket  nicely  the  more  refined 

result  of  the  last  model.  If  the  ionization  threshold, 

!i .  ,  is  increased  bv  a  factor  of  three,  then  the  ionization 
1 

correction  factor  becomes  O.b"  and  the  damage  rate  is 
almost  doubled.  Lowering  the  displacement  to  12  eV 
would  also  almost  double  the  damage  rate.  The  damage 
rate  is  also  sensitive  to  the  fission  cut-off,  L  . 

We  have  used  a  reasonable  number  for  this  parameter, 
however,  it  should  be  realized  that  some  neutron  flux 
exists  tip  to  10  McV ,  however  small,  and  these  high 
energy  neutrons  contribute  additional  displacements  from 
the  nuclear  reactor  by-products  they  produce. 

finally,  wc  emphasize  again  (see  Ref.  1)  that  the 
recoils  produced  from  neutron  capture  and  >  and  r 
emission  contribute  a  damage  rate  of  about  1  x  10  ‘\lisp/ 
cm' /sec  for  this  position  (or  about  3- 10*  of  the  fast 
neutron  damage).  These  damage  rates  will  he  compared  in 
a  later  section  to  the  total  number  of  defects  we  have 
ob>erved  experimentally.  It  will  he  seen  that  only 
3-10’, '  of  the  calculated  displacements  have  been  observed 
e x pe  r  i  men  t  a  11)'. 


It 


Table  3.  Mean  damage  function  and  damage  rate  for  row-2 

reflector  position  for  various  fast  neutron  models. 

Model  v  dVdt 

Uv  ‘  (displ a cement s/neut ron)  (d i spl a cement s/cm'V sec  ) 

Monoenerget  ic  beam  1332  6  x  1014 

E  =  1  MeV 
n 

f(fast)  =  3  x  10  * ~n/cm“/sec 

Cd  ratio  =  10 
No  ionization  correction 

1/En  moderated  spectrum  222.3  1  x  1 0 1 ^ 

Hf  =  1.5  MeV 

;mod  =  3  x  1 0 1 2n/cm2/sec 

Cd  ratio  =  10 
No  ionization  correction 

1/H  moderated  spectrum  181.4  2.07  x  10*'* 

to  fit  reactor  design 
dat  a 

Ef  =  3.2  7  MeV 

1  1  ? 

tmod  =  '-02  x  10  “n/cm“/scc 

4>  tbermal/<timo(.|  =  5.2  5 
!•:.  =  2  7  kcV 


III.  EXPERIMENTAL  TECHNIQUES  AND  APPARATUS  DEVELOPMENT 


In  this  section,  a  description  of  the  apparatus  and  tech¬ 
niques  required  to  irradiate  Si  for  NTD  and  to  evaluate  the 
material  before  and  after  irradiation  is  presented.  Only  new 
equipment  developed  or  used  under  this  contract  will  be  des¬ 
cribed  in  detail.  Previous  equipment  descriptions  are  to  be 
found  in  Ref.  1.  A  brief  review  of  irradiation  facilities  is 
presented,  however,  which  is  relevant  to  the  previous  section. 

A.  Reactor  Description 

The  reactor  at  the  University  of  Missouri  is  a  high 
flux,  light  water  moderated  research  reactor  designed  to 
operate  on  highly  enriched  fuel.  It  has  a  beryllium 
reflector,  a  graphite  moderator  and  forced  cooling  system 
for  pool  and  core.  The  reactor  is  operated  at  1(1  MW , 

120  hours  per  week,  52  weeks  per  year.  A  comparison  of  the 
peak  flux  of  the  MURR  reactor  with  other  research  reactors 
in  the  U.S.  is  shown  in  Eig.  8. 

A  cross-sectional  view  of  the  reactor  core  and  various 
experimental  stations  is  shown  in  fig.  0.  ('no  of  the 
major  advantages  of  MURR  for  Si  irradiation  is  the  easy 
access  to  sample  positions  in  the  pool  and  reflector. 

Presently,  two  2”'  din.  x  30"  long  silicon  sample 

irradiation  positions  with  sample  rotators  are  located 

in  the  graphite  reflector  in  positions  1-2  and  l!-f  (  rmv-  .  i 

in  I  ig.  P .  I  he  average  flux  at  these  positions  is  about 
1  5 

5  x  Id  n/cm  /sec.  These  an  the  positions  !»i  which  tl< 
damige  rate  was  calculated  in  the  previous  section.  i  I  u 
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^ure  S.  Comparison  of  research  reactor  peak  fluxes 
for  U.  S .  reactors. 
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!  idlin'  P . 


i  lux  t  r  :i  p  ,  he  ampo  r  t  s 
fac i 1 i t i es . 


reflector  and  hulk  pool 


i  r  u  i  e 


l1'1.  1  l,er",:i  1  nciil  ron  flux  profiles  in  1  -  f  reflect 

pos it  i '  >  n . 


pneumatic  tubes  used  for  neutron  activation  analysis 
(NAA)  are  also  shown.  Figure  10  shows  the  vertical 
thermal  neutron  flux  profile  for  the  second  row.  This 
profile  should  be  compared  to  the  vertical  profile  which 
we  have  determined  for  the  RSAF  bulk  pool  variable  flux 
facility  (later  section). 

In  addition  to  the  two  2"  dia.  irradiation  positions 
in  row- 2,  two  3"  dia.  x  30”  irradiation  positions  occupy 
row-3.  In  addition,  3"  dia.  material  can  be  irradiated 
in  the  seven  3”  dia.  positions  of  the  S-basket  (Fig.  9). 
For  the  very  low  fluences  required  for  exact  compensation 
of  detector  Si,  a  variable  flux  bulk  pool  facility  has 
been  constructed.  This  facility,  also  on  Fig.  9,  is 
known  as  RSAF.  The  peak  thermal  flux  in  this  facility 
is  about  5x10  n/cm  /sec  in  its  present  position. 

All  of  the  Si  irradiation  positions  have  self- 

powered  neutron  detectors  whose  output  is  integrated 

bv  analog  current  integrators.  Reproducibility  of  the 

1  _  ? 

neutron  fluence  is  generally  better  than  l'. 

H .  EPR  Spectrometer 

Microscopic  defect  structure  is  determined  by  elec¬ 
tron  paramagnetic  resonance  (FPR)  or,  as  it  is  more  com¬ 
monly  known,  electron  spin  resonance  (F.SR).  The  FSR 
instrument  used  for  these  experiments  is  a  conventional 
X-band  (9  (illz)  spectrometer  with  Schot  t  k  y -ba  r  r  i  e  r  diode 
detectors,  and  a  12  inch  Varian  mac. no  l  with  field  control. 


e  a 


The  spectrometer  can  be  operated  either  as  a  homo- 
dyne  system  with  100  kllz  field  modulation,  or  in  a  super¬ 
heterodyne  mode  with  30  Mllz  l.F. 

Accessories  include  liquid  helium  and  nitrogen 
Dewars,  a  variable  temperature  fgas  flow)  cavity,  and 
cavities  with  optical  access.  An  HNDOR  facility  is  also 
avail ab 1 e . 

C .  Liquid  Nitrogen  Variable  Temperature  Dewar 

A  variable  temperature  liquid  nitrogen  Dewar  has 
been  constructed  under  the  present  contract.  It  is  shown 
schematically  in  Fig.  11.  The  Dewar  is  constructed 
entirely  of  stainless-steel  except  for  the  sample  holder 
and  cold  finger  tip  which  are  copper.  The  liquid  nitro¬ 
gen  volume  is  approximately  1.2  l.  The  holding  time  at 
77°K  is  over  24  hrs.  due  to  t  ire  long  narrow  neck  design. 
Optical  window  flanges  with  NaCl  windows  are  provided 
so  that  optical  absorption  and  photoconductivity  can 
be  performed  (only  the  entrance  window  is  shown  in  Fig. 
11).  The  spacing  between  windows  is  sufficiently  narrow 
so  that  this  Dewar  will  fit  all  optical  spectrometers 
now  at  the  reactor  and  physics  department.  The  windows 
are  sufficiently  large'  so  as  to  accept  an  f/S  cone  of 
light  when  close  to  the  spectrometer  entrance  slits. 

F.lectrical  feedthroughs  are  provided  at  the  top  of 
the  Dewar  so  that  temperature  sensor  leads  and  sample 
leads  may  be  brought  out.  it  was  planned  to  use  this 


Dewar  for  DLTS  as  well  as  optical  measurements.  There¬ 
fore,  provisions  have  been  made  to  isolate  the  cold  tip 
by  a  valve  and  stem  system.  A  heater  wire  is  wrapped  on 
the  isolated  cold  tip  while  residual  nitrogen  is  vented 
through  the  stem.  It  is  then  possible  to  provide 
a  temperature  ramp  by  using  a  programmable  power  supply 
and  microprocessor. 

1).  Capacitance  Bridge,  Junction  Capacitance  and 

Dopant  Concentration 

A  Boonton  capacitance  meter  has  been  purchased  and 
has  been  used  for  both  DLTS  and  junction  capacitance  mea¬ 
surements.  Th is  meter  has  been  checked  for  accuracy  by 
measuring  the  junction  capacitance  of  several  commercial 
diodes  as  a  function  of  reverse  bias  voltage  on  this  meter 
and  also  on  a  1  MHz  bridge  for  comparison.  The  balancing 
bridge  can  compensate  for  shunt  conductance,  however, 
the  Roonton  meter  does  not  in  principle. 

Table  4  lists  values  of  C  vs.  V  for  a  small  power 
diode  as  measured  on  our  Boonton  meter  and  a  1  Mllz  bridge 
which  compensates  for  conductance.  It  can  be  seen  that 
the  agreement  is  usually  better  than  1  °0  between  the  two 
instruments,  Several  shunt  resistances  were  placed  in 
parallel  with  the  diode  and  the  capacitance  at  zero  bias 
observed  using  tiic  Boonton.  lor  a  shunt  resistance  of 
10  Ks’,  the  changes  in  capacitance  were  negligible;  for 
1  Ko,  the  change  was  about  a";  however,  no  readings  could 


lie  obtained  for  a  shunt  of  100  i  .  It  therefore  appears 


Table  4.  Comparison  of  reverse  bias  junction  capacitance 
taken  on  1  MHz  bridge  and  Boonton  capacitance 
meter . 


C(pf) 

Boonton 

C.(pf) 

Br  i  dge 

VR  (Volts) 

218 

217.4 

0 

165.4 

16  3.5 

0. 5012 

141.2 

138.5 

1 

115.5 

115 

85.1 

84 

5 

65.4 

65.2 

10 

49.2 

20 

41.4 

30 

36.5 

40 

83.0 

50 

that  the  Boonton  is  insensitive  to  any  shunt  leakage 
conductance  which  might  be  encountered  in  any  practical 
diode  for  small  bias  voltages. 

Junction  capacitance  can  be  used  to  determine  the 
concentration  of  dopant  in  the  lightly  doped’  side  ol  an 
abrupt  junction.  It  can  easily  be  shown  that  the 
junction  capacitance  for  a  reverse  biased  abrupt  junc¬ 
tion  (and  also  a  Schott kv  diode)  of  bias  voltage  V 
is  given  by 


+  V  -  kT)]*2 


C 


a  t  q 


NR/2(Vbi 


(65) 


where  A  is  the  junction  area,  q  is  the  charge  on  the  elec¬ 
tron  in  Coulombs,  NR  the  dopant  concentration  in  the  low 
concentration  side  of  the  junction,  .  the  built-in 
junction  voltage  at  zero  bias  and  =  <r  the  dielec¬ 
tric  constant  of  the  diode  material.  If  is  given  in 
f/cm,  then  C  is  in  farads.  This  equation  can  be  rewritten 
(neglecting  kT)  as 


1/C2 


q  es  NbA 


(06) 


so  a  plot  of  1/C2  vs  V  should  produce  a  straight  line. 
Note  that 


d(l/C2) 

dV 


es  V' 


(on 


Therefore,  the  slope  of  a  1 / C “  vs.  V  plot  determines  the 
dopant  concentration,  NR,  if  the  junction  area,  A,  is 
known . 

When  V  =  0, 


->  2V,  . 

1/C2  =  - — — ^ 

q  es  NrA- 

and  using  this  and  bq.  (66) 

n/c2iv.0 

V.  .  = - • 

Dl  d(l/C")/dV 


(68  ) 


The  built-  in  vo  1 1 age ,  j 

V  using  this  equation. 


,  can  be  determined  from  C  vs. 


The 

W  = 


depletion 


C 


width, 


is  given  by 


(69  ) 


so  that  use  of  Tqs.  (67)  and  (69)  will  yield  the  dopant 
profile  as  a  function  of  distance  from  the  junction  on 
the  low  dopant  concentration  side.  It  should  also  be 
noted  that  I'q.  (67)  is  equivalent  to 

n  = - - - -  ’  ("in 

q  cs  A"  (dC/dV) 

an  often  used  relationship. 

To  gain  experience  with  these  ideas,  several  types 
of  diodes  were  selected  at  random  and  1/C  vs.  Y  measured. 
The  results  arc  shown  in  Tig.  Id.  All  diodes  exhibit  a 
linear  relation  as  expected  up  to  10  volts.  Tor  higher 
bias,  data  points  arc  generally  lower  than  expected. 

The  diode  shown  in  Tig.  13  was  analyzed  for  dopant 
concentration  in  the  base  using  lie).  (6").  From  the  ex¬ 
terior  physical  size  of  the1  diode,  it  was  estimated  that 

"> 

the  junction  area  was  of  the  order  of  0.1  cm" .  Using  this 
area,  was  found  to  be  4.43  x  1 0 ' ' cm  '.  Typical 

dopant  levels  for  high  breakdown  voltage  diodes  are 

11 

usually  in  the  low  10  range.  Ter* a  inly  a  factor  o I 
two  error  in  est  imat  ing  the  iunct  ion  area  was  possible. 


Using  Iiq.  (68),  the  built-in  voltage  was  found  to 


be  0.778  volts.  An  estimate  of  N ^  can  also  be  obtained 

ij 

from  using  the  lol  lowing  relations 

q  V,  .  =  f  -  ci  V  -  q  V 
1  b  1  g  1  p  1  n 

*  Eo  -  kT  In  (Nv/Na)  -  kT  In  (Nc/N])) 

where  1-  is  the  I)  and  gap ,  N(.  and  Ny  the  density  of 
states  in  conduction  and  valence  bands  and  N .  and  Nn 
the  donor  and  acceptor  concentration  on  cither  side  of 
the  junction. 

Assuming  that  the  base  (low)  region  is  n-type,  as 
is  typical  for  power  diodes,  and  that  the  junction  is 
formed  by  diffusing  a  p+  region  into  the  n-type  wafers, 
then  the  Fermi  level  will  be  at  the  band  edge  in  the 
p-region  and  the  last  equation  can  be  approximated  by 

qVb  i  ~  Eg  -  kT  ln  (Nc/Nj)) 

where  =  \'  is  the  base  dopant  concentration.  I  h  i  s 
equation  yields  4.7  x  lO^cm  in  good  agreement  with 
our  previous  estimate  of  4.4  x  1 0 '  'em  \  it  should  be 
noted  that  while  this  second  technique  is  not  as  accurate, 
it  has  the  virtue  of  being  independent  of  junction  area, 
an  unknown  for  a  commercial  diode. 

The  depletion  widths  as  determined  using  Fq .  (09) 

are  shown  on  i ig.  la  for  various  bias  voltages.  The 
lower  slope  in  Fig.  la  at  higher  I'ias  voltages  corres- 

f.n 


(farads 


!■' if.urc  1  7> . 


ponds  to  a  dopant  concentration  of  7.2  x  10* '''em  which 

14-3 

is  close  to  typical  values  of  10  era  ' .  It  theretore 
appears  that  some  of  the  p+  dopant  has  retrodoped  into 
the  base  region  of  this  diode  over  a  distance  of  about 
0.5  um.  This  mixing  of  dopant  regions  is  similar  to 
data  taken  from  Sze  which  is  shown  in  Fig.  14  for  an 
n+  on  n  region.  This  retrograde  doping  also  extends 
over  a  region  of  about  5  um.  We  conclude,  therefore, 
that  the  diode  shown  in  Fig.  13  deviates  from  a  straight 
line  because  of  retrograde  doping. 

H .  Furnace  Tube  Cleaning  Procedure 

When  compensating  silicon  to  high  resistivities  by 
NTD,  very  small  concentrations  of  trace  impurities  can 
be  introduced  into  the  sample  by  diffusion.  Sources  of 
contamination  arc  etchants,  water  supply,  furnace  tube, 
protective  silicon  spacers  (see  Ref,  1  for  anneal  pro¬ 
cedure)  and  the  sample  boat.  Typical  care  in  sample 

handling  usually  introduces  concentrations  of  the  order 
12  -  5 

of  10  cm  ' ,  however,  if  the  anneal  temperature  is  low 
enough,  these  diffusions  can  usually  be  removed  by  a 
second  lapping  after  the  anneal. 

Unfortunately,  our  furnace  tube  apparently  picked 
up  a  contamination  source  while  annealing  a  hatch  of  test 
wafers  for  Rockwell.  We  believe  the  source  of  contamina¬ 
tion  to  have  come  from  the  epoxy  used  to  hold  the  ingot 
during  slicing.  Many  of  t  lie  wafers  in  this  batch  still 
contained  traces  of  this  epoxy  which  we  tried  to  remove 


before  irradiation.  After  annealing  these  samples,  it 
was  found  that  all  our  NTD  resistivities  were  at  first 


j 

j 

1 

too  high  by  orders  of  magnitude  indicating  the  dominance  j 

of  acceptor  impurities.  After  baking  the  furnace  tube 

for  several  weeks  with  flowing  argon  circulating  through  j 

the  tube,  it  was  found  that  the  resistivities  of  annealed  j 

samples  all  approached  5  fl-cm  n-type  indicating  that 

donors  now  dominated  the  impurity  contamination.  After 

trying  to  bake  the  furnace  tube  for  extended  periods 

with  little  improvement,  it  was  decided  to  clean  the 

tube.  The  cleaning  procedure  for  the  spectrosil  quartz 

furnace  tube,  quartz  boat,  and  Si  spacers  was  as  follows: 

QUARTZ  TUBE  CLEANING  PROCEDURE 

1 .  Rinse  in  Xylene 

2 .  Rinse  in  U . I .  -  H ? 0 

3.  Rinse  in  Aquaregia  (  3HC1/1  UNO,) 

4  .  R  inse  in  D .  I  .  -  1I70 

5.  Clean  with  Microsol  detergent 
b.  Rinse  in  1).  1  .  -  11^0 

BOAT  AND  SPACF.R  CM- AN  INC  fPrc-Cleaning) 

1.  Soak  in  Acetone  24  hours,  rinse  in  U.l.  -  11,0 

2.  Soak  in  lit  by  1  Alcohol,  "2  hours,  rinse  in  U.T.  -  11,0 

3.  lit ch  in  HNO^  for  5  min.,  rinse  in  U.l.  -  11,0 

4.  litch  with  commercial  CP-4  for  S  min.,  rinse  in 

U.l.  -  11,0 


o  ! 


I 


BOAT  AND  SPACER  STANDARD  CLEANING 
BEFORE  MACH  ANNEAL 

1.  Etch  in  Caro*  etch  5  min.,  rinse  in  D.l.  -  11,0 

2.  Etch  in  10.'.  Ill  3  mm.,  rinse  in  II.  I.  -  11,0 

5.  Etch  in  1  11,0:1  11,0,:  I  HOI  5  min.,  rinse  in 

D.l.  -  H,0 

4.  Bake  for  48  hrs.  in  flowing  argon  at  ]000°C. 

Various  annealing  tests  were  performed  on  high 
resistivity  wafers  #1  -  #8  to  test  the  cleaning  proce¬ 
dures.  The  anneal  experiments  “  1-8  have  the  following 
histories : 

Before  test  #1: 

1.  Tube  cleaned  as  above 

2.  Tube  baked  at  900°C  in  argon  for  4b. 5  hours 

3.  Sample  etched  in  commercial  CP-4 

4.  Sample  annealed  at  850°C  for  15  min.  without 
spacers  and  boat. 

Before  test  # 2 : 

1.  Tube  baked  at  900°C  in  argon  for  48  hours 

2.  Sample  etched  in  lab-mixed  CP-4 

3.  Sample  annealed  at  850°C  for  15  min.  without 
spacers  and  boat. 

Before  test  # 3 : 

1.  Tube  baked  at  900°C  in  argon  for  24  hours 

2.  Sample  etched  in  lab  CP- 4 

3.  Sample  annealed  in  air  at  850°C  for  15  min. 
without  spacers  and  boat. 


Before  test  #  4  : 


1.  Baked  tube  at  950°C  for  24  hours 

2.  Sample  etched  in  lab  CP-4 

5.  Sample  annealed  at  850°C  argon  at  high  flow  rate, 
no  boat  or  spacer  was  used. 

Before  tests  #5,  >*f  6  and  #7: 

1.  Bake  tube  at  1000°C  for  24  to  72  hours  in  argon 

2.  Etch  sample  in  lab  CP- 4 

5.  Sample  annealed  in  argon  at  850°C  with  boats  and 
spacers . 

Before  test  "8: 

1.  Bake  furnace  tube  for  520  hours  at  1000°C  (boat 
and  spacers  showed  a  greenish  tint  from  some 
unknown  source  of  con  tarn i nat i on - -perhaps  argon 
gas)  . 

2.  Ctch  sample  in  commercial  CP-4 

3.  Sample  annealed  in  argon  at  S50°C  for  15  min. 
with  boat  and  spacers. 

The  resultant  resistivities  for  these  experiments 
are  shown  in  Table  5.  It  is  clear  that  extended  baking 
in  flowing  argon  tends  to  help,  however,  the  best  result 
were  obtained  when  the  wafers  were  protected  in  tire  boat 
with  spacers.  Additional  long  term  baking  for  520  hours 
seems  to  have  caused  an  accumulation  of  contam i na t i on 
on  the  spacers  and  boat.  After  several  months  use  after 
tiie  final  cleaning  above,  N'l'P  resistivities  in  the 
100,  nun  r  cm  range  were  again  possible. 


Table  5.  Listing  of  the  initial,  final  annealed,  and 
annealed  lapped  resistivities  for  the  test 
wafers  used  to  develop  a  furnace  cleaning 
procedure . 


Anneal ing 

Run 

p  ^ (fi-cm) 

p  ^ (fi- cm) 
annealed 

p  ^ (fi-cm) 
lapped 

n 

6255 (P) 

779  (N) 

596 (N) 

H  2 

4902  (P) 

1465  (N) 

10690 (N) 
35250 (N) 

n 

5  30  B (P) 

948  (N) 

4254 (N) 

ft  4 

5712 (P) 

3923 (N) 

17422 (N) 

tt  5 

5849 (P) 

8647 (P) 

— 

tt  6 

5334 (P) 

8305 (P) 

8616(F) 

n 

5918  CP) 

11220 (PI 

10185 (P) 

tt  8 

7144 (P) 

3596  (N) 

20265(F) 

F .  1R  Optical  System  Calibration 

A  considerable  effort  has  been  spent  to  bring  the 
Perkin-Elmer  Model  112  double-pass  IR-prism  spectrometer 
to  best  working  order  so  as  to  obtain  maximum  signal-to- 
noise  and  resolution.  Although  this  is  the  principal 
instrument  used  to  date,  a  Beckman  IR-9  (2-25ym)  has 
been  placed  in  working  order  in  the  physics  department 
and  this  instrument  along  with  the  Cary  offer  double 
beam  capability  from  the  Si  hand  edge  out  to  25  ym.  In 


I'dition,  a  Spcx-0.5  rn  instrument  with  gratings  from 


lym-50ym  has  been  purchased,  however,  it  was  not  used 
under  the  present  contract  because  of  delivery  and 
set-up  time.  A  double  prating  lm  Spex  and  tunable 
laser  Raman  facility  is  also  set-up  and  operating 
in  Physics  and  experimental  samples  were  irradiated 
to  investigate  the  possibility  of  studying  radiation 
damage  by  this  technique.  The  balance  of  this  sub¬ 
section  will  describe  the  placing  in  operation  of  the 
Perkin-Elmer  112. 

The  112  system  has  a  variety  of  prisms  (Glass, 

Quartz,  CaF-,,  NaCl ,  KBr  and  KRS-5)  making  it  useful  as 
a  light  source  or  absorption  spectrometer  from  the 
visible  to  SOym.  The  NaCl  and  KBr  prisms  were  badly 
damaged  by  water  vapor  and  have  been  refigured.  The 
thermocouple  detector  lias  been  replace  1  with  a  new  detec¬ 
tor  which  has  a  KRS-5  window  (l-50ym). 

The  glass,  Cal'7,  NaCl  and  KBr  prisms  have  been  pro¬ 
perly  aligned  and  the  system  focus  fixed  for  each  of 
these  prisms.  The  spectrometer  chopper  contacts  have 
been  replaced  (they  are  Dclco  auto  points  of  1950  vintage). 
The  electronics  have  been  maintained  and  are  functioning 
proper ly . 

Since  the  prisms  have  a  non-linear  dispersion,  the 
drum-number  of  the  spectrometer  must  lie  calibrated  vs. 
wavelength  using  known  absorption  lines.  This  has  been 
done  using  over  45  lines  from  llg-emi ssion ,  C07  absorption, 
and  polystyrene  film  (50um). 
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drum  number  using 


Typical  glowbar  intensity  spectra  taken  with  the 
NaCl  and  CaF?  prisms  are  shown  in  Figs.  15  and  16. 
Wavelengths  marked  on  the  various  air  absorption  lines 
are  in  pm.  It  can  be  seen  from  these  figures  that 
ample  calibration  lines  arc  available  out  to  about 
Spin  so  the  CaF0  prism  is  calibrated  over  most  of  its 
useful  range.  The  calibration  of  the  NaCl  prism  is 
not  as  well  known  between  10-25pm  due  to  a  lack  of  sharp 
absorption  lines  in  this  region,  and  the  calibration  is 
based  mainly  on  polystyrene  absorption  lines.  The  reso¬ 
lutions  and  s  igna  1  -  to-noise  shown  in  Figs.  15  and  16  arc- 
typical  of  data  shown  in  the  Perkin-Flmer  manual.  The 
resolution  has  been  checked  at  selected  wavelengths  and 
is  nearly  theoretical  as  discussed  later.  A  drying 
column  and  flushing  system  has  been  developed  which 
effectively  eliminates  most  of  the  air  absorption. 

A  program  has  been  developed  to  find  the  best  poly¬ 
nomial  least  squares  fit  for  drum  number,  N,  vs.  wave¬ 
length,  X.  A  five  parameter  fit  using  the  equation 


N  '  A  *  ?  *  7  *  r*  *  t*  *  r* 

was  tried.  The  estimated  error  between  fit  and  date 
points  is  generally  less  than  15. 

An  expression  for  the  dispersion  of  Cal  ,,  due  ori¬ 
ginally  to  Ilcrzbergcr  *  ^ ,  has  been  used  to  estimate  re 
solving  power.  This  equation. 


n  =  A  +  BL  +  CL2  +  DL3  +  EL4, 

7  _  1 

where  L  =  (X"  -  0.028)  with  X  given  in  pm,  has  been 
used  to  determine  n(X),  the  index  of  refraction,  and 
dn/dx,  the  linear  dispersion.  Since  the  resolving  power 
is  given  by 

P,  =  4 B  dn/dx 

where  R  =  8  cm  is  the  prism  base  width,  it  is  possible 
to  calculate  the  geometrical  image  width  AX  in  pm  vs. 
slit  width,  w,  from  the  expression  from  Harrison's  spec¬ 
troscopy  text 4  ^ 

AX  b 

w  ”  r  j  p  ? 

In  this  expression,  b  is  the  exit  prism  beam  width  (•  5  cm) 

and  F  is  t he  col  1 imator  focal  length  ('>■  27  cm).  IV  c  can 

use  these  expressions  to  calculate  AX/w  and  compare  these 

with  experimental  data.  The  results  are  shown  in  Table  (i . 

Table  0 .  Comparison  of  Hg  emission  line  half-widths 
per  unit  slitwidth  vs.  calculated  values. 

Ra  t  i  o 

X _ AX/w  (calc.)  AX/w  (exp.) _ 1  exp/ca  1  c  ) 

0 . 5 40pm  161.2  A  7mm  200  A  7mm  1.24 

1.014pm  8  7  8.1  5  A  nun  1000  A  7  mm  1.14 

Since  the  agreement  between  calculated  and  measured 
resolving  powers  is  good,  we  conclude  that  we  have 
achieved  reasonable  optical  alignment  of  the  instrument 
for  the  Cal.,  prism.  Once  this  is  done,  the  other  prisms 
need  only  be  focused,  a  single  adjustment,  which  is  an 


7  ~> 


FWHM  (A; 


Slit  Width  (mm) 


Figure  17. 


Full  width  at  half  max.  of  1.014  uni  llg  line 
vs.  spectrometer  slit  width  using  Cal prism 


extremely  accurate  procedure  since  it  relys  on  the  knife- 
edge  test  which  determines  focus  to  fractions  of  a  wave¬ 
length  . 

We  have  also  checked  linearity  of  slit  opening  by 
finding  the  approximate  area  under  the  1.014pm  Hg  emis¬ 
sion  line  as  a  function  of  slit  width.  The  (peak 
1, 

height  x  FW1W) 2  should  be  proportional  to  slit  width. 

The  experimental  data  for  FIVHM  vs.  slit  width  is  shown 
in  Fig.  17  while  the  (area) 2  vs.  slit  width  is  shown  in 
Fig.  18.  It  can  be  seen  that  the  linearity  of  slit 
opening  is  good,  however,  there  is  a  zero  offset  of 
approximately  0.04  mm  before  the  slit  opens.  This  is  of 
little  practical  concern  except  when  checking  instrument 
resolution  so  we  have  not  adjusted  it. 

A  calibration  of  glowbar  power  vs.  radiant  enittance 
and  blackbody  temperature  has  also  been  determined. 

This  information  will  be  of  use  in  determining  absolute 
light  intensity  at  the  detector  and  in  designing  and 
building  additional  blackbody  sources  and  glowbar 
sources  for  other  spectrometers.  The  blackbody  tempera¬ 
ture  as  a  function  of  electrical  power  into  the  glowbar 
was  determined  from  X  ,  the  blackbodv  peak  wavelength 
(See  Figs.  15  and  16)  for  different  settings  of  t  lie 
glowbar  power  supply.  The  power  to  the  glowbar  is 
measured  with  a  conventional  panel  watt -meter.  The 
blackbody  temperature  can  be  calculated  from 


Radiant  Emittance  at  B.B.  Peak  (W/cm^/pm) 
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ip,ure  19.  Radiant  omittance  at  blackbody  peak  vs, 

ylowbar  electrical  power.  The  photon  power 
H'^Xmax  *s  plotted  in  watts/cm-  on  t  lie  saline 
scale  as  t  lie*  omittance  which  is  in  K/cnr/i.r 


A  T  =  2897  pm- deg 
max 

From  the  radiant  omittance  at  a  ,  K  is  given  by 

max  Amax 

U'  =  1.3  x  10  watts/cm“-pm 

max 

Using  this  equation  and  experimental  data  for  A  at 

different  power  settings,  it  is  then  possible  to  plot 

IV  vs.  glowbar  power.  A  plot  of  this  data  is  shown 

max 

in  Fig.  19.  Note  that  for  a  glowbar  power  of  100  watts 

the  optical  power  emitted  is  about  10  watts/nif ,  Since 

1 

the  surface  area  of  the  glowbar  rod  is  about  10  cm", 
the  electrical  power  is  being  converted  into  IK  optical 
power  at  very  high  efficiency  suggesting  that  the  source 
is  well  designed  and  constructed. 

Figure  20  shows  the  blackbody  temperature  deter¬ 
mined  experimentally  from  A  x  vs.  glowbar  power.  The 
solid  line  is  theoretically  determined  from  the  rela¬ 
tionships  shown  on  the  figure.  It  can  be  seen  that 
only  a  small  gain  in  blackbody  temperature  can  be 
expected  by  further  increasing  the  glowbar  power  beyond 
300  watts.  If  the  power  is  doubled  to  000  watts,  the 
temperature  increases  from  — lb?S°K  to  ~2000°K,  however, 
this  will  produce  very  little  increase  in  optical  power 
at  10  pm.  he  believe,  therefore,  that  the  present 
glowbar  source  is  about  optimum. 

In  Fig.  21,  the  I’erk  i  n  -  F  1  me  r  chopper  amplifier  gain 
was  determined  using  a  10  pV  t es t  input  signal  and 
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by  measuring  amplifier  output  for  various  settings  of 
the  gain  control.  It  can  be  seen  that  the  gain  switch 
changes  the  gain  by  about  4  db/step.  The  maximum  gain 
obtainable  is  about  4  x  10^.  Since  the  detector  sen¬ 
sitivity  is  about  10  pV/pW,  then  1  W  of  optical  power 
on  the  detector  will  produce  an  electrical  signal  output 
of  S  volts.  Since  optical  signals  might  be  expected  to 
be  in  the  nW  range,  output  will  be  in  the  mV  range  when 
the  full  gain  of  5  x  10^  is  used.  I'.’e  conclude,  there¬ 
fore,  that  the  amplifier  has  sufficient  gain  to  perform 
adequately . 

Finally,  we  mention  that  a  microfloppy  disk  has 
been  assembled  to  store  data  and  programs  useful  to  this 
contract.  This  disk  system  interfaces  with  a  variety 
of  6800  microprocessor  systems  on  hand  and  has  been  used 
to  store  analysis  programs  for  both  DLTS  and  for  the 
theories  presented  in  an  earlier  section. 

G.  DLTS  System 

Deep  levels  in  NTD-silicon  have  been  investigated 
using  variations  of  standard  depletion  layer  junction 
transient  techniques.  The  circuit  diagram  for  a  junc¬ 
tion  capacitance  DLTS  apparatus  is  shown  in  Fig.  22. 

All  timing  is  derived  from  the  reference  oscillator 
at  the  top  of  the  figure.  This  oscillator  provides 
a  reference  signal  for  the  lock-in  amplifier  detec¬ 
tor  as  well  as  providing  the  trigger  inputs  for  the 
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pulse  generator  for  diode  bias  and  the  sample  and 
hold  circuit. 

Bias  voltage,  which  forms  a  diode  depletion  re¬ 
gion  and  which  momentarily  collapses  this  depletion 
region,  is  obtained  from  a  pulse  generator  and  varia¬ 
ble  D.C.  power  supply  in  series.  This  arrangement 
will  momentarily  switch  the  reverse  bias,  Vp ,  up  near 
zero  volts  to  collapse  the  depicted  region  and  fill 
the  traps.  An  output  voltage  proportional  to  the 
junction  capacitance  is  fed  into  the  lock-in  signal 
input  through  a  sample  and  hold  circuit.  The  purpose 
of  this  sample  and  hold  is  to  remove  the  large  capac¬ 
itance  pulse  which  occurs  during  the  0.1  psec  bias 
pulse  when  the  depletion  region  is  collapsed.  With¬ 
out  this  circuit,  the  capacitance  pulse  will  easily 
overload  the  input  to  the  lock-in.  The  lock-in  out¬ 
put  is  fed  to  the  Y-channel  input  of  the  X-Y  recorder 
while  the  voltage  from  a  platinum  resistance  thermo¬ 
meter  is  fed  to  the  X-channel  input.- 

If,  instead  of  junction  capacitance,  junction 
current  is  measured,  then  the  experiment  is  generally 
known  as  "current  DLTS"  or  "transient  current  spec¬ 
troscopy"  (TSC) .  The  circuit  diagram  for  this  experi¬ 
ment  is  shown  in  Fig.  25.  This  experiment  is  essen¬ 
tially  identical  to  the  previous  one  except  that  the 
signal  voltage  is  developed  by  passing  the  diode 


current  through  a  load  resistor,  Rp.  The  TSC  experi¬ 
ment  has  certain  advantages  for  detecting  traps  with 
high  emission  rates  since  the  current  output  is  pro¬ 


portional  to  this  emission  rate.  This  experiment  also 
has  the  most  sensitivity  to  concentration  profiling 
very  near  the  junction,  a  region  of  least  sensitivity 
for  DLTS. 

Details  of  the  theory  of  these  experiments  will 
be  presented  in  a  later  section. 

1 1 .  Previous  Px pe_r_i  men t  a  1  Appa  rat  ns 

We  list  here,  for  completeness,  apparatus  developed 
and  described  previously.*  A  detailed  description  is 
found  in  Ref.  1. 

1.  Bulk  Pool  Variable  'flux  Irradiation  facility 

2.  Neutron  Flux  Detection  and  Integration  System 

3.  Neutron  Activation  Analysis  Systems 

4.  Pour-point  and  Thermal  Probes 

5.  Minority  Carrier  Lifetime  Apparatus 

b .  High  Impedance  Van  dcr  I’auw  Resistivity 
and  I  la  1  1  Sy  s  t  cm 

7.  Sample  Preparation  P.quipmcnt  and  Techniques 

.3  .  Annealing  System 

In  addition  to  the  above,  considerable  microcom¬ 
puter  capability  has  been  developed  at  MURK  using 
Motorola  0800's.  These  systems  are  extremely  inexpen¬ 
sive  (-  $1000  per  computer)  and  have'  contributed 


greatly  to  data  analysis  and  theoretical  calculation 
Our  present  experience  with  these  systems  under  this 
contract  is  mainly  to  use  them  for  computing.  All 
of  the  calculations  for  the  first  section  were  per¬ 
formed  on  one  of  these  instruments.  DI.TS  data  rcduc 
tion  was  also  performed  on  these  microprocessors. 


o  r 


IV.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 


A .  Neutron  Activation  Analysis  --  Determination  of 

In  Concentration 

The  neutron  activation  analysis  technique  and 
facilities  have  been  described  in  a  previous  report.'* 

As  typical,  no  significant  (  =»  10  /cc)  concentrations 
of  trace  impurities  have  been  found  in  float  zone 
material  from  any  source,  however,  significant  con¬ 
centrations  of  contaminants  are  found  in  Czochralski. 

The  first  sample,  #011-167"!)",  was  a  sample 
furnished  by  AEML.  This  sample  was  scanned  for  trace 
impurities,  however,  the  only  impurity  to  be  found  with 
significant  concentration  was  the  dopant  In.  The  est i - 

1  (L 

mated  concentration  is  8,0  x  10  atoms/cc  (6.5  ug/g). 
This  concentration  is  somewhat  larger  than  expected  on 
the  basis  of  electrical  measurements  made  at  AEML.  The 
discrepancy  between  electrical  and  physical  determina¬ 
tions  of  In  is  apparently  a  common  experience  of  several 
laboratories  and  suggests  that  In  may  not  be  100"  elec¬ 
trically  active  in  silicon. 

This  sample  was  irradiated  for  40  hours  in  a  quart; 
vial  along  with  four  liquid  mul t i - c 1 ement  standards. 

The  standards  included  Cr,  Co  and  Sc.  No  attempt  was 
made  to  look  for  short  half-life  isotopes.  After  two 
days,  the  standards  were  opened,  diluted  to  volume  and 
a  known  portion  transferred  to  filter  disks  and  dried. 
The  sample  was  then  opened,  etched  and  transferred  to 

y.o 


a  counting  tube.  Standards  and  the  sample  were  then 
counted.  No  peaks  were  detected  except  those  attributed 
to  ^^mIn.  Using  the  cobalt  standard  for  normalization, 
the  sample  was  compared  to  an  indium  standard  irradiation 
of  6  months  prior  to  this  irradiation  for  an  In  concen¬ 
tration  estimate.  This  is  obviously  not  the  best  exper¬ 
imental  procedure  for  a  determination  of  In,  however, 
no  In  standards  were  run  simultaneously  with  this  sample 
since  an  accurate  determination  of  specific  dopants  was 
not  requested  by  APML. 

A  piece  of  Czochralski  silicon  from  a  rather  large 
ingot  (obtained  from  Monsanto)  which  we  have  used  for 
all  optical  and  Hl’R  work  to  follow  (MIJRR  166)  has 
been  analysed  by  NAA.  This  ingot  is  undoped  and  has 
a  resistivity  of  approximately  280  fl-cm  (n-type)  which 
is  rather  good  for  CZ-silicon.  The  net  carrier  con¬ 
centration  is  only  1.8  x  10*\-m 

The  sample  was  etched  and  crushed  in  an  agate 
mortar  prior  to  irradiation.  Three  portions  ( —  ISO  mg) 
of  the  crushed  sample  were  weighed  into  TO- 8  quartz 
vials  and  irradiated  along  with  flux  monitor  standards 
containing  known  quantities  of  antimony,  cobalt,  and 
chromium  for  40  hrs.  in  position  il-1  (4-8).  After  irra¬ 
diation,  standard  solutions  were  transferred  from  the 
irradiation  containers  and  diluted  to  10  ml;  0.8  ml  was 
deposited  on  Whatman  84  1  filters  on  llandiwrap.  After 
drying,  the  plastic  and  filter  were  folded  and  placed 


Table  7.  NAA  for  MURR  166  (Si:Cz 


In  a  counting  tube.  Samples  were  also  transferred  to 
counting  tubes. 


Calculations  were  performed  determining  counts  per 
sec/ii gm-atoinic  for  each  standard  element  and  normaliza¬ 
tion  constants  were  determined  to  correct  previously 
known  CPS/ygm- atomic  for  some  17  elements  to  the  present 
irradiation  and  counting  conditions.  Copper,  gold, 
antimony,  sodium,  chromium  and  cobalt  were  observed  in 
the  sample  aliquots.  Large  bromine  peaks  were  observed 
and  a  calculation  giving  a  lower  limit  for  bromine  was 
performed.  The  results  are  shown  in  Table  7. 

Two  floatzone  samples,  MURR  134  and  156,  have  been 
obtained  from  Rockwell  for  the  purpose  of  determining 
In  concent  rat  ion  and  in  activity.  The  concentration  of 
indium  determined  by  NAA  was  compared  to  a  determina¬ 
tion  made  by  room  temperature  resistivity  measurements. 
Since  In  is  a  deep  acceptor,  it  is  not  completely  ionized 
at  room  temperature  and  corrections  to  the  electrical 
measurements  must  lie  made  to  account  for  this  partial 
ionization.  The  equation  for  the  hole  concentration,  p, 
as  a  function  of  temperature  for  a  two  acceptor  (In  and 
B)  and  compensating  donor  (N  )  model  was  given 
prev i  mis  1 y  as  * 

(P  +  N'|,1(P  +  KjHp  +  k.)  -  NjKjlp  +  K  ,  1  -  K  ( p  +  k }  )  =  ll 

who  re 

_  ’  V  -  i.  ,/KT  (  i  nd  ium  ) 
n  j  -  —  e  i 


i 

L 


(boron) 


K  =  !!l  e*n7/kT 
_  B° 


Ny  =  2(2um*kT/h:)-l/2  . 

1  a  _  |  y  _ 

For  Nj  ~  JO  “cm  ‘  and  N,,  —  3  0  cm  ,  tire  above  equation 
can  be  approximated  by 

(p  +  NpHp  +  K  j )  ( p  +  FU)  -  N  ^  K  -j  (p  +  KJ  =  0 

and,  neglecting  N,j  compared  to  p  and  canceling  the  fac¬ 
tor  (p  +  K0),  we  find  that 

N,  =  N,  (In)  s  p~  +  Klp 

K'i 

Using  m*  =  0.59  m  ,  B  =  4  and  the  activation  energy  for 
In  we  find  that 


=  (Ny  / p )  e'h/kl  =  5,  45S5xl014  rl/2 

where  T  is  in  degrees  Kelvin.  The  carrier 
p,  can  be  determined  from  the  resistivity 
t i on ship 

p  =  (5x10*^)  |  |  cm 


-  1 S50/T 
e 

cone on t  ra t i on 
using  the  re  la 


Using  the  three  equations  above,  the  In  concentration  is 
given  in  terms  of  the  room  temperature  resistivity. 

NAA  was  used  as  above  to  determine  the  In  concen¬ 
tration.  It  yielded  an  In  concentration  of  2 . 50x  1  0  *  ('cm 
for  both  samples.  (An  indium  standard  was  run  with  each 
sample  in  this  case).  The  results  of  the  electrical 
measurements  are  seen  in  Table  S. 


Table  8.  Comparison  of  In  concent  rat i on  determination  by 

NAA  and  room  temperature  resistivity  measurements 


Samp 1 e 

p(ij-cm)  pfem  ''') 

K^(cw 

| J  n 1  by  NAA 

MURR  134 

1.434  8. 71x10* 5 

2. 5  5x10  16 

2.50  x 1 0  1  b 

MURR  136 

1.4  27  8. 87x10  15 

2 .65x10  16 

2. 50x1 0lb 

It  appears  that  the  two  methods  agree  to  within  about 
5'  which  is  within  the  accuracy  of  the  NAA  determina¬ 
tion.  It  is  therefore  not  clear  why  other  laboratories 
find  a  higher  concentration  of  In  by  physical  methods 
than  by  electrical.  hither  incorrect  theory  is  used 
to  fit  electrical  data  or  the  percent  electrical  activa¬ 
tion  of  In  is  very  concentration  dependent.  It  should 
be  noted  that  the  "  ionization  of  In  at  room  tempera¬ 
ture  can  be  determined  by  100  p/N^  to  be  34.2'  for 
MURR  131  and  33.7'  for  MURR  136.  The  average  is  close 
to  34';..  l\’e  find  the  following,  to  be  a  good  annroxi- 
mat ion,  t  he rc  fo  re , 

V  _  P  _  253(5xl01>l 
A  0.34  0.34  P 

2  53  (  5x  1  0  1  ■'  ) 

0.34  p 

o  r 

\v  =  3.  72x  101(>/p 

where  ..  given  in  ii-cm.  This  simple  equation  will  onl\ 
hold,  however,  when  the  room  temperature  is  such  that 
Kj  —  1.52x10* 'Vm  ''  and  when  the  measuring  current  is 
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sufficiently  low  so  as  not  to  heat  the  sample.  The 

last  equation  should  be  accurate  to  within  in':  for  room 

temperatures  between  15  anti  20°C. 

B .  Radioactive  Decay  of  ^'^"'in 

When  in-doped  silicon  is  irradiated  by  thermal 

neutrons,  the  radioactivity  of  the  material  is  governed 

by  the  metastable  isotope  In  (T,  =  50  days).  This 

effect  can  cause  minor  activity  problems  for  the  case 

of  high  phosphorus  concentration  added  by  NIP  and  high 

In  concentration.  These  arc,  in  fact,  the  circumstances 

for  In-doped  CZ  since  In  concentrations  can  reach 

5x10^'  cm  while  the  residual  boron  to  he  compensated 

might  be  a  factor  of  ten  higher  than  in  1Z. 

Five  In-doped  CZ  wafers  have  been  transmutation 

doped  for  West i nghousc .  The  In  concentration  in  these 

17  -5 

wafers  was  approximately  5x10  cm  ' .  These  wafers  were 

irradiated  to  a  fluence  of  1.967x10  n/ciiT  to  add  about 
1  7  -  5 

5.5x10  cm  ‘  of  phosphorus  for  compensation.  The  samples 
were  irradiated  in  a  flux  ol  about  4.2x10  /cm“/scc  in 
position  RSAF  Z0-Z4.  The  samples  were  removed  from  the 
reactor  on  10-7-78.  the  activity  on  10-51 -"8  was 
measured  as  6.77x10  ">pCi/gm  caused  hv  the  decay  of 
while  the  activity  on  11-0-78  was  6.597xl0  'uCi/gm.  The 
''^mln  was  identified  by  the  gamma  emissions  at  100,560 
and  7Z0  keV.  The  half-life  of  ''  '"'in  is  50  days.  In 
addition,  an  activity  of  5.55x10  'SiCi/gm  of  lVo  was 

i' 1 
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Figure  24.  Estimate  of  radioactivity  from  In  as  a 

function  of  days  after  irradiation  for 
various  In  concentrations  and  F  concentrations 
added  by  transmutation.  The  lower  figur 
experimental  datn_for  [In]  =  5  x  lfP7cm~ 

[PI  =  3  x  10^\:m"3  on  which  estimates  arc 
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detected.  A  plot  of  the  decay  of  ^1^mTn  activity  is 

shown  on  the  lower  half  of  fig.  24  and  confirms  the  half- 

life  of  50  days  to  within  counting  statistics. 

It  should  be  noted  that  approximately  5b  half-lives 

or  275  days  are  required  for  this  sample  to  decay  to 

exempt  limits  of  2x10  ^yCi/gm  for  ^^mIn.  It  must  reach 

this  limit  before  it  can  be  sold  as  exempt  material. 

The  limit  for  research  quantities  (not  to  be  sold)  is 

10uCi  total.  This  limit  allows  approximately  1.5  kgm 

of  similar  material  to  be  shipped  for  research  purposes 

25  days  after  irradiation.  for  typical  f-'Z  material  with 

5x10 Tn  cm  ^  and  the  concentration  of  phosphorus  added 
1  0  -  5 

to  be  5x10  “cm  ’ ,  we  can  estimate  the  decay  period  as 
shown  in  the  top  of  fig.  24  to  be  —  160  days  to  reach 
2x10  pCi/gm.  Research  samples  of  FZ  con  be  shipped 
almost  immediately  under  the  10gCi  total  limit. 

C .  Vertical  flux  Profile  in  RSA1:,  Doping  Accuracy  witli 

Industrial  Laboratories,  and  final  Minority  Carrier 
Li f et ime 

The  University  of  Missouri  Research  Reactor  is  now 
working  with  several  industrial  laboratories  to  produce 
compensated  IR  detector  material  and  material  for  other 
detector  applications.  The  demand  for  these  services 
is  far  greater  than  anticipated  and  has  at  times  created 
problems  for  our  researcii  program  in  that  the  calibrated 
positions  in  RSAT  arc  usually  full  or  being  recalibrated. 
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'•’here  is  also  an  additional  burden  on  research  staff 
since  the  high  precision  irradiation  problems  demand 
considerable  attention  to  calibration,  flux  wire 
counting  and  analysis,  and  special  sample  loading. 

While  the  sample  backlog  has  been  as  high  as  4(1  irra¬ 
diations  with  a  total  reactor  time  of  over  HUM)  hours 
this  detector  work  is  not  steady  enough  to  justify 
full  time  technician  help  at  present.  These  industrial 
experiments  therefore  cut  in  the  research  time  o!  the 
staff  to  some  extent.  We  believe,  however,  that  these 
irradiations  are  in  the  best  interests  of  Al'MI.  and  are 
for  the  most  part  relevant  to  the  goals  of  the  present 
cont  met . 

Irradiations  have  been  performed  on  a  fairly  regu¬ 
lar  basis  for  the  following  organisations: 

1.  Rockwell  International  f  1 R  detectors') 

2.  Hughes  -Car  1  shad  (other  detectors') 

5.  Westinghouse  Research  (IR  detectors) 

4.  Or tec  (other  detectors) 

5.  University  of  Chicago  (CaAs  for  Al'OSR  Contract) 

(> .  Sand  i  a  Laboratories  (NTH  research  samples) 

7.  Monsanto  (Research  samples) 

8.  Komatsu  Metal  fleet ric  (Research  samples) 

R.  Motorola  (Research  samples) 

Under  the  previous  contract',  all  experiments  were 
performed  on  single  wafers  irradiated  at  a  fixed  height 


in  the  variable  flux  facility,  RSAf .  The  extremely  high 
doping  accuracy  achieved  previously  (<il",)  was  the 
result  of  accurate  sample  characterization  prior  to 
irradiation  and  the  fact  that  the  vertical  flux  profile 
of  this  facility  could  be  ignored  as  compared  to  the 
thickness  of  the  sample  wafers. 

Such  a  procedure,  although  demonstrating  the  accu¬ 
racy  of  the  neutron  flux  counting  system,  is  undesirable 
for  large  scale  production  because  of  the  small  sample 
volumes  per  irradiation  (~lgm).  It  is  important,  there¬ 
fore,  to  transfer  this  technology  from  wafers  to  ingots. 
This  will  be  the  subject  of  the  remainder  of  this  sub¬ 
sect  ion. 

The  first  experiment  which  was  performed  was  a 
determination  of  the  vertical  flux  profile  in  RSAF  using 
Co  doped  A1  flux  wire. 

Figure  25  shows  the  flux  profile  obtained  from  flux 
wires  for  the  total  length  of  RSAF.  Thirty  flux  wires 
have  been  irradiated  at  various  heights  and  for  identical 
times.  The  activity  which  is  proportional  to  flux  is 
then  used  to  determine  the  profile.  it  should  be  noted 
from  Fig.  15  that  the  original  position  wc  have  used 
for  all  the  experiments  of  the  previous  contract  (RSAF 
12-16")  is  extremely  linear.  A  second  position,  RSAF 
20-24,  has  also  been  calibrated  for  detector  material 
irradiations.  The  ratio  of  integrator  calibrations  for 


Height  (In.) 


these  two  positions  have  been  compared  w' i  t  h  flux  data 
taken  from  hip.  24  and  is  summarised  in  fable  9 . 


Table  9.  Comparison  of  calibrations  of  two  RSAF  positions 
bv  integrator  counts  and  flux  wire  data  of 
T'ig.  2A\ 


T  NT  1  ill  RATO  R  g.Al.l  ORATION  FLUX  RATIO 

POSITION’  RATIO  (n/cm'/scc)  I  ROM  lie'..  2S 


RSAF  12- lb  1.472  x  I  0  1 

RSAF  20-24  .  , . .  , 

1 . a  a  4  x  10 


-  1.18f. 


Tt  should  lie  noted  that  the  flux  ratio  for  these  two 
positions  agrees  to  within  .06"  for  the  two  different 
determ inat ions . 

The  flux  profile  in  RSAF  12-20  his  also  been 
checked  using  NTD  silicon  wafers,  These  wafers  were 
spaced  in  a  10"  sample  can  using  dummy  silicon  wafers 
and  annealed  after  irradiation  at  8S0°C  in  argon.  The 
results  of  these  experiments  are  shown  in  1 igs.  26  and 
27.  The  solid  line  in  these  figures  is  the  best  lit 
line  to  the  flux  wire  data  in  Fig.  2S.  Although  there 
is  some  scatter  in  data  points,  presumably  due  to 
sample  contamination  problems  during  annealings,  the 
profiles  obtained  in  this  way  agree  with  the  flux  wire 
da  t  a  . 

Lxpcriments  have  also  been  run  for  Rockwell  and 
Hughes  for  stacks  of  wafers  irradiated  i r  RSAF  and 
flipped  end  for  end  at  SO"  of  total  flucnce  so  as  to 
obtain  the  best  profile.  A  set  of  five  12  wafers,  S  1 
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RSAF  12-20 
MURR  168-178 
(Monsanto  Wafers) 


0-  T  I  1  i 
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Phosphorus  Added  (ppb) 

i<;ure  2().  \ertieal  flux  profile  in  RSA!  Jet  e  r:u  i  ne  J  fr 
final  Nil)  wafer  resistivities  (*!URR  Ins -l~s 


RSAF  12-20 
MURR  185-195 
(Monsanto  Wafers) 


Flux  Wire  Position 
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Flux  (Arbitrary  Units) 

Vertical  flux  profile  in  RSAI  determined  from 
final  \TI>  wafer  resistivities  (Ml'KR  1S5-IP5). 


through  S-5,  were  irradiated  for  Hughes -Ha  r  1  shad  .  The 
first  purpose  of  this  experiment  was  to  examine  the  uni¬ 
formity  in  RSAf  20-24  over  a  -1"  range  to  see  if  the  uni¬ 
formity  was  good  enough  to  do  4"  long  ingots.  The  test 
wafers  were  again  spaced  at  equal  intervals  in  the  sam¬ 
ple  can  using  Si  spacers.  The  second  purpose  of  this 
experiment  was  to  check  on  doping  accuracy  relative  to 
Hughes  resistivity  measurements  which  were  made  in  an 
industrial  environment  by  technicians  using  heavily 
utilized  resistivity  probes.  The  anneals  were  also  per¬ 
formed  at  Hughes  since,  ultimately,  an  IR  detector 
program  would  depend  on  in-house  annealing. 

Four  point  probe  measurements  were  made  at  MIJRR 

and  at  Hughes  before  and  after  the  irradiation  and 

annealing.  The  wafers  were  irradiated  to  a  flucnce 

1 7  1 

determined  from  integrator  counts  of  1.235  x  10  '  n/cm“. 
This  corresponds  to  an  added  phosphorus  concentration  of 
0.4  143  ppb  =  2.0  7  x  10  * cm  \  The  concentration  of 
phosphorus  added  was  also  determined  from  initial  resis¬ 
tivity,  p  before  NT  I)  and  final  resistivity,  p  ,  after 
\TD  and  annealing  using  the  equation' 


253  89.54 


(ppb) 


where  the  phosphorus  concentration,  N  ,  is  given  in  ppb 
(1  ppb  =  5  x  10  1 cm  ^ )  . 

Data  obtained  from  MURR  resistivity  measurements  and 
IIUGIII.S  measurements  are  shown  in  Table  Id.  No  effort  was 


10’ 


sistivity  measurements. 


made  to  compare  resistivity  measurements  prior  to  this 
experiment . 


Ke  believe  these  results  are  rather  typical  of 
what  can  be  expected  on  a  production  basis  and  arc 
significant  in  several  respects: 

1.  If  we  assume  that  samples  S- 1  and  S-5  were  con¬ 
taminated  slightly,  the  MURR  data  indicates  that 
the  doping  target  was  achieved  to  better  than  1". 

2.  Annealing  of  NTH  ingots  for  detector  applications 
can  easily  be  achieved  in  an  industrial,  non¬ 
research  laboratory  environment.  If  the  results 

of  single  wafer  anneals  arc  this  good,  whole  ingots 
should  be  even  more  successful.  In  fact,  several 
recent  experiments  working  with  Hughes  has  produced 
compensated  p- type  samples  with  ROPIK  p  JO , OPO  p-cm 
to  an  accuracy  of  better  than  A  (compare  to  our 
earlier  results  in  the  previous  contract  in  making 
high  resistivity  p-type  material')  and  with  a  final 
minority  carrier  lifetime  of  about  550  nsec . 

A.  The  vertical  doping  profile  can  be  held  uniform  ovei 
a  length  of  1"  to  better  than  11.  except  perhaps  for 
some  flux  fringing  at  the  Si -water  interface  at 
each  end  of  the  sample  can. 

4.  To  hit  target  to  within  SI,  relative  to  industrial 
resistivity  measurements  and  annealing  with  no 
special  non- routine  calibration  measurements  is 
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wafers  (B-2  t  li  rough  B  -11). 
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extremely  encouraging  and  of  sufficient  accuracy 
to  insure  success  of  the  1R  detector  program. 

A  similar  set  of  experiments  were  attempted  for  the 
Rockwell  group,  however,  the  results  were  not  as  encoura 
ging  due  to  the  accidental  contamination  mentioned  ear¬ 
lier.  In  this  set  of  experiments,  all  data  was  obtained 
by  MURR.  A  rather  interesting  doping  problem  was  found 
in  this  set  of  wafers.  One  end  of  the  wafer  stack 
(apparently  taken  in  succession  from  a  single  ingot) 
showed  a  rather  steep  doping  gradient.  The  hole  concen¬ 
tration  vs.  wafer  number  is  shown  in  Fig.  2  S .  Beyond 
wafer  B - 1 1 ,  uniformity  from  wafer  to  wafer  was  very  good 
The  radial  profiles  of  these  wafers  was  also  curious. 
These  profiles  are  shown  in  Fig.  29  after  NTD  doping 
and  annealing.  All  wafers  between  B-2  and  B-10  (The 
region  of  axial  non-uniformity)  showed  extremely  low 
resistivity  at  the  edges  of  the  wafer,  and  a  peak  in 
resistivity  near  the  center  which  got  progressively 
broader  down  the  length  of  the  ingot.  The  data  in  Fig. 
29  is  uncorrected  for  image  edge  effects  in  the  probe 
measurements,  however,  these  effects  tend  to  raise  the 
apparent  resistivity  as  the  probes  approach  the  edges. 

We  know  that  N'TD  doping  did  not  contribute  to  this  non¬ 
uniformity  (sec  NTD  doped  profiles  in  Fig.  57  of  Ref.  1) 
however,  a  similar  effect  has  been  seen  previously  in 
a  C’.a  doped  wafer  (see  Fig.  SB  in  Ref.  1). 


Resistivity 


;  i 

100  f 

:  t 


100  O - cm 


1200  Q-cm  — ■/ 


Irradiated  and 
Annealed 


960  O-cm— _ 


^r970  Q  cm 


1  t  t  i  ;  i  i  i  i  i  i  i  |  i  i  i  i  i  1 

0.5  1.0  1.5 

Probe  Position  (in.) 
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gure  31.  Doping  uniformity  ns  n  function  o  I'  wnfet 
hoi ght  position  for  Rockwell  wafers. 
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Figure  30  shows  a  normal  profile  before  N'TD 
for  a  sample  wafer  from  the  uniform  region  of  the 


wafer  stack  (position  B-12).  Samples  B-2,  4,  b,  8, 

10  were  used  to  determine  NTD  vertical  doping  unifor¬ 
mity  as  was  done  for  the  Hughes  group. 

Figure  3L  shows  graphically  the  fluenccs  obtained 
from  the  initial  and  final  resistivities  of  five  addi¬ 
tional  wafers  used  to  determine  the  vertical  flux  pro¬ 
file  in  RSAF.  The  flucnce  was  calculated  from  the 
cquat ion 


N  k 

=  _11  =  _ 
k  k 


W  "J 


-  1  18  2 
where  =  (ep  )  and  K  =  3.355  ppb  per  10  1  n/ctn" 

and  p  and  n  ^  are  the  initial  and  final  resistivities 

before  and  after  NTD.  Wafers  were  all  annealed  at 


850°C  in  flowing  argon  as  usual. 

The  fluencc  profile  obtained  in  this  way  is  shown 
in  fig.  31.  Although  the  linearity  is  excellent,  an 
apparent  doping  gradient  exists  from  top  to  bottom. 

The  unusual  radial  profiles  led  us  to  suspect  that  a 
mobility  gradient  might  exist  along  the  ingot.  We 
therefore  made  Hall  measurements  on  the  two  end  samples, 
B-2  and  B-10,  ami  redetermined  the  irradiation  fluenccs 
on  the  basis  of  these  Hall  data.  The  triangles  shown 
on  fig.  31  are  from  Hall  data  and  agree  with  the  probe 
resistivity  data  to  within  the  accuracy  of  our  Hall 
measurements  (  •  5  1  .  '  furthermore,  the  flucnce  at  mid 
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can  position  (0  mm)  was  about  70"  of  that  expected  from 
the  integrator  counts.  These  facts  suggested  that  the 
can  was  irradiated  at  the  wrong  sample  position  and  was 
not  flipped  end  for  end  as  requested. 

It  was  ultimately  determined  that  the  sample  was 
irradiated  in  position  RSAF  22-26  instead  of  position 
RSAF  12-16  because  of  the  accidental  insertion  of  an 
extra  10"  spacer  can  in  RSAF.  The  gradient  was  caused 
by  not  flipping  the  can  as  requested.  The  profile 
shown  in  Fig.  St,  then  is  characteristic  of  the  flux 
profile  in  RSAF  22-26. 

As  a  final  topic  of  this  subsection,  we  discuss 
briefly  the  effects  of  radiation  damage  and  annealing 
on  minority  carrier  lifetime  in  NTD  silicon.  Two  1" 
long  x  1"  dia.  pieces  from  Hughes  ingot  FO1771101  were 
cut  for  minority  carrier  lifetime  measurements  before 
and  after  annealing.  One  sample,  MURR-98,  was  irradiated 
to  a  thermal  fluencc  of  5.5  x  10  n/cm"  =  0.1174  ppb 
P-added  while  the  other  piece,  MURR-99,  was  not  irra¬ 
diated.  Roth  samples  were  then  annealed  at  850°C. 

Before  irradiation,  both  samples  had  initial  resistivi¬ 
ties  of  about  5450  fi-cm  p-type.  The  lifetime  of  MURR-98 
was  initially  630  usee  while  that  of  MURR-99  was  820  yscc. 
After  irradition  and  annealing,  MURR-98,  the  irradiated 
sample  had  a  final  resistivity  of  1290  :  -cm,  about  87 
lower  than  expected.  Since  the  resistivity  of  these  1" 
long  pieces  varied  from  end-to-end  hv  67,  t  li  i  s  final 


resistivity  was  close  to  that  expected  from  the  fluence. 
The  lifetime  of  this  sample  dropped  from  030  ysec  to 
about  100  nsec,  rather  close  to  values  reported  by  others 
for  NTD-Si. “ 

he  believe,  however,  that  lifetime  is  being  limited 
by  introduction  of  fast  diffusing  impurities  during 
annealing  and  not  by  radiation  damage.  The  unirradiated 
sample,  MURR-99,  which  was  annealed  simultaneously  with 
MURR-98  had  a  change  in  resistivity  to  16123  S2  -  cm  (p-type) 
and  a  final  lifetime  of  50  yscc. 

We  emphasize  again  the  Hughes  result  of  500  pscc 
after  NTD  doping  at  MURR,  and  annealing  at  Hughes. 

TOPSIL  has  also  reported  very  high  lifetimes  for  selected 
samples  and  has  suggested  that  lifetime  is  limited 
almost  entirely  by  annealing  contamination^  In  fact, 
they  suggest  that  the  lifetimes  in  harder  neutron  spectra 
of  light  wafer  reactors  such  as  MURR  are  even  higher 
than  in  soft  spectra  irradiated  material. ^  This  obser- 
vat  ion  has  also  been  made  by  the  Oak  Ridge  group. “ 

Because  of  the  difficulty  of  isolating  large  samples 
in  our  ~1"  I.D.  annealing  furnace  from  the  quartz  fur¬ 
nace  tube,  a  known  source  of  contamination,  a  study  lias 
been  made  to  sec  how  the  observed  lifetime  varies  with 
sample  dimensions.  The  sample  used  for  this  experiment, 
MURR- 102,  was  successively  cut  in  half  and  the  lifetimes 
measured  as  a  function  of  thickness.  The  results  for 


both  pieces  102A  and  102B  arc  shown  in  fig.  52.  (hie 
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Figure  32. 


Lifetime  in  pscc  vs.  sample  thickness  for 
MURR  102. 


would  expect  to  see  problems  with  surface  recombina¬ 
tion  when  the  minority  carrier  diffusion  length  given 
by 


f  =  V/  D  t 
n  v  n  n 

is  of  the  order  of  the  smallest  sample  dimension.  lor 

'J 

D  =  (kT/e)p  =  3b . 4  cm“/sec  and  for  t  =  500  ysec, 
n  v  n  n  ’ 

then  l’  =  1.35  mm.  Tliis  measured  lifetime  should  tlien 
represent  the  bulk  lifetime  for  samples  whose  smallest 
dimension  is  thicker  than  about  three  diffusion  lengths 
or  4  mm . 

The  results  of  Fig.  32  show  that  both  pieces  behave 
in  an  anomalous  way  abruptly  at  this  thickness,  however, 
the  results  of  these  experiments  arc  difficult  to  inter¬ 
pret  beyond  this.  It  is  concluded,  therefore,  that 
lifetimes  must  be  measured  on  ingot  sized  pieces  for 
reliable  results.  Relative  measurements,  t/i  ,  might  be 
possible,  however,  on  smaller  pieces.  Because  of  the 
expense  of  performing  these  experiments  on  whole  ingots 
(  ~ $  1  , 000  -  2 , 000/ i ngot )  and  because  a  new  annealing  fur¬ 
nace  and  tube  of  large  diameter  are  required  (  —  $6,000'), 
we  do  not  intend  to  pursue  these  experiments  further 
without  a  specific  financial  commitment  to  do  so.  he 
believe  that  silicon  producers  are  in  the  best  position 
to  do  these  experiments  because  of  the  cost  of  the 
ina  t  e  r  i  a  1  . 
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While  the  experiments  of  the  previous  sections 
have  been  related  to  the  development  and  extension  of 
NTH  technology,  the  experiments  in  this  subsection  and 
those  to  follow  will  he  devoted  exclusively  to  a  studs' 
of  the  defects  created  by  radiation  damage  and  their 
anneal  i  rig  properties. 

Although  a  considerable  body  of  radiation  damage 
knowledge  exists,  neutron  damage  studies  in  silicon 
have  been  motivated  by  radiation  hardening  of  devices. 
These  experiments  arc  not  particularly  relevant  to  NTH 
material  for  several  reasons.  first,  the  hardening 
studies  are  usually  performed  on  doped  material  while 
the  NTD  staring  material  for  most  devices  is  undoped. 

Since  it  is  well  known  that  defect  production  rates  and 
migration  energies  are  dependent  on  defect  charge  states, 
the  differences  in  l-'ermi  level  position  for  doped  material 
and  NTD  starting  material  are  expected  to  produce  dif¬ 
ferences  in  defect  production  and  annealing.  Further¬ 
more,  in  radiation  hardening  experiments,  thermal  neutron 
effects  arc  always  eliminated  by  Cd  or  B  sample  shielding 
in  order  to  minimize  activity  problems  for  the  experi¬ 
mentalist.  The  most  essential  component  of  our  study, 
the  thermal  neutron,  has  in  the  past  been  systematically 
excluded  from  the  radiation  damage  field. 

We  must  formulate  a  new  radiation  damage  picture 
for  ourselves,  then,  which  is  quite  different  from  the 

1  '  ’ 
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traditional  studies  and  experiments.  In  a  soft  neutron 

spectrum,  nearly  as  many  displacements  are  produced  by 

1  ’ 

gamma  and  beta  recoils  as  by  fast  neutron  knock-ons.  ’  “ 
The  nature  of  the  damage  is  also  quite  different.  The 
fast  neutron  event  is  quite  rapid  creating  cascades  of 
defects  in  extremely  short  time  periods.  The  thermal 
neutrons,  however,  create  radioactivity  which  produces 


recoil  displacements  of  a  few  atoms  per  recoil  from 
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at  a  slow  rate  until  the  radioactivity  becomes  insigni¬ 
ficantly  small  after  many  half-lives.  Tor  In-doped 
material,  this  cooling  off  period  is  of  the  order  of 
weeks  as  we  have  seen  previously. 

Uc  must  also  expect  to  observe  considerably  more 
detail  in  annealing  studies  in  NTI)  material  because  this 
initially  undoped  high  resistivity  material  will  be 
extremely  sensitive  to  rapid  changes  in  the  position  of 
the  lerini  level  which  does  not  occur  in  doped  material. 

A  much  greater  sensitivity  to  small  delect  concent ra 
tions  is  therefore  realised  in  our  experiments.  he 
might  mention  the  additional  complication  that  the 
annealing  will  be  sensitive  to  the  amount  of  phosphoric; 
added  by  transmit  at  i  cm ,  and  hence,  thermal  neutron  I  luenm 
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A  careful  inventory  of  defect  production  rates  and 

defect  annealing  stages  is,  therefore,  required.  We 

have  mentioned  the  TOPS  I L  and  Oak  Ridge  results  which 

suggest  that  a  hard  neutron  spectrum  produces  a  higher 

1 

final  minority  carrier  lifetime.  This  result  is 

totally  unexpected  by  intuition  and  extrapolation  from 
previous  damage  studies  but  emphasizes  the  need  to  re¬ 
assess  our  knowledge  of  neutron  damage  as  related  to 
NT1)  silicon.  In  this  and  the  following  subsections, 
a  beginning  of  this  task  is  presented. 

In  all  the  isochronal  data  to  follow,  anneals 
were  made  in  an  argon  atmosphere  for  IS  min.  at  each 
successive  anneal  temperature.  Resistivities  were 
measured  with  four  point  probes  while  carrier  type  was 
determined  by  a  thermal  probe.  Thermal  neutrons  were 
shielded  from  the  samples  in  some  experiments  to  isolate 
the  fast  neutron  damage  by  using  a  boron  rich  sample 
container.  Beta  recoil  damage  was  isolated  by  irra¬ 
diating  samples  in  the  pneumatic  tube  facility  at  MURR 
which  provided  fast  sample  retrieval.  These  samples 
were  immediately  annealed  for  15  min.  at  S5H”f  to  remove 
last  neutron  damage  and  then  allowed  to  sc  1 f - redamage 
from  the  remaining  beta  activity. 

An  examination  of  the  dependence  el  isochronal 
annealing  on  thermal  neutron  dost'  is  shown  in  I  ig.  55. 
All  three  samples  were'  cut  from  the  same  liign  purity 
float  zone  ingot  (5J5u  cm,  p-type),  and  ro».<.-i\ed  thcr 
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mal  doses  of  1.07  x  10  n/cm“  for  MIJRR  75  (0); 


3 .  5  ('  x 


lt)^n/cm“  for  MIJRR  ”6  (A);  and  1.5  7  x  lo'\i/cnr  for 
MlIRR  92  (LI).  In  order  to  perform  these  irradiations 
over  three  orders  of  magnitude,  it  was  necessary  to  use 
positions  of  different  flux  and  different  thermal  - to- 
fast  ratios.  The  cadmium  ratios  of  those  positions 
used  for  sample  irradiation  were  50:1,  12:1,  and  10:1 
respec  t i ve ly . 

All  three  samples  of  Tig.  53  were  p-type  at  the 
start  of  the  anneal  and  converted  to  n-type  at  some 
higher  temperature.  In  fact,  a  thermal  probe  check  of 
50  ingots  irradiated  to  a  variety  of  doses  always  indi¬ 
cated  p-type  after  irradiation  in  contrast  to  Ref.  12. 
The  type  conversion  annealing  temperature  depended  upon 
the  neutron  dose  and,  therefore,  was  dependent  upon  the 
amount  of  transmuted  phosphorus  added.  Sample  (0)  which 
was  irradiated  to  balance  the  boron  concentration  with 
phosphorus,  never  fully  type  converted,  but  attained 
mixed  conduction  type  above  800°C.  Samples  (A)  and  ( f J ) 
type  converted  at  725°C  and  at  l2S°f  respectively.  This 
suggests  that  the  transmuted  phosphorus  does  not  become 
electrically  active  until  higher  annealing  temperatures 
and  that  this  temperature  is  lower  for  larger  phosphorus 
c  oncen t  ra t  i  on  s . 

At  about  POUT  an  extremum  in  the  resistivity  can 
he  observed  in  fig.  35.  IV  e  ideal  i  I  y  this  extremum  with 
the  production  of  acceptors 
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resulting  from  damage  plus 


annealing  and  will  call  these  defect  structures  "O00°l!- 

acceptors."  Since  the  initial  resistivity  was  5250  -cm 
1  ?  _ 

(p  =  2.4  x  10  “cm  '),  the  resistivity  minimum  in  samples 
(0)  and  (.a)  corresponds  to  an  acceptor  concentration  of 
■v  0.3  x  10*'\m  or  an  increase  of  about  (>  x  lO^cm 
This  concentration  is  at  least  an  order  of  magnitude 
higher  than  any  electrically  active  impurity  concentra¬ 
tions  in  the  starting  material.  Sample  (□'  is  not  as 
easily  analyzed  since  it  is  n-typc  at  600“C,  but  if  we 
assume  that  the  acceptor  has  balanced  the  donor  concen¬ 
tration  at  600°C  to  raise  the  resistivity  to  the  intrin¬ 
sic  value,  we  need  only  approximate  the  resistivity  wo 
would  expect  if  the  acceptor  were  absent.  A  straight 
line  drawn  from  the  resistivity  value  at  475°C  to  the 
value  at  675°C  yields  a  resistivity  value  at  000°C  of 
about  90  Q-  cin  or  a  donor  concentration  of  about  5  x 
10  ‘’em  Since  the  "000°C- acceptor"  has  balanced  these 
donors,  the  concentration  of  the  "(>00c  C  -  accept  or"  is 
also  about  5  x  lO^cin  \  This  rather  crude  analysis  has 
shown  that  the  ”000 °C - accept o r"  is  not  strongly  depen¬ 
dent  on  a  neutron  dose.  In  samples  from  other  sources, 
we  have  observed  that  the  concentration  of  ”o00°C- accep¬ 
tors”  is  an  order  of  magnitude  higher  suggesting  that 
this  defect  complex  is  sample  dependent  rather  than  dose 

dependent . ^  Shallow  acceptors  resulting  from  annealing 

[  1 

to  these  temperatures,  in  concentrations  of  between  10 
and  lo'Vm  have  been  observed  by  other  authors.'^  ' 
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It  should  also  be  noted  that  the  isochronal  annealing 

for  these  float  zone  samples  shows  considerably  more 

annealing  structure  than  has  been  previously  reported 
1 2 

for  float  zone. 

The  question  of  cadmium  ratio  was  examined  by 
irradiating  three  samples  in  the  same  position,  with 
and  without  thermal  neutron  shielding.  Figure  34  shows 
the  results  for  samples  cut  from  the  same  ingot  of  float 
zone  silicon  (5500  p-cm,  p-type).  This  ingot  differs 
from  that  used  for  the  experiments  shown  on  Figure  33. 
Samples  (0)  and  (□)  ,  which  were  shielded  from  thermal 
neutrons,  remained  p-type  during  the  entire  anneal,  while 
the  sample  which  was  not  boron  shielded  type  converted 
from  p-type  to  n-type  between  700°C  and  725°C  apparently 
due  to  transmuted  phosphorus.  Two  qualitative  observa¬ 
tions  can  be  made  concerning  the  anneals  of  Fig.  34. 
First,  the  major  annealing  features,  other  than  a  slight 
shift  in  the  curves,  do  not  seem  to  depend  upon  the 
presence  of  thermal  neutrons  during  irradiation.  Second¬ 
ly,  the  temperature  at  which  the  resistivity  minimum 

occurs  is  lower  for  smaller  fast  neutron  doses.  Again 
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the  "600°C-acceptor"  concentration  is  ^  6  x  1 0  'em  '  and 
is  apparently  due  to  fast  neutron  damage.  It  should  be 
noted  that  a  study  of  damage  due  to  recoil  effects  would 
be  difficult  in  the  presence  of  the  fast  neutron  damage 
since  this  damage  tends  to  dominate  the  annealing  charac¬ 


ter  i  s  t  i  c  s  . 


Figure  35  shows  the  comparison  of  isochronal 

annealing  of  a  Czochralski  silicon  sample  (50  :-cm,  n- 

type),  thermal  neutron  dose  =  3.80  x  10  ln/cin“,  to  the 

heavily  irradiated  float  zone  sample  shown  previously 

in  Fig.  32.  1'he  two  anneals  arc  surprisingly  similar 

both  in  the  type  conversion  point  from  p-type  to  n-type 

at  400°C  and  in  the  rise  in  acceptor  concentration  to 

about  5  x  10 1  \'in  at  600°C.  IVe  can  directly  compare 

1 

t  he  sc  results  to  those  of  Kharchenko  ct.  a  I  .  “  IVltile 
the  results  for  Czochralski  silicon  agree  quite  well, 
those  for  float  zone  are  very  different.  Furthermore, 

their  float  zone  sample  was  initially  n-type  after  irra- 
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diation  while  ours  are  p-type. 

In  measuring  the  resistivity  after  each  annealing 
step,  a  rise  in  the  resistivity  as  a  function  of  time 
was  noted.  A  measurement  of  this  rise  was  taken  for  the 
heavily  irradiated  float  zone  sample  of  Fig.  33.  Data 
for  the  first  60  minutes  alter  cooling  the  sample  to 
room  temperature  is  shown  in  Fig.  36  and  is  arranged  in 
order  of  annealing  temperature  (this  figure  does  not 
indicate  the  absolute  change  in  resistivity  between 
annealing  steps,  however,  all  relative  changes  arc  to 
the  same  scale).  A  large  effect  can  be  noted  between 
375°C  and  500°C  which  corresponds  to  the  type  conversion 
peak . 

IVhile  this  effect  may  have  something  to  do  with  the 
initiation  of  phosphorus  electrical  activity,  such  a 
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conclusion  cannot  be  drawn  without  direct  examination 
of  energy  levels  in  the  band  gap.  The  mechanism  for 
this  resistivity  instability  is  not  understood  at  this 
time.  We  have  observed  this  effect  in  samples  using 
soldered  contacts  as  well  as  probes  and  in  etched  sam¬ 
ples  as  well  as  lapped.  it  should  also  be  noted  that 
this  effect  is  very  large  from  375  to  500°C  but  is  neg¬ 
ligible  in  the  annealing  range  of  550-650°C  even  though 
the  resistivity  is  similar  in  magnitude  in  these  two 
annealing  peaks.  This  effect  might  be  due  to  carrier 
trapping  at  defects  with  extremely  small  capture  cross 
section,  to  very  slow  room  temperature  migration  of 
defects  freed  from  complexes  during  the  anneal,  or  even 
due  to  electrical  activation  of  donor  oxygen  complexes 
which  are  normally  only  observed  in  Czochralski.  further 
experimentation  will  be  required  to  identify  the  proper 
mcchan ism. 

figures  37  and  38  arc  isochronal  anneals  of  two 

float  zone  silicon  samples  which  had  been  beta  recoil 

damaged  by  the  technique  described  previously.  The 

sample  of  fig.  37  (1200  n-cm,  p-type)  was  irradiated  to 
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a  thermal  neutron  dose  of  3.84  x  10  n/cm“  making  the 
boron  concent  rat i on  approx  imately  equal  to  the  concen¬ 
tration  of  transmuted  phosphorus.  In  Fig.  38,  the  sam¬ 
ple  (5250  ;;-cm,  n-typej  had  a  dose  of  1.02  x  101  n/cuf  , 
making  the  phosphorus  to  boron  ratio  about  10:1.  In 
contrast  to  previous  samples  which  were  p-tvpc  as  a 
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result  of  fast  neutron  damage,  both  of  the  samples  were 
n-tvpe  after  the  recoil  damage  was  completed.  Since 
both  samples  are  n-type  at  the  start  of  the  anneal,  we 
can  make  an  estimate  of  the  donor  defect  concentration 
due  to  beta  recoil  events  by  assuming  that  the  mobility 
has  not  drastically  changed.  The  starting  material  in 
fig*  37  has  NA  -  1  x  10^‘cm  and  after  beta  recoil 
damage  n  1.8  x  1012cm'3,  so  N’u  -  N'A  +  n  =  1.18  x  1013cm 
In  Fig.  38,  Na  =  2.4  x  1012cm'J>,  n  *  8. 93  x  1012cm'\ 
so  Ny  -  1.14  x  10  ‘cm  .  Therefore,  a  donor  concentra¬ 
tion  of  about  lO^cm  *  seems  to  be  produced  by  the  beta 
recoil  damage  in  both  samples.  Other  features  of  these 
anneals  are  not  as  directly  comparable  and  apparently 
depend  upon  either  the  boron  to  phosphorus  ratio  or  dose. 
We  can  sec,  however,  the  appearance  of  an  acceptor  or 
the  disappearance  of  the  initial  donor  in  annealing 
between  300°C  and  400°C  in  both  samples. 

IVc  believe  that  these  p  recoil  experiments  arc 
the  first  concrete  experimental  attempts  to  isolate  the 
effects  of  rccoi  1  damage  from  those  of  fast  neutron 
damage.  Although  the  concentrations  of  recoil  defects 
which  can  be  produced  are  too  low  for  FI'R  or  opt  ical 
techniques  to  be  experimentally  successful,  Pl.TS  should 
be  useful  once  correlations  of  PITS  energy  levels  have 
been  made  with  known  defect  structures. 


In  summary,  we  have  surveyed  the  isochronal  annealing 

of  both  float  zone  and  Czochralski  N'TD  silicon  under  a 

variety  of  irradiation  conditions.  A  discrepancy  with  some 
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of  the  previous  literature  has  been  found  for  NTD  float 
zone  both  as  to  the  carrier  type  after  irradiation  but  before 
annealing  and  in  the  amount  of  annealing  structure  and  final 
annealing  temperatures. 

We  have  found  that  the  dominate  fast-neutron  damage- 

pi  us- anneal ing  effect,  in  terms  of  the  concentration 

of  electrically  active  defects,  is  the  production  of  the 

''600°C  acceptor."  Although  these  acceptors  have  been 

reported  by  others  in  float  zone  material,  we  have  been  able 

to  show  that  this  defect  also  occurs  in  Czochralski,  that  its 

concentration  is  reasonably  independent  of  dose,  and  hence 

phosphorus  concentration,  and  that  the  concentration  of  those 

defects  appears  to  be  sample  dependent.  Since  this  defect 

occurred  in  high  purity  silicon  in  concentrations  of 
IS  -  3 

5  x  10  "cm  or  higher,  it  is  clear  that  these  defects  do  not 
involve  impurities  other  than  perhaps  oxygen  or  carbon. 

We  have  also  observed  an  instability  in  the  resistivity 
as  a  function  of  time  after  anneal  which  is  most  pronounced  in 
the  temperature  range  of  375-500°C.  We  have  suggested  several 
possible  mechanisms  for  this  resistivity  instability,  however, 
further  experimentation  will  be  required  to  determine  which, 
if  any,  of  these  models  is  true. 

Finally,  we  have  been  able  to,  for  the  first  time, 
isolate  the  damage  and  annealing  effects  of  beta  recoil 


damage  from  the  other  damage  mechanisms  in  the  NTD  process. 

This  damage  mechanism  clearly  produces  donors  before  annealing. 

By  comparing  two  samples  with  different  initial  doping 

conditions,  we  have  been  able  to  show  that  the  "beta  recoil 

donor  defect"  concentration  is  of  the  order  of  10  cm  , 

again  about  an  order  of  magnitude  larger  than  the  largest 

of  the  electrically  active  impurities  in  detector  grade  material. 

The  discrepancy  between  our  annealing  data  recovery 
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temperatures  850°C)  and  those  of  Kharchenko  et  al, 

(<  500®)  led  us  to  investigate  different  possibilities.  It 

would,  of  course,  be  very  desirable  to  obtain  annealing 

temperatures  below  500°C  to  keep  contamination  to  a  minimum. 

In  a  previous  report1  we  had  offered  the  hypothesis  that  the 

high  doping  levels  of  Ref.  12  could  have  perhaps  driven  the 

material  amorphous.  Although  still  a  possibility  since  we 

do  not  know  the  reactor  spectrum  of  the  experiment  in  Ref.  12, 

subsequent  calculations  by  H.  Stein  suggest  that  this  is  not 
2 

likely.  Also,  our  data  in  Fig.  33  and  Fig.  35  for  the 
heavily  irradiated  samples  fail  to  show  this  lowering  of  anneal 
temperature  to  any  significant  degree. 

We  therefore  felt  that  irradiation  temperature  might  be 
a  significant  factor.  A  piece  of  FZ  fMURR  159)  was  irradiated 
at  Argonnc  National  Laboratory  in  the  CP-5  reactor  in  their 
newly  constructed  hot  irradiation  facility.  The  sample 
temperature  during  irradiation  was  496°C,  greater  than  any 
temperatures  expected  from  normal  gamma  heating.  This 
sample  was  subsequently  isochronally  annealed.  The  results 
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are  shown  in  Fig.  3!'.  The  annealing  is  similar  to  our 

previous  results  above  the  irradiation  temperature .  The 

"600°C  defect"  minimum  is  still  seen  followed  by  a  type 

conversion  peak.  Below'  500°C,  slight  annealing  structure 

associated  with  radioactivity  induced  displacements  is  just 

barely  detectable.  We  must  conclude,  then,  that  elevated 

irradiation  temperatures  do  not  explain  the  Russian  data. 

We  will  now  briefly  discuss  recent  experiments  on  annealing 

unirradiated  silicon,  primarily  for  the  purpose  of  oxygen 

detection.  It  has  been  known  for  many  years  that  CZ-silicon 
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contains  oxygen  in  the  5  x  10  to  1  x  10  cm  range  and 
that  this  oxygen  can  be  activated  electrically  by  heating  at 
450°C  for  periods  of  hours. ^  By  annealing  to  still  higher 
temperatures ,  these  oxygen  containing  donors  disappear.  Kaiser, 
Frisch  and  Reiss  have  shown  that  the  activation  process  is 
fourth  order,  the  initial  rate  of  donor  formation  being  given 
by  the  equation 

(dA4/dt)t(Uo  =  kfAT 

where  is  the  concentration  of  oxygen  donors  and  A  j  the 

concentration  of  "atomically  dispersed  oxygen."^  They 

suggest  the  value  of  k^.  at  T  =  450°C  is  about  1.2  x  10  ^ 
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cm  sec  .  If  one  anneals  at  450°C  for  longer  and  longer  times, 
the  concentration  of  donors  produced  reaches  a  maximum,  A^*, 
which  is  given  by  the  approximate  equation^ 

•V  *  fW • 


They  give  an  approximate  value  at  450°C  for  k^/k 
“38  6 

5  x  10  cm  .  Finally,  on  heating  to  higher  temperatures, 
they  found  that  the  donors  disappear  again  with  an  approx i mat e 1 \ 
first  order  activation  energy  of  2.8  eV,  i.e. 

A4  Ct  D  ^  A4VKt 

where 

K  =  K  e-2'8  eV/kT  . 
o 

The  constant,  K  ,  is  determined  approximately'*'^  from  the 
fact  that  K  =  10  **sec  when  1000/T  =  1.44,  or 

K  =  (lO^sec'^e2-8^8*6  x  1  O'  5)  C694 . 4) 

=  4  x  lO^sec  k . 

It  should  then  be  possible  to  determine  electrically  the 
concentration  of  oxygen  by  studying  the  anneal  properties  of 
unirradiated  material.  It  must  be  mentioned  that  the  above 
relationships  are  only  approximate  kinetics  equations  which 
have  been  linearized  over  a  restricted  temperature  range. 

Also,  these  authors  warn  us  that  this  technique  is  not  likely 
to  work  with  FZ  silicon  since  the  concentration  of  other 
oxygen  traps  in  FZ  might  be  of  the  same  order  of  magnitude 
as  the  concentration  of  oxygen.  ^  This  last  hypothesis  has 
apparently  never  been  tested  experimentally. 

We  therefore  decided  to  run  preliminary  isochronal  anneals 
in  FZ  to  see  if  an  increase  in  donor  concentration  could  be 
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found  near  450°.  Since  we  expect  very  little  change  in 
donor  concentration  in  FZ  ,  high  resistivity  material  was  a 
necessary  reo^i rement .  Kaiser,  Frisch  and  Reiss  also 


report  that  several  trapping  levels  have  been  found  associated 
with  oxygen  donor  activation.16  We,  therefore,  investigated 

the  minority  carrier  lifetime  as  well,  both  with  silicon 
filtered  light  <  Ey. )  and  unfiltered  xenon  light.  The 

results  of  this  experiment  are  shown  in  Fig.  40  and  in  more 
detail  (below  600°C)  in  Fig.  41.  The  material  used  in  this 
experiment  was  undoped  p-tvpe  5000  0-cm  FZ  (MIJRR  179). 

It  can  be  seen  in  these  figures  that  an  increase  in 
donor  concentration  which  reaches  a  maximum  during  15  min. 
isochronal  anneals  at  550&C  is  observed.  Although  this 
maximum  occurs  at  a  higher  temperature  than  the  450°C 
maximum  discussed  in  Ref.  16,  this  is  to  be  expected  since 
the  time  to  reach  a  maximum  at  450 °C  is  of  the  order  of 
100  hrs.  This  maximum  in  resistivity  is  accompanied  by  a 
decrease  in  minority  carrier  lifetime  of  almost  809o  as  expected. 
The  higher  temperature  peak  at  90n°C  in  Fig.  !1  is  believed 
to  be  due  to  furnace  contaminat  im  and  has  been  seen  previously 
in  un irradiated  p-type  material.1 

A  working  hypothesis  is  that  the  maximum  donor  concentration 
should  be  reached  at  a  somewhat  higher  temperature  than 
450°C  during  15  min  isochronal  anneals.  This  means  the  rate 
constants  given  in  Ref.  16  arc  not  unite  correct  at  our  maxi 
mum  donor  t  ciimc  rat  u  re .  for  lack  ol  better  constants,  however, 
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The  peak  near  5  00°C  is  believed  to  be  due 
to  oxygen  activation  while  ttie  D00°C  peak 
is  nrohablv  due  to  furnace  contamination. 


we  will  use  them  as  a  rough  estimate  of  oxygen  concentration 
in  FZ.  Therefore,  from  Fig.  41, 

p  =2.5x10  cm 

r  o 

12  -  3 

p  =2.2xl0icm 

^max 

*  11-3 

#  oxygen  donors  =  =  3  x  10  cm 

Using  the  previous  relations,  then 

(Aj-)5  ■  ...  -  3  x  X0U  -  6  x  1048 

(k,/k  )  5  x  10 

k  f  p 

o  16-3 

A-j  =  initial  oxygen  concentration  -  1.8  x  10  cm 

This  is  not  an  unreasonable  concentration  of  oxygen  in  FZ  and 
is  just  at  the  limit  of  detection  by  optical  absorption  at 
9  pm  because  of  interference  of  a  lattice  band  with  an 
absorption  coefficient  of  the  same  size  as  the  absorption 
coefficient  due  to  this  oxygen  concentration.  To  complete 
our  analysis,  the  decay  of  donors  should  be  complete  by 
600°C  =  893 °K  as  we  observe.  That  this  is  true  is  easily 
seen  from 

K  ■  K  e‘Z-8  eV/kT 
0 

,,  in14  -K  -2. 8/(8. 6  x  10'5)  (893) 

=  (4x10  sec  ) e 

=  5.86  x  10  “sec 


There  fore  , 


\(t)  '  A4*e'Kt 

•  (1.8  x  10n)e'(S-86  x  10‘2H15  X  60) 

=  a . 8  x  1 o1 1 )  (1.25  x  10'23) 

^  0 

as  expected.  If  we  investigate  this  decay  of  donors  at  500°C, 
we  can  see  why  the  isochronal  annealing  is  still  increasing. 

At  500°C  =  773°K,  e  ^  =  0.832  and  therefore  very  little  decay 
is  occurring,  however,  the  donor  concentration  is  growing 
as  the  fourth  order.  It  is  impossible  to  improve  on  these 
calculations  without  a  knowledge  of  the  temperature  dependance 
of  the  production  factor,  k^-.  Our  only  recourse  is  to  do 
isothermal  anneals  at  450°  for  hundreds  of  hours.  It  is 
clear,  however,  from  the  above  analysis  that  the  increase  in 
resistivity  between  400  and  550°C  in  Fig.  -)  ]  is  not  inconsistent 
with  oxygen  donor  production.  We  believe  that  this  is  the 
first  time  that  oxygen  concentrations  have  been  found  in  F2 
by  this  method.  By  using  NTD  techniques  to  compensate  to 
higher  resistivities,  greater  sensitivity  could  be  obtained 
and  the  possibility  of  detecting  oxygen  in  FZ  to  very  low 
concentrations  could  result.  The  annealing  of  the  NTD  would 
have  to  exclude  the  400-600°C  temperature  range,  hbwevcr,  so  as 
not  to  activate  the  oxygen  prior  to  analysis. 


F. .  Defect  Detection  by  Optical  Absorption  and  F.PR 

We  will  discuss,  in  this  section,  the  eight  defect 
structures  which  we  have  been  able  to  observe  in  XTD 
silicon  to  date.  Although  this  list  of  point  defects 
does  not  exhaust  all  known  defects  which  are  observed 
after  room  temperature  irradiation,  it  covers  most  of 
them.  With  the  exception  of  interstitial  oxygen,  all 
the  defects  are  thought  to  contain  vacancies.  The 
silicon  interstitial  centers  have  not  as  yet  been  observed. 
We  will  count  the  number  of  vacancies  contained  in  these 
centers  at  the  end  of  this  secion  and  compare  this  to 
the  number  of  displacements  we  calculated  in  the  theory 
section. 

Figure  41  shows  the  atomic  configurations  of  three 
of  the  four  defects  we  have  observed  with  optical 
absorption,  namely,  the  oxygen  interstitial  at  9  pm, 
the  A  center  vacancy  oxygen  complex  at  12  pm  and  the 
divacancy  center  at  1.8  pm  and  3  urn.  The  other  observed 
defect  is  an  absorption  band  at  710  cm  ^  which  results 
upon  annealing  and  is  one  of  the  higher  order  bands 
fHOR's')  observed  by  others. 

Our  first  attempt  to  observe  radiation  damage 
with  room  temperature  optical  absorption  was  a  quick 
experiment  on  unctched  02  10  mil  thick  wafers.  The 
results  are  shown  in  Fig.  1  .  The  divacancy  band  at  1.8  am 
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is  clearly  seen  after  a  thermal  fluence  in  row-2  of 
19  2  1 7 

8.89  x  10  'n/cm  .  An  unirradiated  FZ  transmission 

spectrum  is  shown  for  comparison.  An  absorption  near 

3.3  um  was  also  observed  and  originally  thought  to  he 

one  of  the  3  pm  divacancy  bands,  however  subsequent  etching 

removed  this  feature.  It  is,  therefore,  apparently 

related  to  surface  states,  a  rather  surprising  effect 

to  be  found  with  optical  absorption.  This  first 

experiment  allowed  us  to  determine  sample  thickness 

(2-3  mm)  and  fluences  for  subsequent  experiment^. 

Figures  11  and  13  show  the  before  and  after 

irradiation  near  band  edge  absorption  of  FZ  and  CZ 
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irradiated  to  a  thermal  flucnce  of  1.56  x  10  n/cm". 

These  samples  were  mechanically  polished  and  etched. 

Both  the  shifting  to  longer  wavelengths  of  the  band 
edge  and  the  1.8  pm  divacancv  band  can  be  observed.  The 
3.3  pm  surface  absorption  is  not  present  in  these 
properly  prepared  samples.  \To  other  defect  bands  were 
seen  in  these  lightly  irradiated  samples  except  the  oxygen 
interstitial  band  at  9  um  in  the  CZ  samples  as  shown 
in  Fig.  1  (>. 

The  solid  line  in  Fig.  In  shows  the  oxvgen 

18 

interstitial  band  after  the  1.56  x  10  n/cm"  irradiation 
while  the  dotted  band  shows  a  similar  s nmole  after 
annealing  to  4",5°C  for  15  min.  A  shift  to  dtnr tor- 
wave  lengths  of  the  absorption  peak  is  seer..  The 
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absorption,  in  general,  occurs  at  a  wavelength  shorter 

1  ft 

than  9  pm  because  the  data  were  obtained  at  77°K. 

The  double  peak  structure  of  Ref.  18  could  not  be 
observed  in  these  irradiated  samples. 

Figure  -17  show's  the  behavior  of  the  9  pm 
interstitial  oxygen  band  at  77°K  as  a  function  of 
annealing  to  475°C  and  then  500°C.  After  the  first 
anneal,  the  formation  of  the  71C  cm  ^  HOE  is  observed, 
however,  no  other  of  these  bands  could  be  found.  The 
oxygen  band  shifts  to  its  shortest  wavelength  after 
annealing  to  475°C  and  also  recovers  somewhat  in 
intensity.  A  subsequent  anneal  to  500°C  causes  this 
band  to  loss  intensity  and  shift  to  longer  wavelengths 
again.  This  is  shown  in  detail  in  Fig.  1,8. 

We  believe  the  sequence  of  events,  here,  is  as 
follows.  After  irradiation  but  before  annealing, 
oxygen  interstitials  are  lost  by  attracting  mobile  defect 
and  the  intensity  of  the  9  pm  band  decreases.  After 
annealing  to  475°C  for  20  min.  these  oxygen  defect 
complexes  dissociate  to  form  other  defect  complexes, 
the  HOB's,  and  to  return  oxygen  interstitials  to  the 
9  pm  absorption.  On  heating  to  500°C,  however,  the 
oxygen  interstitials,  as  we  have  previously  seen 
coagulate  into  oxygen  donor  complexes,  hence,  the  loss 
of  9  pm  absorption  again  as  expected . ^ ^  ^  The  slight 
shifts  in  frequency  arc  unreported  to  our  knowledge  and 
of  unknown  origin. 
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This  hypothesis  is  justified  from  the  isochronal 

annealing  in  Fig.  11  which  shows  that  maximum  oxygen 

donor  formation  occurs  during  isochronal  annealing  at 

500°C  or  above.  We  must  now  ask,  what  are  the  defects  which 

remove  oxygen  interstitials  from  the  9  ym  band  during 

irradiation  but  before  annealing?  The  obvious  answer 

is  that  oxygen  interstitials  are  lost  during  irradiation 

19  20 

to  A-center  vacancy-oxygen  complexes.  ’  We  therefore 
investigated  the  room  temperature  optical  absorption 
near  12  ym  as  is  shown  in  Fig.  49.  Although  no  clear 
absorption  band  is  seen  at  this  neutron  dose,  an  increase 
in  transmission  is  seen  at  shorter  wavelengths  and  a 
slight  decrease  is  found  at  longer  wavelengths  near 
12  ym.  On  annealing  to  475°C,  these  irradiation 
induced  changes  reverse  indicating  the  loss  of  these 
vacancy  oxygen  centers  and  the  return  of  additional  oxygen 
interstitials  to  the  9  ym  absorption. 

Because  of  the  above  analysis,  we  felt  that  the 
features  shown  in  Fig.  49  near  12  ym  were  indeed  related 
to  the  formation  and  annealing  of  A-centers.  To  prove 
this,  a  second  CZ  sample  from  the  same  ingot  was 
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irradiated  to  a  higher  thermal  dose  of  1.34  x  10  n/cm“ 
in  row-2.  The  raw  data  for  the  optical  absorption  near 
12  ym  is  shown  in  Fig.  SO  for  this  sample.  The  room 
temperature  data  is  shown  in  the  upper  spectral 
trace  while  this  same  absorption  is  shown  at  77°K  in 


I  I  1 

10.0  9.5  9.0 

Drum  Number 

I  iyure  49.  lirst  hint  of  A-contcr  ( vacancy - oxyren 

complex)  at  12  pm  as  irradiated.  Anneal  in 
to  4~5°(!  return1-  the  spectinn  near  12  .,r 
to  it>  nre  -  i  r  rad  i  a  t  i  on  heha  v  i  o  r . 


Relative  Intensity 


the  lower  trace.  The  A-center  band  is  now  clearly 
seen  for  the  low  temperature  data.  Furthermore,  a 
plot  of  the  absolute  transmission,  I/I  ,  also  shows 
the  12  pm  absorption  for  the  room  temperature  trace 
as  shown  in  the  inset.  The  strange  shape  of  the 
room  temperature  absorptions  is  caused  by  an  underlying 
lattice  absorption  band  which  decreases  in  intensity 
as  the  sample  temperature  is  lowered. 

Do  these  A-centers  steal  oxygen  interstitials 

from  the  9  pm  band?  They  do  as  shown  in  Fig.  SI 

where  the  decrease  in  this  band  is  clearly  seen.  The 

lower  curve  is  before  irradiation  while  the  upper  curve 

is  after  irradiation.  Using  this  figure,  we  can  estimate 

the  concentration  of  A-centers  formed.  The  change  in 

absorption  coefficient  for  the  9  pm  band  is  0.22  cm 

This  corresponds  to  a  decrease  in  oxyger  interstitials 

of  about  1  x  101  cm  using  the  data  of  Kaiser  and  Keck 

to  determine  the  concentration  of  oxygen.  The  oxygen 

concentration  before  irradiation  is  also  determined  to 
17  -  X 

be  6  x  10  cm  .  We  can  also  estimate  the  concentration 
of  A-centers  by  assuming  that  the  oscillator  strength 
for  A-centers  and  9  pm  oxygen  interstitials  is  the 
same.  Doing  this  then,  and  using  the  data  in  the 

insert  of  Fig.  49,  the  difference  in  absorption  coefficient 
between  maximum  and  minimum  is  about  1  cm  ^  which  also 
yields  about  1  x  10  cm  A-centers.  We  therefore  sec 


that  1  out  of  6,  or  about  171  of  the  oxygen  interstitial:-- 
available  are  saturated  with  vacancies.  Do  these 
oxygen  interstitials  getter  vacancies  and  deplete  the 
concentration  of  other  vacancy  centers?  To  answer  this 
question  we  must  investigate  the  defects  seen  by  F.PR. 

All  of  the  optical  absorption  data  were  obtained  on 
samples  cut  from  two  similar  resistivity  ingots.  The 
Monsanto  FZ  ingot  (MURR  167)  ,  and  the  Monsanto  CZ  ingot 
(MURR  166).  Both  ingots  are  approximately  250  P-cm 
n-type  before  irradiation.  We  felt  it  important  to 
match  Fermi  level  positions  in  the  RPR  samples  since  the 
occupancy  of  the  defect  levels  will  affect  to  some  extent 
that  which  is  observed. 

Preliminary  irradiations  were  performed,  however, 
on  FZ  samples  from  the  Wacker  ingot  which  we  have 
described  previously.*  This  is  a  very  high  resistivity 
p-type  sample  (p  >  10,000  P-cm)  and  is  substantially 
different  from  the  optical  absorption  samples.  After 
irradiation  to  a  thermal  neutron  fluence  sufficient  to 
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produce  25  P-cm  NTD  material  (1  -  1  x  10  n/cm"),  the 

vacancy- phosphorus  P-center"  was  found  on  cooling  the 
Wacker  sample  below  room  temperature.  The  I'.PR 
spectrum  for  the  I:  center  is  shown  on  Fig.  Also 

shown  on  this  figure  is  an  atomic  model  of  this  defect. 
The  l;PR  resonance  is  cha  rac  t  er  i  2  ed  by  the  two  strong 
hvperfine  doublets  indicated  by  dashes  or,  this  figure 


with  the  1001  abundant  ^P,  by  the  mean  g  value  of 
2.007  and  by  the  dramatic  temperature  dependance  which 
is  related  to  motion  of  the  unpaired  spin  electron 
around  the  three  possible  silicon  bonds  next  to  the 
P-vacancy  pair.  At  a  suitable  low  temperature 
(approximately  180° K)  this  resonance  is  motion ally  narrowed 
and  the  E-center  resonance  becomes  much  larger  than  any 
other  resonances  making  the  E-center  observable. 

Although  an  exact  determination  of  this  defect  con¬ 
centration  has  not  been  made,  it  is  clear  that  an 
appreciable  fraction  of  all  NTD  phosphorus  (1.8  x  lO^cm 
is  bound  in  E-centers  after  irradiation  but  before 
annealing.  This  is  likely  to  be  due  to  the  P  recoil 
vacancies  which  get  immediately  trapped  in  this  center. 

The  room  temperature  EPR  spectrum  of  this  sample 

is  dominated  by  P-3  centers  which  are  linear  chains 

2  2 

of  four  vacancies.  This  seems  to  be  the  dominate 
room  temperature  defect  and  will  be  discussed  later. 

After  annealing  this  sample  to  450°C,  the  E.PR  spectrum 
is  dominated  by  the  P-1  five  vacancy  cluster  shown  in 
Fig.  f>3  for  the  magnetic  field  along  the  <lriQ>.  This 
P-1  center,  is  identified  by  the  g  value  at  2.009  and  the 
six  hyperfine  lines  with  “‘Si.  This  defect  is  also 
the  dominant  defect  in  neutron  damaged,  phosphorus 
doped  float  zone,  however,  it  is  generally  thought  to 
anneal  away  at  100°C  lower  temperatures  then  we  observe 
here.  This  P-1  center  appears  to  grow  in  N'T!''  silicon 
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on  annealing  while  the  P-3  center  decreases  at  around 

n  a 

200  C,  the  temperature  at  which  divacancies  anneal." 

It  is  likely,  therefore,  that  the  4-vacancy  P-3  centers 
capture  vacancies  from  dissociated  divacancies  and  col¬ 
lapse  into  the  5  vacancy  P-1  center.  After  annealing  to 
830°f,  the  only  1: PR  signal  to  he  seen  at  room  temperature 
(!  ig.  81)  is  a  characteristic  surface  signal  at  g  =  2.0035 
that  has  appeared  because  the  sample  was  purposely  crushed 
to  randomize  orientation.  On  cooling  to  100°K  the  free 
carrier  EPR  signal  at  g  =  1.999  becomes  visible  because 
reduced  scattering  at  low  temperature  has  sharpened  the 
resonance  as  expected.  Since  the  line  shape  of  this 
signal  is  known  from  theory,  it  is  straight  forward  to 
convert  this  intensity  to  relative  concentration.  It  is 
hoped  eventually  to  obtain  further  information  about 
h’  Nil)  production  from  this  resonance. 

We  now  return  to  an  investigation  of  the  FZ  and  CZ 

samples  (MURR-167  and  168)  which  had  been  irradiated  to 

IQ 

a  thermal  neutron  fluence  of  1.54  x  10  "n/cm".  The  CZ 
sample  of  this  pair  is  the  sample  in  which  the  con¬ 
centration  of  A-centers  was  determined  and  the  decrease 
in  oxygen  interstitials  observed.  Since  a  large  con¬ 
centration  of  vacancy- oxygen  centers  was  observed 
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fa,  i  x  10  cm  ')  which  accounted  for  17$  of  the  total 
concentration  of  oxygen  interstitials  originally  in  the 
CZ  sample,  one  might  assume  that  the  vacancies  for  these 
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figure  54.  Id’R  spectrum  of  surface  states  at  room 

temperature  and  conduction  electrons  from 
MU  phosphorus  at 


A-centers  grew  at  the  expense  of  other  vacancy  containing 

defects.  Watkins  has  reported  on  the  difficulty  of 

observing  E-centers  in  CZ  material  because  of  the 

?1 

scavanging  of  vacancies  by  oxygen.  This  hypothesis 

was  in  fact  checked  by  irradiating  FZ  and  CZ  simultaneously 

and  observing  divacancy  and  P-3  centers.  The  resultant 

P-3  center  resonances  for  the  samples  irradiated  to  a 

19  2 

thermal  fluence  1.39  x  10  n/cm  is  shown  in  Fig.  S3. 

The  P-3  four  vacancy  linear  chain  dominates  the 

EPR  spectra  of  both  the  CZ  and  FZ  samples  at  room 

23 

temperature  after  irradiation.  This  P-3  center  is 
characterized  by  a  pronounced  angular  dependence  and 
a  complicated  structure.  The  fine  structure  splittings 
(S  =  1)  for  a  [110]  field  orientation  are  very  clear 
and  show  the  basic  12  line  pattern;  8  along  [110]  and 
an  outer  three  (2  down  field  and  one  up  field].  Basically, 
there  are  six  inequivalent  orientations  for  this  defect, 
however,  with  II  along  [110]  there  are  only  four  observed 
since  two  become  degenerate  and  we  have  the  intensity 
pattern  (1,2, 2,1)  shown  at  the  top  of  the  figure.  Since 
each  has  S  =  1,  there  are  two  transitions  for  each 
orientation  accounting  for  the  12  line  pattern. 

An  investigation  of  Fig.  SS  indicates  that  the 
intensity  of  the  P-3  resonances  for  the  CZ  and  FZ 
samples  are  similar.  On  correcting  these  intensities  for 
differences  in  sample  mass,  it  is  found  that  the  P-3 
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Figure  55.  PPR  spectrum  of  P-5  four  vacancy  chain 
dominate  defect  concentration  in  room 
temperature  MI)  silicon  irradiated  to 
Mi,  =  1.54  x  10^n/cnr-  in  both  (’5  and 


intensities  are  identical  to  within 


II!  The  concen¬ 


tration  of  P-3  centers  has  been  estimated  from  F.PR 
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signal  strength  to  be  about  3  to  4  x  10  cm  ,  i.e., 
a  concentration  similar  to  the  concentration  of  A-centers 
in  the  CZ  sample.  Since  the  EPR  and  optical  absorption 
have  been  done  on  the  same  samples,  then  it  is  clear 
that  P-3  centers  have  not  lost  vacancies  to  A-centers 
as  expected.  This  is  surprising  and  must,  of  course, 
be  checked  at  lower  fluences;  however,  this  experiment 
suggests  that  the  P-3  centers  are  intrinsic  to  the 
fundamental  cascade  process  and  do  not  coagulate  from 
a  general  vacancy  gas  created  during  irradiation.  This 
is  obviously  an  important  conclusion  and  could  lead  to 
an  understanding  of  the  differences  in  reactor  neutron 
spectra  as  they  affect  the  NTD  lifetime  problem. 

What  about  other  vacancy  centers?  In  this  heavily 
irradiated  pair  of  samples  no  E-centers  are  found  in 
either  the  CZ  or  FZ.  Since  the  concentration  of  phosphorus 
in  these  samples  produced  by  NTD  must  be  about 
2  x  lO^cm  J  and  the  B  recoil  vacancy  concentration  a 
similar  number,  F-centers  can  not  contribute  enouj ' 
vacancies  to  account  for  the  observed  A-center  con¬ 
centration. 

Perhaps  the  divacancy  concentration  contributes  to 
A-center  formation  in  CZ  material.  We  run  into  a 
fundamental  experimental  difficulty  here,  as  shown  in 
Fig.  3f).  For  these  heavy  irradiations,  the  hand  edge 
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Figure  56.  Near  band  edge  absorption  for  CZ  sanmle 
irradiated  to  1.5-1  ,x  JO*-1  thermal  n/cr- 
at  room  temperature  and  at  "’“K.  The  l.S  uni 
di vacancy  band  is  lost  due  to  band  edge 
shift  in;’,  however,  the  5  mn  divacancy  bands 
now  appear. 
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shift  is  so  large  compared  to  a  sample  irradiated  to 
18  2 

1.51  x  10  n/cm  that  no  estimate  of  divacancy  con¬ 
centration  can  be  made  from  room  temperature  optical 
absorption  on  these  thick  samples.  On  cooling  to  77°K, 
however,  the  three  divacancy  bands  near  3  urn  are  seen 
for  the  first  time  in  our  samples.  It  is  perhaps 
possible  to  compare  the  3  pm  bands  in  the  CZ  and  FZ 

samples,  however,  the  3  pm  bands  are  known  to  be  very 

24 

Fermi  level  dependent.  It  is  therefore  not  clear 
that  the  procedure  is  useful. 

We  are  now  in  a  position  to  compare  the  displace¬ 
ment  theory  of  the  first  section  with  our  count  of 
vacancies  observed  experimentally.  We  have  made  a 

rough  estimate  of  divacancy  concentration  in  these 

2 

samples  from  data  obtained  by  the  ORNL  group.  These 

data  are  shown  in  Fig.  57  by  the  open  circles  while 

17  -  3 

the  cross  represents  our  estimate  of  4  x  10  cm  '  for 
the  approximate  divacancy  concentration. 

From  the  irradiation  time  and  using  the  displacement 

-T  J  *7 

rate  in  Table  3  of  2  x  101  di sp/cm  /sec ,  the  total 

1 9 

number  of  displacements  is  calculated  to  be  ^  4  x  10  cm 
The  number  of  vacancies  in  each  of  the  observed  defects 
is  listed  below. 


[V-V]  divacancies  =  4  x  lO^cm  ^ 

1  n 

[P-3]  4  vacancy  =  3  to  4  x  10  cm 
[I’-l]  5  vacancy  <  3  x  lO^cm  ^  = 


8  x  10 


17 


=  1 . 4  x 
17 


vacancies/cm' 

,n18  .  .3 

10  vacancies/cm 

vacanci  es/cm'^ 


104 


1.4  x  10 


io’8  1019 

Fluence  Above  1  Mev  (n/cm^) 

Estimate  of  concentration  of  di  vacancies 
sample  irradiated  to  1  .  x  10^  thermal 
neutron/cm1-  from  ORX'I,  data. 


[E]  vacancy- phosphorus  <  2.25  x  lO^cm  ^  =  2.25  x  10^ 
vacancies/cm"3 

1  ^  *  3  17  3 

[A]  vacancy- oxygen  =  1  x  20  cm  =  1  x  10  vacancies/cm 

1  8 

Total  vacancy  concentration  observed  =  2.44  x  10 
vacnnc  ics/ cm"’ 

]  q  -  3 

Total  calculated  displacements  =4x10  "cm 

Total  vacancies/total  displacements  =  —  -  -  — 

4  x  10iy 

or  6 . 1  %  . 

We  conclude  that  about  6%  of  the  total  displacements 
have  been  observed  in  these  experiments  by  optical 
absorption  and  EPR.  We  do  not  know  the  number  of 
displacements  which  have  been  lost  by  annealing  at 
room  temperature,  however,  61  remaining  is  not 
unreasonable.  It  would  be  of  interest  to  repeat  these 
same  experiments  at  lower  fluences  and  at  higher  flucnces 
to  see  if  effects  of  overlapping  cascades  could  lead  to 
a  saturation  of  the  defect  production  rate.  We  believe 
that  our  approach  is  sound  and  that  a  pursuit  of  this 
subject  can  yield  fundamental  knowledge  about  the 
displacement  cascades  and  recoil  events  which  is  not 
well  understood  at  this  time. 

I : .  Raman  Scattering  in  Neutron  Irradiated _S i 

In  the  previous  subsection,  a  defect  bookkeeping 
showed  that  only  about  6'1  of  the  total  displacements 


are  observed  as  vacancy  related  point  defects.  Where 

are  the  rest  of  the  displacements?  Raman  scattering 

has  suggested  that  most  of  the  remaining  displacements 

are  located  in  highly  disordered  material  which  appears 

in  many  ways  to  be  amorphous. 

It  is  well  known  that  irradiating  silicon  with 

neutrons  produces  damage  in  the  crystal.  Theoretical 
7  S 

models'"  predict  that  recoil  atoms  from  gamma  ray  or 
particle  emissions  either  due  to  the  capture  of  ther¬ 
mal  neutrons  or  due  to  the  elastically  and  inclast icall 
scattered  fast  neutrons  produce  clusters  of  atomic 
displacements  in  a  crystal.  At  a  sufficiently  high 
neutron  dose  the  clusters  appear  to  be  amorphous  - 1  ike 

in  nature,  as  evidenced  by  the  appearance  of  new  de- 

2  6 

feet  induced  vibrational  modes."  The  extent  of  amor- 

phization  depends,  however,  on  the  damage  density. 

A  damage  density  of  6  x  lO^eV/cm"'  at  80° K  is  shown 

to  produce  an  amorphous  layer  by  ion  implantation."' 

We  have  performed  experiments  with  samples  irradiated 

at  300°K  at  the  University  oT  Missouri  Research  Reactor 

facility  (MIJRR),  to  thermal  neutron  flucnces  of  3.5  to 
19  2 

7  x  10  ‘ n/cm  .  lor  silicon  the  energy  available  for 

atomic  displacements  is  principally  due  to  Tast 

neutrons,  in  light  water  reactors."'1  The  fluence  of 

fast  neutrons  (energy  >  188  eV l  in  our  experiments  is 

1 ()  ’ 

therefore  0.7  to  1.3  x  1 0  n/cn". 


from  t  lie  reactor 


spectrum  at  MURR  we  have  estimated  the  mean  Si  recoil 

energy  to  be  ^  22  keV.  For  the  highest  fluence  the 

concentration  of  displacements  is  calculated  to  be 
20  -  3 

3.6  x  10  cm  ,  corresponding  to  a  cascade  density 

1  O  _  7 

of  2  x  10  cm  .  Consequently,  0.71  of  the  crystal 

volume  is  damaged.  Stein^  has  estimated  that  the 

21  x 

damage  density  for  a  22  keV  recoil  is  ^  10  eV/cm  . 

The  ratio  of  this  damage  density  to  that  required  to 

produce  totally  amorphous  material  is  0.6®,  which 

is  in  good  agreement  with  the  damaged  fraction  of  the 

crystal  volume  we  obtain  for  our  reactor  spectrum. 

We  irradiated  a  number  of  Si  samples  to  a  fluence 
19  1 

of  7  x  10  ntj1/cm“.  Before  and  after  irradiation, 
the  resistivity  of  each  sample  was  measured  (Table  11). 

We  have  begun,  to  our  knowledge,  the  first 
studies  of  defect  induced  vibrational  modes  in  neutron 
irradiated  Si  using  Raman  Scattering.  We  observe 
features  in  the  Raman  spectra  due  to  first  order 
processes  that  arise  from  the  destruction  of  transla¬ 
tional  invariance  in  the  lattice  with  peaks  that 
correspond  to  maxima  in  the  vibrational  density  of 
states  (VDOS)  averaged  over  the  Brillouin  Zone  (RZ) . 
These  new  peaks  are  very  similar  to  those  seen  in  the 
Raman  spectra  of  amorphous  Si''  ,  and  we  believe  that 
they  arise  from  amorphous  - 1 i kc  zones  that  are  formed 
in  the  crystal  due  to  irradiation.  Some  of  the  first 
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Table  11.  Resistivity  of  Raman  samples  before  aiid 
after  irradiation  (<t>  =  7  x  1 0 ■*  --n^/ cm- )  . 


Resistivity  (ft-cm)  and  type 


Sample 

Number 

Before 

irradiation 

After 

irradiation 

HD  #  1 

2  x  10'2  (p) 

3.5  x  103 

HD  #  3 

4  x  10'3  (p) 

6.8  x  103 

HD  H  4 

5  x  10'*  (p) 

1.1  x  in4 

HD  #  6 

4  x  10  3  (n) 

1  x  103 

HD  ft  7 

2  x  10'3  (n) 

2.0 

HD  ft  8 

9  x  10'3  (n) 

2.9  x  104 

HD  ft  9 

1  x  10'3  (p) 

8.6  x  102 

CZ  ft  2.1c 

3.4  x  102  (n) 

1.6  x  in6 

FZ  It  2.1 

2.5  x  102  (n) 

1.9  x  10° 

order  processes  that  we  observe  have  been  observed 
before  in  infrared  absorption  spectra.  However, 

Raman  spectroscopy  offers  the  advantage  of  observing 
both  crystalline  and  amorphous  - 1 i ke  features  in  the 
same  spectrum,  while  only  the  amorphous  - 1 i ke  features 
are  observable  in  the  IR  due  to  the  selection  rules. 

The  Raman  scattering  spectra  were  obtained  at 

O  + 

300°K  using  the  6471  A  radiation  from  a  Kr  laser, 
a  double  monochromator  with  holographic  gratings,  and 
a  cooled  photomultiplier  with  photon  counting  clcctroni 
figure  58  shows  the  Raman  spectrum  of  Czochralski  Si 


Neutron  Irradiated 
6471A,  HH+HV 

- Before  Irradiation  I 
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Figure  5S.  Stokes  Raman  spectra  of  Si  before  and  after 
neutron  irradiation.  The  insert  shows  a 
schematic  of  the  BI  for  a  crystal  with  2,  4 
and  In  atoms  per  unit  cell. 


(340  2- cm  before  irradiation,  sample  #  CZ  2.1c) 

before  and  after  irradiation  to  a  fluence  of 
10  2 

7  x  10  " n/cm  .  Before  irradiation,  the  Raman  spec¬ 
trum  of  Si  is  dominated  by  the  only  allowed  first 
order  process,  the  zone  center  optic  phonon  at  a 
frequency  of  520  cm  ^ .  Other  weaker  modes  are  over¬ 
tone  processes  such  as  the  second  order  transverse 
acoustic  mode,  2TA(X)  at  300  cm  *.  Upon  irradiation, 
two  new  peaks  appear,  at  w.  05  and  400  cm  *,  superim¬ 
posed  on  an  enhanced  background  scattering  below  the 
frequency  of  the  optic  phonon.  The  new  features 
appear  in  a  Tj  scattering  geometry  and  are  first 
order  in  nature.  They  are  observable  both  in  the 
float  zone  (25(1  T-crn,  sample  >'  l-'Z  2.1)  and  boron 

10  Of)  7 

doped  (6  x  10  -  1  x  10"  holes/cm'  ,  samples  H  111)  1, 

3,  4  and  0)  silicon.  The  intensity  of  the  peaks 
scales  with  the  neutron  dose,  as  we  find  from  a  com¬ 
parison  of  Raman  spectra  taken  for  sample  «  CZ  2.1c 

and  sample  CZ  i‘  10'  1,  which  was  irradiated  to  a  fluence 

]  q  •> 

ot  3.6  x  10  '  nt|,/cin“.  In  the  insert  of  Fig.  58  we 
show  schematically  the  origin  of  these  new  peaks. 

In  fully  crystalline  Si,  which  has  two  atoms  per 
unit  cell,  the  BZ  would  have  acoustic  (LA  and  TA )  and 
optic  (TO  and  TO)  branches  as  in  panel  fa),  Fig.  58. 

The  optic  mode  (520  cm  '  in  Si)  at  the  zone  center 
would  be  the  only  allowed  first  order  process  f£  -  0) 


which  is  Raman  active.  If,  for  instance,  a  defect 
was  introduced  on  the  average  in  every  alternate  cell 
so  that  the  new  unit  cell  is  double  the  original,  the 
new  BZ  is  halved,  and  is  equivalent  to  folding  the 
BZ  of  (a)  back  into  the  zone  center  once  ("panel  (b)  , 

Fig.  58).  New  Raman  active  modes  now  appear  at  the 
zone  center.  By  extending  this  argument,  in  fully 
amorphous  Si,  where  translational  invariance  is  com¬ 
pletely  lost,  the  entire  BZ  is  folded  into  the  zone 
center.  The  narrow  zone  center  optic  phonon  mode  at 
520  cm  1  disappears,  and  the  spectrum  obtained  con¬ 
sists  of  broad  peaks  representative  of  the  vibra¬ 
tional  density  of  states  averaged  over  the  entire  BZ. 
This,  schematically,  is  what  is  observed  in  amorphous 
Si.  In  our  samples  we  see  both  crystalline  and  amor¬ 
phous-like  features,  which  we  believe  is  due  to  the 
fact  that  we  have  amorphous- 1  ike  zones  in  a  mostly 
crystalline  matrix. 

The  comparison  between  the  spectra  obtained  in 
our  samples  and  in  amorphous  Si  is  shown  in  Fig.  59. 

We  plot  the  reduced  Raman  spectrum  of  amorphous  and 

neutron  irradiated  Si  (crystalline  spectrum  subtracted). 

2  9 

As  defined  by  Shukcr  and  Gamon  ,  the  reduced  Raman 
intensity  I^(w) 
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deconvolutes  from  the  measured  Raman  intensity 
the  contributions  arising  from  the  harmonic  oscilla¬ 
tor  part,  Bose  Einstein  distributions  and  dependence 
on  exciting  frequency  The  resulting  spectrum 

I R (to)  contains  only  effects  due  to  matrix  elements, 
which  are  slowly  varying,  and  more  importantly,  the 
vibrational  density  of  states  (VDOS) .  The  crystalline 
VDOS  obtained  from  a  shell  model  fit  to  neutron  scat¬ 
tering  data  is  also  shown  in  Fig.  59  (dashed  lines). 

We  note  that  the  broad  peaks  in  neutron  irradiated 
Si  do  indeed  correspond  in  frequency  to  peaks  in  the 
VDOS.  We  therefore  attribute  the  low  frequency  peak 
(the  95  cm  1  peak  in  I^(w)  is  the  125  cm  1  peak  l^(w) 
due  to  the  form  of  the  deconvolut ing  function)  as 
arising  from  averaged  transverse  acoustic  (TA)  modes, 
the  central  background  from  longitudinal  acoustic 
(LA)  modes  and  the  high  frequency  (490  cm  ^)  peak 

from  optic  modes  LO  and  TO).  This  assignment  is 

?  8 

similar  to  that  made  in  amorphous  Si. 

An  interesting  feature  of  the  amorphous  zone  for¬ 
mation  by  neutron  irradiated  Si  is  that  the  zones 
arc  distributed  uniformly  through  the  bn  1 k  of  the 
sample.  Neutrons  have  a  penetration  depth  of  many 
centimeters  in  Si,  and  the  damage  and  amorphous  zones 
produced  arc  truly  within  the  bulk  and  not  confined 
to  surface  layers  of  the  order  of  a  micron  in  thkk- 

1  -  ' 


ness  as  in  the  case  of  ion  implanted  Si.  At  a  suffi¬ 
ciently  high  neutron  dose,  therefore,  it  should  be 
possible  to  irradiate  Si  into  an  amorphous  state  and 
make  bulk  amorphous  Si. 

Much  of  the  damage  produced  by  irradiation  can 
be  annealed.  We  have  studied  the  Raman  spectrum  as 
a  function  of  isochronal  annealing  for  one  sample. 

Our  results  arc  shown  in  Fig.  60,  for  sample 
#  CZ  2.1c.  The  anneals  were  performed  in  a  furnace 
for  15  minutes  at  each  temperature  in  an  atmosphere 
of  flowing  nitrogen  gas.  The  structures  induced  by 
irradiation  decrease  in  intensity  (Fig.  01)  with  in¬ 
creasing  anneal  temperature  and  anneal  out  completely 
at  600°C,  when  the  Raman  spectrum  is  identical  to 
that  of  non - i rrad iated  Si.  While  most  of  the  annealing 
activity  takes  place  around  250°C,  the  peaks  linger 
on  and  disappear  into  the  noise  only  around  550  to 
600°C.  This  kind  of  annealing  behavior  is  typical 
of  amorphous  Si  made  by  other  techniuues. 

Tn  addition  to  a  decrease  in  intensity,  the  low 
frequency  peak  shifts  to  higher  energies  as  we  anneal 
the  sample  (Fig.  62).  The  peak,  which  is  at  125  cm’1 
in  the  reduced  Raman  spectrum  prior  to  anneal  (i.c., 

30°C  anneal)  shifts  up  to  150  cm’1  before  it  disappears. 
W'e  attribute  this  to  the  fact  that  as  we  approach 
crystallinity,  contributions  at  the  zone  center  arise 
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Figure  ()1.  'Hie  intensity  of  the  amorphous  - 1  i  kc  TA  and 
optic  peaks  produced  by  neutron  irradiation 
as  a  function  of  anneal  temperature. 
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62.  The  frequency  shifts  of  the  Inch  and  lou  fre¬ 
quency  amorphous  - 1 i kc  Raman  peaks  as  a  function 
of  anneal  temperature. 
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mainly  from  the  edges  of  the  BZ,  corresponding  to  a 
few  folds  in  the  BZ,  which  would  tend  to  shift  the 
peak  to  higher  frequencies.  In  other  words,  the 
peak  is  similar  to  the  overtone  2TA  structures  in 
crystalline  Si  (c-Si)  but  is  at  half  the  frequency, 
being  a  folded  back  single  phonon  mode.  The  main 
contribution  to  the  2TA  structures  in  crystalline 
Si  arises  from  the  2TA(X)  mode  at  300  cm  ^ .  The 
single  phonon  mode  due  to  a  folded  TA(X)  branch  would 
yield  a  frequency  of  150  cm  1 ,  a  number  in  good  agree¬ 
ment  with  the  150  cm  1  peak  we  observe  after  a  large 
degree  of  annealing  lias  taken  place.  The  high  fre¬ 
quency  (490  cm  1)  peak  does  not  shift  with  anneal  tempera¬ 
ture,  as  far  as  we  can  tell  within  our  present  signal-to- 
noise  levels.  Since  the  optic  branch  has  considerably 
less  dispersion  than  the  acoustic  branches,  the  smaller 
number  of  BZ  folds  upon  annealing  arc  expected  to  pro¬ 
duce  only  slight  shifts  in  the  peak  frequency.  Detecting 
these  small  changes  will  require  a  more  detailed  study. 

The  Raman  spectrum  of  intrinsic  silicon  is  charac¬ 
terized  by  an  intense  zone  center  optic  phonon,  at  a 
frequency  of  520  cm  * ,  and  other  weaker  two  phonon 
overtone  bands.  Donor  for  acceptor)  levels  close  to 
the  conduction  (.valence)  bands  are  produced  by  doping 

with  impurities.  Discrete  electronic  Raman  transi* 

3d  -  52 

tions'  '  between  even  parity  states  of  the  impurity 
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levels  can  be  observed  at  impurity  concentrations  of 

1016  to  10^  cm*\  At  heavier  doping  levels  O  10^ 

cm  ) ,  the  discrete  impurity  states  broaden  into  an 

impurity  band,  the  material  is  rendered  degenerate, 

and  the  Fermi  level  moves  into  the  conduction  (or 

valance)  bands  producing  a  free  carrier  solid  state 

plasma.  In  polar  semicondr  tors  such  as  GaAs  ,  GaP, 

etc.,  the  dominant  effect  of  the  free  carrier  plasma 

is  the  mixing  of  the  plasmon  with  the  longitudinal 
33 

optic  phonon.  In  a  heavily  doped,  non-polar  semi¬ 
conductor  such  as  silicon,  the  large  number  of  free 
carriers  participate  in  low  frequency  inter-  and  intra¬ 
band  transitions  that  arc  Raman  active,  and  can  be  ob¬ 
served  by  visible  laser  Raman  spectroscopy. 

Interband  transitions^4  that  occur  between 
neighboring  bands  (  Tj  -+  Fj<  conduction  bands  or  the 
two  uppermost  Tg  valence  bands)  form  an  electronic 
continuum  that  overlaps  the  frequency  of  the  zone 

center  optic  phonon  for  carrier  concentrations  of 
19  -  3 

10  cm  and  higher.  These  transitions  also  have 
the  same  symmetry  as  that  of  the  optic  phonon,  con¬ 
sequently  a  discrete-continuum  interaction  takes 
place,  producing  a  typically  asymmetric,  broadened 
and  shifted  Fano  type  lineshape  for  the  optic  phonon 
described  by  the  expression 

7  7 

I  (u)  (q  *■  ()““/)  +  t  “ 


1  SO 


where 


e  =  (w-Wg  -  Aw) /T 

and 

q  =  (— E—2-  +  Aw  )/r 
‘e  ° 

The  frequency  shift  AwQ ,  from  the  f requcncy  w0  in  in¬ 
trinsic  Si,  and  broadening  represent  the  real  and 
imaginary  parts,  respectively,  of  the  self  energy 
of  the  phonon  due  to  a  deformation  potential  inter¬ 
action  with  the  electronic  continuum. ^  The  asymmetry 
parameter  q  reflects  the  signs  and  relative  strengths 
of  the  optic  phonon  T  ,  electronic  continuum  Tg  and 
the  electron  phonon  interaction  through  the  res¬ 
pective  matrix  elements.  We  have  measured  and  ealeu- 
lated  the  three  Fano  parameters'  ' ‘  for  n-Si  and 
find  excellent  agreement  between  theory  and  experi¬ 
ment.  The  more  complex  valence  bands  involved  in 

T  O 

p-Si  make  calculations  somewhat  difficult.' 

Intraband  Transitions'^  “  ’  *  ''  11  m  an  i  f  cs  t  them¬ 

selves  in  the  form  of  a  low  freuuency  tail  in  the 
Raman  spectrum,  extending  from  near  zero  to  about 
300  cm  '  in  both  n-  and  p-Si.  The  tail  is  attri¬ 
buted  to  intervallev  fluctuations  arising  from  the 
six  ellipsoidal  conduction  band  valleys''  '  ^  in  n-Si. 
In  p-Si  we  have  found  that  the  complex  anisotropic 
valence  bands  ’  can  be  well  approximated  by  eight 


ellipsoid  < 1 1 0 >  "valleys."  Such  a  model  accounts  for 
the  experimentally  observed  selection  rules  and  other 
stress  and  frequency  dependent  behavior. 

In  addition  to  free  carrier  effects,  I ocal i zed 
vibrational  modes  associated  with  light  dopant  atoms 
appear  at  frequencies  higher  than  that  of  the  host 
optic  phonon. These  local  modes  also  interact  with 
the  interband  continuum  and  consequently  take  on  an 
asymmetric  lineshape  similar  to  that  of  the  optic 
phonon.^3  ^ 

These  three  major  features  that  appear  in  the 
Raman  spectra  of  heavily  doped  Si  arc  all  affected  by 
irradiation.  The  most  dramatic  change  appears  in  the 
optic  phonon.  The  removal  of  free  carriers  by  radia¬ 
tion  induced  damage  removes  the  interband  electronic 
continuum  responsible  for  the  Tano  type  lineshape, 
and  the  optic  phonon  becomes  narrow  and  symmetric  as 
in  intrinsic  Si.  The  intraband  tail  is  absent  (Tig. 

6  3 ,  samp  1 e  HD  %  9 ) . 

Furthermore,  the  local  modes  of  boron  disappear 
in  irradiated  Si(B)  suggesting^  ^  that  boron  is 
removed  from  its  usual  substitutional  site  and  is 
probably  also  surrounded  by  defects  that  prevent  it 
from  vibrating  at  its  normal  local  mode  frequency. 

An  investigation  of  this  boron  local  mode  intensity 
decrease  and  the  increase  in  t  lu  amorphous  Raman 
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Figure  (o.  The  Raman  spcct  nun  in  t  lie  region  of  the  optic 
phonon  in  sample  HP  *  !'  before  and  al  ter  i  rra 
tiiation.  Note  the  marked  asymmetry  in  the 
lineshape  doe  to  the  interaction  with  flee 
carriers,  which  disappears  alter  irradiation. 


peaks,  can  possibly  yield,  for  the  first  time,  a 
direct  measure  of  the  size  of  the  neutron  produced 
displacement  cascades.  This  information  is  of  con¬ 
siderable  interest  from  the  point  of  view  of  funda¬ 
mental  radiation  damage  theory  because  direct  experi¬ 
mental  evidence  of  cascade  size  and  overlap  is,  to 
date,  very  meager. 

The  free  carriers  removed  by  irradiation  can  be 
recovered,  however,  in  a  controlled  fashion  by  con¬ 
ventional  or  laser  annealing.  The  concentration  of 
free  carriers  after  each  anneal  is  determined  by  elec¬ 
trical  measurements.  The  inter-  and  intraband  transi¬ 
tions,  and  the  interactions  with  the  phonons  and  local 
modes  as  a  function  of  carrier  concentration  are 
presently  being  studied. 

The  resistivity  of  the  irradiated  sample  CZ  2.1c 
was  also  studied  as  a  function  of  anneal.  The  measure¬ 
ments  were  done  using  the  four  point  probe  technique. 

A  plot  of  resistivity  as  a  function  of  anneal 
temperature  is  shown  in  Fig.  6l.  Upon  irradiation, 
the  resistivity  of  the  sample  increased  from  310  o-cm 
to  'v  10&  ft-cm  due  to  the  well  known  trapping  of  free 
carriers  by  the  defects  induced  by  irradiation.  The 
resistivity  stayed  in  the  10b  n-cm  region  for  anneal 
temperatures  <  bOO°C,  then  dropped  to  1.5  !:-tm  after 
the  700°C  anneal,  where  most  of  the  defects  are  believe^ 
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to  have  annealed  out.  This  temperature  agrees  well 
with  the  disappearance  of  the  amorphous - i i ke  peaks 
in  the  Raman  spectrum.  Based  on  the  fiuence  and  the 
production  rate  of  phosphorus  due  to  transmutation  of 
^Si,  tiiis  is  close  to  the  expected  final  resistivity 
of  '  0.7  fi-cm.  Since  we  used  a  highly  polished  and 
etched  surface  (that  was  made  for  the  Raman  Scattering 
experiments)  for  the  resistivity  measurements,  we  do 
not  know  for  sure  if  the  small  peaks  in  the  resistivity 
(e.g.,  the  peak  at  400 °C!  anneal)  arc  related  to  the 
annealing  of  specific  defects. 

Cl .  PITS  in  Neutron  Irradiated  Si 

In  this  final  subsection  we  will  attempt  to  answer 
the  primary  question,  "How  do  the  defects  observed  in 
the  previous  experiments  behave  electrically?"  The 
PI.l'S  and  TCS  techniques  are  ideal  experiments  to  deter¬ 
mine  defect  energy  levels  in  the  hand  gap  and  investi¬ 
gate  their  annealing.  Because  of  the  high  sensitivity 
of  this  technique,  however,  it  should  be  noted  that 
the  very  low  concent ra t i ons  of  defects  in  this  experi¬ 
ment  discourage  defect  interactions  which  are  expected 
at  the  much  higher  fluences  used  in  the  previous  fl'R 
and  optical  experiments. 

In  order  to  understand  the  implications  of  carrier 
capture  and  thermal  emission  from  deep  defect  levels 
in  the  depletion  region  of  a  reverse  biased  diode,  it 


is  first  necessary  to  explain  these  kinetics  in  a 
normal  semiconductor.  The  advantages  of  using  a 
diode  depiction  layer  will  then  become  obvious. 

We  will  begin  by  describing  the  carrier  kinetics 
of  trap  emptying  which  lead  to  the  conventional  equa¬ 
tions  for  junction  capacitance,  depletion  width  and 
junction  current  as  a  function  of  time.  Ke  will  then 
derive  signal  processing  equations  for  single  and 
double  gate  boxcars  and  lock-in  amplifiers.  finally 
we  will  present  our  results  for  a  lock-in  amplifier 
detection  system  on  Si  diodes  which  have  been  neutron 
irradiated. 

Figure  65  represents  schematically  the  four 
possible  processes  for  filling  and  emptying  traps 
thermally.  We  will  neglect  optical  and  Auger  processes 
The  analysis  presented  is  taken  from  Ref.  50  which 
appears  to  be  the  definitive  paper  on  the  subject. 

Other  useful  papers  are  Refs.  51,  52,  5"  and  54. 

The  rate  equations  for  the  four  processes  a,  b, 
c,  d  in  l  ip,.  (-.5  are: 

-dn/dt  a  -  b  =  cn  npT  -  ennT  I'M! 

-dp/dt  =  c  -  d  =  cp  pnT  -  eppT  ("21 


dn  j/dt  - 

(a  b)  -  (c  - 

d)  = 

(cnnFT  -  cnnT^ 

- 

(CppnT  -  eppT) 

!cnn  >  e.^p.j  - 

fcpP 

+  °n1nT  1  " 

nT  =  PT  * 

n'l  ' 

"  1 

I 


The  parameters  cn,  cp  are  the  electron  and  hole  cap¬ 
ture  coefficients  while  en,  ep  are  the  thermal  emission 
rates.  The  equation  for  dpy/dt  can  be  obtained  from 
Eq.  (73)  and  the  constraint,  Eq.  (74).  Eliminating  Py 
in  Eq.  (73)  by  using  Eq .  (74)  and  realizing  that  in 
the  depletion  width  of  a  diode,  n  =  p  =  0,  we  obtain 
the  important  equation: 

dnT/dt  =  -(en  +  ep)nT  +  epNT  .  (75) 

This  equ  ion  and  Eq.  (74)  for  py  =  Ny  -  ny  are  the 
basis  of  DLTS  as  well  as  TCS  theory. 

These  kinetics  equations  can  be  solved  once  ini¬ 
tial  conditions  are  defined.  We  have  the  following 
set  of  initial  conditions  for  various  experimental 
situations  assuming  the  bias  pulse  either  completely 
empties  or  fills  the  trap  of  interest. 

I.  N-type  Depletion  Region 


A. 

Majority  carrier  trap: 

t  =  o, 

nT  (0)  =  N y 

B. 

Minority  carrier  trap: 

t  =  0, 

Py ( 0 )  =  0 

P- 

type  Depletion  Region 

A. 

Majority  carrier  trap: 

o 

II 

■M 

Py(0)  =  Ny 

B. 

Minority  carrier  trap: 

O 

II 

4-> 

ny(0)  =  0 

These  initial  conditions  can  be  considered  as 
working  definitions  of  majority  and  minority  traps  in 
n-  and  p-type  material.  If  the  initial  conditions  are 
not  entirely  satisfied,  then  the  traps  are  said  to  be 
unsaturated  with  carriers. 
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The  solutions  to  these  differentia]  equations  are 


now  straightforward.  For  case  I. -A,  we  find  that 


nl(t:)  =  N't 
and  using  Eq.  (74), 


en  +  eP 


e’  ^en  +  ep^)tj 


Pq’(t)  =  Nj  -  n-p  - 


nT  ei 


e  -i-  g 

n  cp 


(1  -  e 


(en  +  en)^ 


)  .  I 


I 


Exactly  the  same  equations  exist  for  the  other 
three  cases  with  a  suitable  substitution  of  the  appro¬ 
priate  carrier  concentration  and  emission  rate  sub¬ 
scripts.  These  results  are  summarized  below: 

N-TYPE  -  MAJORITY  TRAP 
Use  Eqs.  (76)  and  (77) 

N-TYPE  -  MINORITY  TRAP 

Use  Eqs.  (76)  and  (77)  with  the  substitution 
nT(t)  - 9-  pT(t) 

Pj  (  t  )  nq-  (t  ) 


II 


III. 


P  n 

i.e. ,  nT(t) 


ep  NT 


en  +  ep 


(1  -  e 


* (en  +  ep) I 


P-TYPE  -  MAJORITY  TRAP 


Same  Equations  as  II. 

IV.  P-TYPF.  -  MINORITY  TRAP 
Same  Equations  as  I. 

Equations  (76)  and  (7”),  therefore  represent  all 
possible  situations  with  appropriate  interpretation. 


1  P(l 


To  determine  the  temperature  dependence  of  the 


thermal  emission  rate,  we  define  the  probability  that 
a  trap  is  occupied,  using  Fermi  statistics  as 

P(Et)  =  nT/NT  =  [ l  +  g'1  e  (ET  "  V^V1.  C78) 

In  equilibrium,  the  total  emission  rate  must  equal 
the  total  capture  rate,  therefore 

emission  rate  =  enP(ET) 
capture  rate  =  cn [ 1  -  P(E-j')] 


or 


en  P  =  cnU  '  Pi- 

Substituting  Eq.  (78)  into  this  last  equation  we 
obtain 


e 


n 


-n 


g-l  e(ET  - 


But  defining  the  capture  coefficient  in  terms  of 
the  capture  cross  section,  cn  =  n  <vn>  on  = 

Nc  e  ~^EC  '  Ep)/kT  <yn>  CTn-  therefore  (assuming  g=l), 


e„  -  nc  on  <vn>  -  BtJ/W 


r(.n 


and  by  similar  arguments 


ep  -  Nv  op  <vp>  e'CFT  '  Evi/kT 


lb  l 


(  80  ) 


2 


The  thermal  carrier  velocities  are  given  by  <vn> 

=  3kT/me*,  etc. 

We  shall  now  consider,  for  simplicity,  only 
majority  carriers  in  n-type  material.  The  equations 
can  easily  be  generalized  to  minority  carriers  and 
to  p-type  material. 

The  actual  diode  bias  scheme  used  in  a  p+n  DLTS 
experiment  is  shown  for  a  diode  at  the  top  of  Fig.  66. 
The  diode  depletion  width  is  modulated  by  a  rapid 
change  in  the  diode  bias  as  shown  at  the  bottom  of 
this  figure  while  the  actual  junction  capacitance  is 
shown  schematically  at  the  center  of  Fig.  66. 

The  diode  is  originally  in  a  reverse  biased  con¬ 
dition  (T)  so  that  a  depletion  region  is  formed  and 
the  traps  have  had  time  to  empty.  At  ( 2 )  ,  the  reverse 

bias  voltage  is  collapsed  which  reduces  the  depletion 
region.  Because  electrons  rapidly  rush  into  the  con¬ 
duction  band  of  the  depleted  region,  the  traps  will 
now  tend  to  fill  very  rapidly  because  the  concentra¬ 
tion,  n,  of  electrons  is  large  in  Eq.  (71).  The 
junction  capacitance  rises  very  rapidly  at  (2)  as  a 
result  of  the  collapse  of  the  depletion  region  since 
Cft)  =  cA/W(t).  We  begin  our  capacitance  transient 
at  Q)  .  We  have  designated  this  time  t  =  0.  Because 
extra  charge  is  left  in  the  depletion  region  which 
neutralizes  the  charge  of  the  ionized  donors,  the 
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junction  capacitance  is  lower  at  t  =  0  by  an  amount 


AC  than  the  capacitance,  C  ,  before  the  bias  pulse. 

We  shall  show  that  the  amplitude  of  this  capacitance 
transient,  AC,  is  proportional  to  the  trap  concentra¬ 
tion,  M.j. .  During  the  rest  of  the  period,  P,  in  region 
(^) ,  the  trap  carrier  concentration  decays  toward  CQ 
with  a  time  constant,  t,  which  we  will  now  show  to  be 
the  reciprocal  of  the  carrier  thermal  emission  rate. 

The  time  dependence  of  the  junction  capacitance 
is  given  by  the  usual  equation  for  a  junction  of 
area  A, 


C(t)  = 


eA 

wTtj 


(81) 


where  the  depletion  width,  W(t),  is  given  for  a  one 
sided  abrupt  junction  as 


l  /  2e  (VBi  +  V) 

K(t>  -V  oNfCtl- 


(82) 


and  where  Nj(t)  is  given  for  majority  carrier  traps 
in  n-type  by  *  n^. (t).  Thus, 


C (t )  =  A  |/q  e  (Nr  -  nT(t)) 

y  2(vBi  +  v) 


This  last  equation  can  be  written  as 


c(t)  -  A 

nx 

s  Co  fl  -  Ifrf 


o  *  ktt) 


nT^ 


r.'  i 


A 


u 

where  Co  =  A[qe  Np/2(vB^  +  V)]  is  the  junction  capa¬ 
citance  an  infinite  time  following  the  bias  transient. 
Tli is  can  also  be  written 

a£  -  c°  ~  c:(tl  -  1!  (tm 

C0  -  c0  -  :n„  •  (R-'1 

Using  fq.  (70),  and  assuming  that  en^->Cp  for  majority 
traps  in  n-type,  which  follows  for  a  level  in  the 
upper  half  of  the  band  gap  from  bus.  (79)  and  (SO), 

n-j-(t)  =  Nje  where  1  / r  ~  en 


o  r 


AC 
C  . 


c-t/> 

2Nn 


(  S 1  I 


I'll  is  is  the  basic  equation  for  the  1U.TS  signal.  Since 
AC/C0  a.  10  for  the  Boonton  bridge,  ;i  doping  level 
\p  =  10^('cm  1  will  yield  a  trap  concentration  Nj-  a> 

10  “cm  lowering  \’p  will  increase  this  sensitivity. 

Notice  that  this  equation  at  t  =  0  determines  the 
trap  concentration 


Instead  of  measuring  junction  capacitance,  it  is 
possible  to  measure  junction  current.  Ibis  transient 
current  spectroscopy  (TUN)  is  useful  for  traps  with 
high  emission  rates  or  for  concentration  ('ref  i  liny- 


near  the  inaction. 


The  total  junction  current  due  to  trap  emission 
is  made  up  of  three  parts,  an  electron  trap  emission 
current,  a  hole  trap  emission  current  and  a  Maxwell 
Displacement  current  given  by^ 


ip(t) 

id(tj 


r 

A  /  q  (dn/dt)c!x 

J0 

/W 

A  /  q  (dp/dt } dx 

y0 

,  w 

A  J  q  (dnT/dt)(J)dx. 


The  total  junction  current  is,  therefore,  the  sum  of 
these  three  terms.  Using  Eqs.  (71)  -  (74)  it  can  be 
shown  that 


if1)  =  xn  +  ip  +  id  =  (ennT  +  ep?T^ •  (85) 


This  is  the  basic  equation  of  TCS  theory  and  is 
probably  sound,  however,  some  minor  approximations 
have  been  included  by  assuming  a  one-sided  abrupt 
junction,  neglecting  the  Debye  length,  etc. 

Again,  assuming  majority  traps  in  n-type  material 
and  using  Eqs.  (76)  and  (77),  wc  find  that 


)  ,  1/t 


(87) 


This  equation  contains  a  steady,  time  independent 


minority  leakage  current  proportional  to  e^  and  a  time 


dependent  majority  transient  ene 


t/T 


In  the  analysis 


of  most  TCS  experiments,''’"’  '’<)  it  is  assumed  that  the 
first  term  can  be  neglected  relative  to  the  second. 
This  is  clearly  not  true  for  an  infinite  time  after 
the  transient  since  the  second  term  vanishes.  Never¬ 
theless,  Tq.  (87)  is  usually  written'’^ 


i  -t/T  1  , 

2  t  »  T  “  n 


(88) 


Note  that  this  TCS  signal  equation  is  similar  to  the 
DLTS  signal  hq .  (84)  except  for  the  1/t  term.  Since 
1 / t  is  the  thermal  emission  rate,  the  TCS  signal  inten¬ 
sity  is  proportional  to  this  rate.  This  justifies  our 
previous  statement  that  TCS  can  be  used  to  advantage 
for  those  traps  with  high  emission  rates. 

We  shall  now  do  a  signal  analysis  of  DLTS  and 
TCS  signals  (laps.  (84)  and  (88))  using  both  boxcar  and 
lock-in  amplifier  detection.  A  signal  processing  dia¬ 
gram  is  shown  in  Tig.  (i 7  for  a  DLTS  lock-in  amplifier 
experiment.  The  signal  waveforms  arc  shown  on  this 
diagram.  If  a  boxcar  amplifier  is  used,  this  diagram 
is  essentially  the  same,  however,  the  form  of  the 
reference  signal,  er(t),  is  different. 

The  timing  diagram  used  in  the  following  analysis 
and  experiments  is  shown  in  Tig.  (>8.  The  reference 
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Figure  68. 


Timing  diapriin  for  DI.TS  cxncr  imcnt 


signal  for  a  double  boxcar  amplifier  is  shown  at  the 

bottom.  A  single  boxcar  reference  signal  would  be 

identical  except  that  the  second  gate-on  at  t?  would 

be  missing.  Note  that  either  the  lock-in  or  double 

boxcar  references  will  remove  any  D.C.  component  to  the 

signal  while  a  single  boxcar  will  not.  According  to 

Eq.  (88),  TCS  experiments  could  make  use  of  a  single 

boxcar  amplifier  provided  that  the  leakage  term  ep  is 
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truly  negligible  in  Eq.  (87).  If  a  D.C.  component 
exists  such  as  C0  in  a  DLTS  experiment,  then  the 
slight  temperature  dependence  of  this  D.C.  component 
will  cause  the  baseline  of  the  spectra  to  shift  as  a 
function  of  temperature.  For  a  diode  doped  on  the  low 
side  to  101()cm  \  carrier  freczeout  from  300°K  to  7?°K 
causes  a  20 change  in  C  . 

We  will  begin  by  analyzing  the  DLTS  signal  using 
lock-in  detection.  As  shown  in  Fig.  67,  the  capaci¬ 
tance  meter  outputs  a  voltage  proportional  to  the 
diode  capacitance.  The  signal  voltage,  es  =  Ac  +  K, 
where  A  is  an  amplitude  equal  to  CQNj/2Np  times  the 
capacitance  to  voltage  conversion  factor  of  the  Boonton 
bridge,  and  B  is  proportional  to  the  D.C.  component 
of  the  signal,  CQ .  From  Fig.  67,  the  demodulator  mul¬ 
tiplies  this  signal  times  the  square  wave  reference 
signal  (we  arc  using  the  lock-in  in  the  broad  band 
mode)  given  by 

200 


er  (t) 


+1,  0  <  t  <  P/2 
-1,  P/2  <  t  <  P 


where  P  is  the  period  of  the  reference  and  also  the 
time  between  bias  pulses.  The  output  eQ  is  obtained 
by  taking  the  average  of  this  product  over  a  complete 
period.  Therefore, 


e„  = 


P  / 
& 


p  i  er(t)  es(t)  dt 
P/2 

p  .  (fle't/T  *  B) dt 

-  e-v/ur 


/'< 

P/2 


p-  /  (Ae‘t/T  +  B1DT 


(89) 


Notice  that  the  D.C.  component,  B,  is  removed  from 
the  first  half-period  by  the  negative  going  second 
half-period.  Prom  Figs.  67  and  68,  however,  the  A.C. 

-  t  /  <r 

signal,  Ac  ,  is  smaller  in  the  second  half-period, 
therefore,  we  obtain  a  net  positive  average  signal  for 
the  whole  period.  The  double  boxcar  operates  on  this 
same  principle. 

Since  we  wish  to  know  the  emission  rate  corres¬ 
ponding  to  a  DLTS  peak  we  note  that,  since 

I  -  cvNe-^kT 


then  dx/dT  /  0.  Therefor 
the  output  signal  for  a  p 
expression  yielding  this 

de°  dT  =  o 

TT 


c,  to  produce  a  maximum 
art i cul ar  temperature , 
max i mum , 


i  n 
the 


is  equivalent  to  de0/dx  = 
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0. 


To  find  the  signal  at  Dl.TS  max .  peak  height,  we  set 


de0/di  =  0  and  find  that 


Tmax  =  0*3981  =  0.398/f  (90) 

where  f  is  the  lock-in  reference  frequency.  Note  that 
this  expression  holds  for  every  peak  in  our  spectrum 
for  a  fixed  frequency. 

Substituting  this  relation  into  e0 ,  we  obtain 
the  peak  signal  height  in  terms  of  amplitude  A  as 


=  0.2036  A  =  e(mcter) 

°max  — j — y-f — 

Since  our  lock-in  is  calibrated  in  RMS  sinewave  units, 
but  measures  average  values,  the  factor  1.11  relates 
the  measured  output  to  the  meter  reading. 

From  the  above,  our  amplitude  is  given  by  A  = 

4.24  e(meter).  From  this,  AC  can  be  determined  as 


AC 


4.424 


[pf/volt  out] [mV  in/volt  out ][volts/inch][ peak  hgt .  in 

inches] 

capacitance  lock-in  recorder 

meter 


or 


AC  =  const  x  peak  hgt.  in  inches. 

To  make  our  production  rates  quantitative,  we 
must  investigate  the  frequency  response  of  our  system. 
We  have  done  this  for  several  DLTS  peaks  found  in 
neutron  irradiated  silicon  as  shown  in  Fig.  69. 


20  2 


J 


When  the  amplitudes  of  the  peak  heights  are 
normalized  to  the  maximum  peak  height  at  mid- frequency 


range,  the  frequency  responses  of  the  various  trap 
levels  are  quite  similar.  We  have  plotted  this  fre¬ 
quency  response  in  Fig.  69  and  find  that  our  peak 
heights  roll-off  at  both  the  high  and  low  ends  but 
are  a  maximum  between  ^  10  Hz  to  100  Hz.  Both  roll¬ 
offs  are  about  6  db/decade.  The  low  frequency  roll¬ 
off  is  just  the  lock-in  amplifier  low  frequency  res¬ 
ponse  while  the  high  frequency  roll-off  is  related 
to  our  gating  off  the  capacitance  signal  for  an 
appreciable  fraction  of  a  complete  period.  We  can 
conclude,  therefore,  that  quantitative  measurements 
can  be  made  of  amplitude,  aC,  using  this  frequency 


response  curve. 

The  analysis  of  the  DLTS  signal  is  quite  similar 


for  double  boxcar  detection  and  has  been  presented  by 
Lang.^  In  this  case,  the  output  signal  is  given  by 


in  the  limit  that  At  ->  0.  Again,  to  obtain  the 
emission  rate,  we  set  de0/di  =  0  and  obtain^4 

Tmax  =  Cti  -  t2 ) / [ In ( t j /t 2 ) ] .  (92) 


This  expression  is  analogous  to  Eq.  (90)  for  the 
lock-in.  In  essence  then,  in  both  experiments  we 
are  comparing  out  exponentially  decaying  signal, 
es(t),  to  a  fixed  time  gate,  either  the  period  P 
for  a  lock-in,  or  the  expression,  (tj  -  t£)  x  con¬ 
stant,  for  a  double  boxcar.  Equations  (90)  and  (92) 
then  tell  us  the  value  of  our  trap  emission  rate, 

Tmax>  at  the  DLTS  peak.  For  a  fixed  time  gate, 
either  P  or  tj  -  t2»  this  emission  -ate  is  known 
and  is  a  constant  for  each  peak  in  the  spectrum. 

By  changing  our  time  gate  (in  our  case  the  period) 
we  obtain  a  new  value  for  -rmax  for  each  peak  which 
occurs  at  a  shifted  temperature  from  the  peak  posi¬ 
tion  at  the  previous  gate  setting.  We  are,  in  effect, 
obtaining  a  table  of  values  of  imax  vs.  T  for  each 
peak.  Each  gate  setting  (period)  yields  another  pair 
of  coordinates.  Since 

—  =  o N v  e  ,  (9.i) 

7 

a  plot  of  In  tT  vs.  1/T  will  yield  a  straight  line 
slope  proportional  to  aE  and  an  intercept  proportional 
to  o.  These  two  parameters,  the  depth  of  the  energy 
level  and  the  majority  carrier  capture  cross  section 
characterize  our  defect. 

Typical  DLTS  spectra  for  room  temperature  irra¬ 
diated  Si  arc  shown  in  Fig.  7D  for  a  single  period  P. 

2  OS 


tCk) 


A  series  of  such  spectra  at  different  periods  yields 
the  Arrhenius  plots  shown  in  fig.  71  from  which  Ah 
and  a  can  he  determined. 


To  obtain  quantitative  results,  we  must  obtain 
the  pre-exponential  factor  temperature  dependence 
for  Eq.  (93).  Expressed  in  explicit  form  for  both 
n-  and  p-types 

°n  (p)  1//Tn(p)  °c (h )  vn(p)  Nc(v)  c 

where  aE  =  Ep  -  E-r  for  majority  traps  in  n-type  and 
&E  =  Ej  -  Ely  for  majority  traps  in  p-type.  The 
carrier  thermal  velocity  is  given  by  (cgs  units) 


The  density  of  states  is  given  by 

_  .  [  2tt  me  (h)  kT 
Nc(v)  -  -  1 - p - 

-  4.851  x  1015  fn>S  ( hi /m0>  3/2  T'V3  • 

(Combining  v  and  N  and  dropping  subscripts  we  obtain 

c  =  1  / t  =  oKT2  e'Al;/kl  (94) 

where  K  =  3.2  57  x  1  0“'*  (m  /m0)  . 

Using  mc  =  1.08  m0  and  =  0.59  mD  we  obtain 


MIS 


the  constants 


(95) 


8.518  x 

so21. 

n  -  type 

1  .  92  2  x 

«02i. 

p-type 

Once  AE  is  determined  from  the  Arrhenius  plot, 
it  is  useful  to  solve  Eq.  (94)  for  the  capture  cross 
section  which  becomes,  using  the  above  factors, 


tsj 

CO 

Xi. 

ho 

X 

2  0 "  2  2  gM:21<T/tT2 

( cm" ) 

5.204  x 

lO-22  0«i:/kT/tT2 

(cm2 ) 

(96) 


We  now  proceed  to  analyze  the  output  signals  for 
a  TCS  experiment  using  Eq.  (88)  realizing  that 
es(t)  =  i(t)R^  where  is  our  load  resistor  in  Fig. 
28.  We  will  begin  by  assuming  that  D.C.  leakage 
in  Eq.  (89)  is  negligible."’5  This  makes  it  possible 
to  do  the  experiment  with  a  single  channel  boxcar  as 
was  done  by  Wessel.5"’  We  shall,  therefore,  begin  with 
this  signal  analysis. 

A  single  channel  boxcar  consists  of  a  gated  am¬ 
plifier  which  turns  on  at  a  delay  time  t^  after  the 
bias  pulse  trigger  and  signal  averages  over  a  gate 
width  At  after  turn-on.  It  is  usually  assumed  that 
the  gate  width  is  a  small  fraction  of  the  total  decay 
time.  For  the  TCS  experiment,  the  signal  es(t  )  = 
i ( t ) R]  and,  using  Eq.  I  88),  we  out  a  in 
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Assuming  a  negligible  gate  width  (At-»0) 


To  make  a  measurement  of  the  emission  rate,  1  /i, 
we  need  an  expression  which  determines  the  TCS  peak 
position  as  a  function  of  temperature,  T.  Note  that 

deo  di 

dT~  3T” 


is  the  condition  for  a  maximum  in  e  (T)  ,  however,  since 
t (T)  has  no  maximum,  this  condition  becomes  dc0/dT  =  0. 
Taking  this  derivative  we  find  that  the  TCS  peak  occurs 
when  t  =  tj.  We  therefore  obtain  a  measurement  of 
the  emission  rate  en  =  1/t  =  1/t^  at  the  peak  tempera¬ 
ture.  By  generating  a  set  of  spectra  with  different 
delay  times,  the  capture  cross  section  and  trap  depth, 
F.y,  can  be  determined  from  a  plot  of  In(tjT-)  vs. 

1/T.  The  above  analysis  is  from  Ref.  55. 

Several  features  are  clearly  not  correct.  This 


result,  Tmax  =  tj,  is  based  on  F.q.  (88)  which  neglects 
the  leakage  current  due  to  Cp  in  bq .  (87).  When  we 
take  this  into  account,  we  obtain 


^-ANTRj.  (2e  ,  e 


Now  a  different  result  occurs  since  both  Cp  and  en  are 
temperature  dependent.  The  leakage  current,  if  signi¬ 
ficant,  can  therefore  shift  the  peak  positions. 


This  can  be  averted  by  using  a  double  boxcar 
amplifier  which  measures  at  delays  tj  and  t£  and  takes 


the  difference  between  the  two  signals,  A-B.  The 
time  independent  leakage  current  (a  DC  level)  obviously 
vanishes.  Therefore  the  double  boxcar  analysis  yields 


= 


qiv'ANTRL 
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^l/1  -  e'1^/1) 


(98) 


and  the  peak  position  can  be  obtained  by  solving  the 
equation  de0/dt  =  0  which  gives'^ 

'  i  „  e-(t2  -  tj)/. 


t2 


(99) 


To  find  the  peak,  we  need  to  solve  the  above 
equation  numerically  for  -rft}/^). 

Since  a  lock-in  is  also  a  signal  averager,  it 
can  also  be  used  for  I)LTS  and  TCS.  The  lock-in  output 
when  operated  broadband  is  given  by 
T 

^o  =  T  /  er  )  es  (t)  dt  (10  0 ) 

J0 


where  er  = 


+  1 
-1 


0  <  t  <  T/ 2 
T/2  <  t  <  T 


(101) 


if  we  set  the  phase  to  be  zero.  Similar  expressions 
can  be  derived  from  narrowband  operation  by  using  the 
Fourier  components.'^ 

In  a  sense,  the  lock-in  is  working  exactly  like 
a  double  boxcar  with  t^  =  0,  t?  =  T/2  and  At  =  T/2. 


The  gate  width  is  no  longer  negligible,  however,  the 
process  of  taking  the  difference  between  the  first 
and  second  half  periods  eliminates  the  leakage  term 
due  to  ep.  We  can  therefore  use  F.q.  (88)  for  the 
analysis.  We  obtain,  from  Fqs.  (100)  and  (101  I,  and 
the  signal  function 


(t)  -  i  e"t/i 


the  result 


e0  -  1  [,  ■ 

where  P  is  the  period.  Tor  DI.TS  we  had 
e0  f  [1  -  e‘r/2x]2. 


(10  2) 


(10a) 


The  character  of  these  two  equations  is  quite  different. 
Note  that  P.q .  (102)  has  no  max.  as  t  is  varied  but  Fq  . 
(103)  does.  Equation  (102)  tells  us  that  a  TCS  spec¬ 
trum  with  peaks  should  not  occur1  Figure  72,  however, 
shows  this  not  to  be  true  experimentally.  The  electron 
trap  peaks  Fj,  F?  and  F^  are  clearly  related  to  the 
DLTS  peaks  Fj,  [■-,  and  F-  obtained  previously.  There¬ 
fore,  something  is  clearly  wrong. 

Analysis  of  the  TCS  -  Lock-in  Problem 

From  the  above,  something  must  be  wrong  with  es(t) 
or  l(t)  or  ep.  Several  possibilities  exist: 

1)  fiateoff  Problem  -  flic  signal  is  gated  off 


during  the  bias  pulse  for  a  time  At  to  prevent 


overload.  The  output  should,  therefore,  be 
calculated  as 

Fo  =  'p  /  eserdt. 

At 

This  same  problem  exists  for  DLTS  and  has  been 
analyzed.  We  obtain  the  same  function  for  eD 
as  before  in  the  limit  At<<P.  We  have 
checked  this  effect  experimentally  and  have 
reduced  the  bias  pulse  width  until  the  peaks 
no  longer  shift. 

2)  Time  Dependence  of  Depletion  Width  -  The 
depletion  width  is  clearly  a  function  of  time 
because  this  is  the  only  time  dependence  found 
in  the  DLTS  experiment,  see  F.qs.  (82)  and  (84). 
This  means  that  W(t)  should  remain  under  the 
integral  contrary  to  popular  analysis  as  used 
above  and  in  References  55  and  56. 

We  have  obtained  the  proper  current  ex¬ 
pression  as  follows 

i(t)  =  -j-i  w(t)  [ 2ep  ♦  ene  t/T]  ;  1/t  =  en  . 

But  K(t>  -  ^ 

where  C(t)  -  C(-)  (1  -  IJJg). 

Therefore,  using  (83),  it  can  be  shown,  ne¬ 
glecting  leakage,  that 

i(t)  .  a.Aw(«) 


214 


N 
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N-p  e  t  /  t 

2TT 


*D 


(  104  ) 


where  W(°°)  = 


eA 


is  time  independent. 

Although  the  most  reasonable  solution  to 
our  problem,  it  should  be  noted  that  Eq .  (104) 
reduces  to  Eq.  (88)  where  Nj<<Np.  From  DLTS  we 
know  that  Nj  -c  10^  to  10^  and  %  10^  from 
breakdown  measurements. 

3)  Bias  Modulation  from  Load  Resistor  -  After 
the  bias  pulse,  the  depletion  width  also  changes 
because  the  bias  voltage,  V,  is  in  series  with 
the  load  resistor  voltage  i(t)R^.  Ke  can  also 
calculate  this  effect  in  an  approximate  way 
since 


wet)  = 


2  c  [  V 


Bi 


+  V 


i (t )RlJ 


q(Nn  -  nT(t)l 


1. 

'2 


W(-)  ( 


1  - 


(t)«I- 


2 (V  +  Vr i ) 


1) 

Substituting  into  i(t)  and  solving  for  i(t)  we 
obtain 
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This  last  equation  includes  all  effects  ex 


cept  leakage,  however,  it  is  completely  unmanage- 


able  mathematically.  It  is  possible  to  check  this 
effect  experimentally,  however,  since  varying  Rp 
should  produce  a  peak  shift  if  the  effect  is  im¬ 
portant. 

4)  Unknown  Theory  Problems  -  The  three  effects 
above  can  all  be  controlled  experimentally  by 
varying  At,  Ny  and  RT  to  the  point  that  they  can 
be  made  negligible.  The  analysis  of  these  effects 
is  based  on  the  fundamental  current  equation 


~9—  (enny  +  eppT) 


qhAN-p  /  2enep  ^  en(en  *  ep)  -X./i 
~T~  \en  +  ep  en  +  ep 


where  1/t  =  en . 


It  is  perhaps  possible  that  some  significant 

error  has  been  made  in  this  theory  in  Ref.  SO. 

In  fact,  an  alternate  equation  to  (105)  has  been 
5" 

given  as 

ito  =  q(x2  -  Xj)  [p  ennT  +  ( 1  -  1)  eppT] 

where  x ^  -  Xj  is  an  effective  generation  region 
smaller  than  IV  and  where  P  is  the  ratio  between 
the  excited  charge  within  x ->  -  X|  and  the  mea¬ 
sured  charge  in  the  external  circuit.  The  fern: 
of  this  equation  is  similar  to  ft).  (I  OS)  and 
although  the  quantitative  measurement  of  peak 


height  might  be  altered,  the  time  dependence  is 
expected  to  he  the  same. 

Although  the  TCS  experiment  shown  in  Fig.  72  is 
not  completely  understood,  one  aspect,  namely  spatial 
variation  of  defect  concentration  is  evident  and  the 
results  reasonable.  When  the  bias  pulse  amplitude  is 
varied,  the  width  of  the  modulated  depleted  zone  is 
also  modulated.  In  Fig.  72,  since  the  reverse  bias  is 
-10V,  those  bias  pulses  which  approach  OF,  are  sam¬ 
pling  traps  which  are  physically  located  near  the 
junction  in  addition  to  those  traps  at  greater  distances 
from  the  junction  which  are  observed  for  smaller  bias 
voltages.  It  is  clear  from  fig.  72  that  the  trap  li 2 
has  a  high  concentration  nearer  the  junction  than  Ej 
or  E3.  This  is  evident  because  of  the  increase  in 
amplitude  of  F. 2  at  high  bias  pulse  voltages  relative 
to  the  other  peaks. 

The  same  experiment  using  junction  capacitance 
docs  not  show  this  effect!  Lang  has  shown,  however, 
that  the  sensitivity  of  spatial  vaiiation  is  greatest 

C  Q 

for  TCS  near  the  junction  and  lowest  for  PITS. *  Wc 
must  therefore  conclude  that  the  defect  responsible 
for  level  Ft  has  migrated  during  irradiation  or  after 
toward  the  junction  as  a  result  of  the  electric  field 
of  the  junction.  Since  the  electric  field  in  the  n- 
type  depleted  side  of  the  junction  is  negative,  the 


pile-up  of  E2  at  the  junction  suggests  a  net  positive 
charge  state  for  this  defect. 


For  DLTS  to  be  completely  useful,  the  question  of 
resolution  should  be  investigated.  If  two  energy  levels 
of  similar  capture  cross  section  are  closely  space,  a 
theoretical  model  for  the  DLTS  line  width  is  required 
to  deconvolute  these  closely  spaced  peaks.  We  have, 
therefore,  attempted  to  fit  the  experimental  DLTS  spec¬ 
trum  for  neutron  irradiated  n-type  Si  (solid  line  in 
Fig.  73)  with  the  theoretical  function  for  lock-in  out¬ 
put  signal 


g_Ej/2ti^2 


(106) 


where  the  subscript  indexes  the  il  peak  of  ampl itude 
,  and  where  the  i1"  emission  rate  is  given  by 


°iveNce 


-E^/kT 


A  program  to  calculate  Fq.  (106)  as  a  function 


of  the  Pt  thermometer  resistance  R(T)  was  used  to  cal¬ 


culate  the  theoretical  curve  indicated  by  the  dots  in 


Fig,  73.  It  is  clear  that  the  experimental  line  widths 
are  greater  in  Fig.  73  than  the  calculated  line  widths. 
The  difference  between  experimental  and  calculated  is 
shown  in  the  bottom  of  Fig.  73  and  appears  to  be  a 
first  derivative  of  the  experimental  data. 

The  above  procedure  yielded  excellent  results  for 
peak  height  and  peak  position.  A  1"  shift  of  the  acti- 


CALCULATED 


Comparison  of  calculated  PITS  spectrum  (dots) 
vs.  experimentally  observed  spectrum.  Lower 
trace  is  the  difference  between  upper  snectra 


vat  ion  energy,  K  ^ ,  was  quite  apparent  when  comparing 
the  measured  and  calculated  results.  In  fact,  the 
experiment  was  repeated  to  obtain  r  XT*-  vs.  1/T  plots. 

The  least  squares  correlation  was  greater  than  99.8". 

This  precision  is  necessary  to  obtain  the  proper  peak 
positions.  The  small  remaining  least  squares  error  is 
the  result  of  measuring  peak  positions  on  the  X-Y  re¬ 
corder  plots. 

We  have  investigated  by  computer  calculation  various 
reasons  for  the  experimental  line  widths.  The  difference 
curve  suggests  that  amplifier  band  pass  might  be  affecting 
the  line  widths.  The  experimental  band  pass  is  limited 
in  several  respects.  We  must  first  consider  the  band 
pass  of  the  amplifiers  in  the  capacitance  meter.  The 
preamplifiers  before  the  lock-in  modulator  must  also 
be  considered.  Finally  after  modulation  (usually  called 
demodulation  in  lock-in  terminology)  the  effects  of  the 
low  pass  integrating  filter  must  lie  considered, 

The  band  passes  before  modulation  (product  with 
reference  signal)  can  be  analyzed  by  the  technique  of 
Laplace  transforms.  If  o  s  ( t )  =  Ac  t/7  is  the  signal 
before  low  or  high  pass  filtering,  then  the  transform 
of  this  signal  is  given  by 

CO 

Ms)  =  =  /  es^c'st  dt  =  f1""1 

0 

where  s  is  a  complex  variable  and  where  we  have  defined 
the  thermal  emission  rate  ]  /  =  u>  for  convenience. 

2  2n 


The  effects  of  amplifier  roll-off  can  be  analyzed 
by  the  Bode  plot  approximation  to  the  amplifier  as 
follows.  The  idea  is  to  find  an  approximate  amplifier 
transfer  function,  T(s),  in  the  complex  space  which 
represents  the  amplifier  response  effects  to  the  input 
function.  If,  for  example,  and  ojj j  are  the  low  fre¬ 
quency  -Sdb  points  of  the  amplifier  frequency  response, 
and  if  the  amplifier  has  a  6  db/octave  roll-off  at 
both  ends  of  the  response  and  a  mid-band  gain  of  G, 
then  the  exact  transfer  function  for  the  approximate 
Bode  plot  representation  of  the  frequency  response  is 
given  by 

T(s)  '  vr*  ULKs  nos, 

which  is  the  convolution  of  the  low  frequency  transfer 
function  s / ( s  +  wL)  and  the  high  f requency  transfer 
iunction  1  /  (s  +  u>||).  Faster  roll-offs  (12  db/octave, 

18  db/octave,  etc.)  can  easily  be  approximated  by  the 
product  of  more  terms  of  the  form  l/(s  +  wjj)  or  s/(s  + 
The  overall  effect  of  the  amplifier  response  on 
the  signal  function  is  then  given,  in  the  complex 
space,  by  the  convolution  of  the  response  function 
and  the  transfer  function  as 

E0(s)  =  T(s)  *  Es(s) 

_ G  s 

(s  +  ( S  ■*  U>]jyTs  +  b\) 


(109) 


The  function  HQ ( s  )  is,  however,  just  the  Laplace 
transform  of  the  amplifier  output  in  the  time  space. 
Therefore  the  output  function  for  our  exponential  in¬ 


put  is  given  by 

eD(t)  =  r1  {  E0(s)|  . 

The  Laplace  inverse  is  easily  found  for  Eq.  (109), 
since  it  may  be  broken  down  by  partial  fractions  (the 
advantage  in  using  approximate  Bode  plots).  The  result 
is  a  sum  of  three  exponential  functions  with  time  con¬ 
stants  r  =  1/w,  =  1/wp  and  x j.j  =  l/wq. 

This  method  can  also  be  used  to  analyze  the  effect 
of  the  high  pass  or  low  pass  of  the  amplifier  separate¬ 
ly.  For  example,  for  the  high  pass  case 

G 

M5)  =  (s  *  uH)  ( s  +  '  »') 

and 

vt)  -  n  i  e0cs)| 

=  ( _ £ — \  (e'U)t  -  e**0!!1)  .  (110) 

\UH  '  “/ 

Likewise,  for  the  low  pass  case  we  have 

eoO>  ■  (qph:)  C“c‘“t  -  •  oil) 

This  function  is  the  actual  time  dependent  function 
which  is  multiplied  by  our  square  wave  in  the  lock-in 


demodulator  and  then  averaged  over  a  period.  There¬ 
fore  the  actual  lock-in  output  is  given  by 
P 

=  J  J  e0(t)er(t)dt 
0 

if  we  assume  that  the  final  low  pass  filter  of  the 
lock-in  is  a  perfect  integrator.  [It  is  clear  from 
Eqs.  (110)  and  (111)  that  the  lock-in  output,  , 
contains  additional  powers  of  w  or  t  which  can  perliap 
explain  why  the  TSC  experiment  discussed  previously 
actually  has  peaks.]  We  have  calculated  the  DLTS 
spectrum  shown  in  Fig.  73  using  the  output  functions 
given  by  Eqs.  (110)  and  (111).  Neither  the  low  fre¬ 
quency  or  high  frequency  roll-off  appears  to  have  any 
effect  on  the  experimental  line  widt hs  sufficient  to 
explain  the  observed  effects. 

Although  not  conclusive,  because  of  the  approxi¬ 
mations  used,  we  feel  that  the  effects  of  band  width 
should  not  be  investigated  further  at  this  point. 

It  is  interesting  to  note,  however,  that  near  the  end 
of  the  frequency  response,  slight  shifts  in  peak  posi 
tion  are  tc  be  expected  which  could  alter  the  activa¬ 
tion  energies  and  cross  sections  slightly.  It  is  re¬ 
markable  to  us  that  this  analysis  has  not  been  inves¬ 
tigated  by  others  until  now.  Clearly  these  methods 
can  be  expanded  to  boxcar  amplifiers  as  well. 

As  a  final  point,  the  question  of  the  final  low 


pass  filter  integration  must  be  addressed.  If  the 


r 


|  integration  is  not  perfect,  then  the  lock- in  output 

j  signal  could  conceptually  lie  represented  as  a  mixture 

|  of  the  integrated  output  and  its  derivative  (non-inte- 

i 

grated  output).  In  fact,  such  a  guess  would  perhaps 

£ 

S  explain  the  difference  curve  of  Jig.  "3.  IV e  have  not 

|  pursued  this  problem  from  the  laplace  transform  point 

of  view,  however,  we  have  investigated  it  exper i ment a  1 1 
|  by  running  D1.TS  spectra  with  different  lock-in  time 

i 

|  constants.  No  effect  on  line  width  has  been  observed 

I 

J 

in  these  experiments  suggesting  that  an  analytical 
approach  is  not  justified. 

finally,  we  have  been  led  by  the  above  investiga¬ 
tion  to  the  following  question.  Can  a  DLTS  spectrum 
be  obtained  at  constant  temperature  by  varying  the 
frequency  or  gate  times?  Such  an  experiment,  which 
we  call  frequency  scanned  I'h'lS  (TSUbTS),  has  several 
advantages  from  a  radiation  damage  point  of  view.  If 
the  temperature  is  fixed  or  varied  only  over  a  small 
temperature  range,  defects  can  be  prevented  from 
annealing  making  the  technique  generally  useful  for 
low  temperature  experiments.  Secondly,  it  is  concep¬ 
tually  possible  to  build  a  room  temperature  HITS  appara 
tus  and  eliminate  the  cryogenics. 

Stimulated  by  the  above,  we  have  calculated  the 
I'l.TS  spectrum  of  fig.  "is  in  the  frequency  domain  at 
three  fixed  temperatures.  flic  results  arc  shown  in 


Fig.  74.  Although  the  peaks  appear  in  reverse  order 
in  the  spectrum,  they  are  easily  identified.  Unfor¬ 
tunately,  the  range  over  which  the  frequency  must  he 
scanned  in  order  to  cover  half  the  band  gap  of  Si  is 
extremely  large  (10  ^  Hz  to  10(l  Hz)  and  beyond  present 
day  electronic  capabilities.  It  is  interesting  to 
note,  however,  that  the  FIVHM  line  widths  on  a  log- 
frequency  scale  are  independent  of  both  temperature 
and  peak  height.  This  is  the  ideal  situation  we  have 
been  looking  for  in  resolving  closely  spaced  DLTS 
peaks.  The  FSDLTS  technique  can,  therefore,  find 
application  in  simple  deconvolution  of  DLTS  peaks  if 
the  temperature  of  the  sample  is  fixed  so  as  to  cause 
the  DLTS  peak  to  occur  within  the  band  pass  of  the 
amplifiers  used  to  detect  the  signals.  This  technique 
might  also  be  useful  in  analyzing  the  effects  of  am¬ 
plifier  band  pass  on  ideal  signals  in  that  a  deviation 
from  ideal  line  shape  at  t lie  low  or  high  frequency 
sides  of  the  peak  should  quickly  indicate  at  what  fre¬ 
quencies  amplifier  band  pass  effects  become  important. 

The  above  analysis,  although  lengthy,  is  necessary 
in  order  to  put  the  following  observed  data  into  proper 
perspective.  IVe  have  observed  a  variety  of  defect 
levels  in  neutron  irradiated  silicon  as  a  result  of 
room  temperature  irradiation  and  annealing  to  fiSO'T 
(seven  electron  traps  in  n-type  and  fifteen  hole  traps 
in  [)-  type  )  . 


5.32  *  10“13  cm 
3.50  *  10“,s  cm 


i 


Samples,  In  general,  are  the  col  lector- to-hase 
junctions  of  easily  obtained  commercial  transistors. 
This  was  convenient  in  a  general  survey  study  such  as 
this  because  of  the  large  number  of  samples  run.  loi 
the  p+n  (n-type)  data,  a  2N2270  was  used  while  a 
2N2904  was  used  for  the  n+p  (p-type)  samples.  All 
traps  listed  are  majority  carrier  traps,  i.c.,  elec¬ 
tron  traps  in  n-type  and  hole  traps  in  p-type.  Several 
Schottky  barrier  and  diffused  junction  diodes  were 
fabricated  on  phosphorus  and  boron  doped  silicon  in 
order  to  check  that  the  traps  obtained  were  not  im¬ 
purity-dopant  related.  Agreement  was  found  in 
general  with  the  results  obtained  on  transistors. 

We  have  been  fortunate  in  that  preliminary  reports 

5<) 

of  a  similar  experiment  by  C.  II.  Barnes'*  but  using  a 
boxcar  amplifier  system  have  been  available.  F.xccl- 
lent  agreement  on  the  room  temperature  energy  levels 
is  obtained  with  Barnes  data  which  gave  us  considerable 
confidence  in  our  more  extensive  results.  The  only 
other  neutron  irradiation  data  at  room  temperature  is 
indeed  sketchy  and  is  available  only  in  unpublished 
form  ("Al-'CRl.-TR-  76-  0024)  .  Our  data  and  those  of 
Barnes  arc  generally  not  in  agreement  with  the  data 
in  this  technical  report.  fable  12  lists  the  room 
temperature  production  rates  for  electron  and  hole 
traps  ( f.  |  ,  \.  ,  F-,  11 1  ,  1U.  !!-,  II  which  are  typ'uallv 


Table  12.  Production  rates,  capture  cross  sections  and 
energy  levels  of  electron  and  hole  traps  in 
neutron  irradiated  Si.  Data  of  C.  E.  Barnes 
is  taken  from  Ref.  59. 
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observed.  Several  p-tvpe  samples  had  a  different 
hole  trap  spectrum  (H|',  I  i  ,  11-'  and  11  j  '  )  ,  however, 
modest  annealing  to  higher  temperatures  made  the 
spectra  of  the  two  types  of  samples  similar.  IV  e 
therefore  believe  the  effects  to  be  related  to  impurity 
defect  interactions. 

All  samples  were  irradiated  in  RSA!:  12-16  in 
boron  walled  cans  to  reduce  activation  of  the  Au  in 
the  devices  by  thermal  neutrons.  Several  samples  were 
irradiated  with  thermal  neutrons  as  well,  with  little 
difference  in  the  spectra  observed.  fluences  for  the 
production  rate  measurements  of  Table  12  ranged  from 
about  2  x  id15  to  2  x  1 0 1  (’nt j-,/cm"  for  n-type  material 
and  from  2  x  1  (I  to  2  x  10  Tt^/cm"  for  p-tvpe 
material.  The  fast  neutron  fluences  are  estimated  to 
be  about  1/3(1  of  the  thermal  neutron  fluences  in  RSAi'. 
The  production  rates  in  Table  11  were  based  on  this 
estimate.  Tt  should  lie  noted  that  the  agreement  with 
ratios  of  production  rates  with  the  data  of  Barnes  is 
good,  however,  we  report  a  considerably  lower  absolute 
rate  for  all  levels  in  p-tvpe.  Reasons  for  these  dif¬ 
ferences  are  probably  related  to  the  nature  of  the 
fast  neutron  component  in  these  two  reactors. 

Isochronal  annealing  of  the  PI. IS  spectra  in  both 
n-type  and  p-type  material  has  been  completed  to  Suit'd'. 
Diode  failure  occurred  at  (din'd'  in  the  2X2..'"'  n-tvpe 


transistor  and  at  S50°C  in  the  2N2904  p-tvpe  transis¬ 
tor.  Figure  75  shows  the  n-type  spectra  after  each 
15  min.  anneal  (the  anneal  temperature  is  near  the 
spectrum  while  the  DLTS  scan  temperature  is  at  the 
bottom  of  this  figure).  A  total  of  seven  electron 
trap  levels  were  observed.  The  absolute  concentrations 
of  these  levels  is  shown  in  Fig.  76  as  a  function  of 

anneal  temperature.  The  thermal  flucnce  on  this  p+n 

15  2 

diode  ivas  1.8  x  10  n/cm  .  The  additional  phosphorus 

11  -  5 

produced  by  this  irradiation  is  below  3  x  10  cm  '  or 
about  an  order  of  magnitude  below  the  lowest  defect 
concentration  observed. 

The  defect  energy  levels  and  capture  cross  sec¬ 
tions  are  shown  in  Table  13  for  the  p+n  junction  of 
transistor  T - 1 5  (2N22  70). 
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Figure  75.  HITS  electron  traps  in  n-lype  ncutrrr  irrn- 

d  i  ate  d  Si  ns  u  function  of  i  sqchrcr.u  i  nnr.cals 

non  n?_,  <tt],  =  i.s  x  m1  n-i/cmn . 


CONCENTRATION 


Figure  76. 


Concentration  of  electron  trap  levels  vs. 
isochronal  annealing  temperature  (data  from 


It  can  be  observed  in  Fig.  75  that  the  peak 
position  of  the  H ^  peak  shifts  slowly  to  the  right 
as  annealing  temperature  is  increased.  F.ach  shift 
is  accompanied  by  a  drop  in  defect  concentration  (Fig. 
76).  We  have  designated  the  shifted  peaks  E-(I), 
H^(II)  and  F 3  (1 1  I )  .  No  energy  level  data  was  obtained 
for  however,  F^(II)  =  0.42  eV ,  2.5  x  10  '  ''em" 

which  is  within  experimental  error  of  the  values  for 
=  0.42  cV,  5."  x  10  ''t’cm“.  The  shifts  in  peak 
{•-.(ill)  are  caused  for  the  most  part  by  changes  in  the 
value  of  the  carrier  capture  cross  section  as  shown  in 
Ta h 1 c  11. 

Table  14.  Ionization  energies  and  capture  cross 
sections  for  the  Infill)  level. 

Anneal  Energy  of  Capture 
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0. 02)eV 
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slightly  higher  than  that  of  1-  = 

in.42  • 

0.02) cV.  Thi 

s  suggests  that  peak 

F^  is  composed 

o  f  a  t 

least  two  very  closely  spaced  energy  levels  which 
anneal  at  different  temperatures. 


New  peaks  which  are  produced  by  annealing  n  r 
the  I-4 ,  1:5,  h^  and  by  levels.  The  F.^  peak  lias  been 
observed  by  Barnes.^  This  peak  grows  rapidlv  from 
250  to  350°C  and  then  anneals.  The  bj  peak  appears 
as  the  peak  anneals  suggesting  a  relationship  be 
tween  these  two  defects.  There  is  also  some  sugges¬ 
tion  from  the  shape  of  the  peak  that  h 5  might  also 
be  a  double  level.  At  still  higher  temperatures 
(500- 5 50 °C)  the  b^  and  by  levels  appear. 

Figures  77  and  78  show  the  DLTS  spectra  for  iso- 
chronally  annealed  p-type  material  (transistor  T-b, 

2 N 2 9 0 4  ) .  This  sample  was  irradiated  to  a  thermal 
flucnce  of  1.8  x  10  ln./cm",  a  factor  of  ten  higher 
than  the  flucnce  on  the  n- type  sample.  Two  separate 
reference  gate  frequencies  are  shown  (10  Ilz  for  big. 

77  and  500  Hz  for  Fig.  78).  The  50  0  II z  scan  is  in¬ 
cluded  to  show  peak  II, fable  15  lists  the  energy 
levels  and  capture  cross  sections  for  the  p-type 
hole  traps. 

The  isochronal  annealing  of  the  p-type  sample  is 
shown  in  Figs.  70  and  SO.  The  data  has  been  eparatcu 
into  two  figures  for  clarity.  Of  the  levels  shown, 
the  original  levels  1! ->  and  II-  or  11^  have  been  observed 
previously  by  Barnes'  a.s  well  as  !!()  and  perhaps  II-.. 

The  annealing  of  this  sample  is  truly  complicated 
and  deserves  much  further  study.  The  beha\  inr  of  ll7 
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Table  15.  bnergy  levels  (by  -  by  I  and  carrier 

capture  cross  sections  for  hole  traps 
in  p-type  material. 
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is  remarkable.  The  concentration  of  this  defect  drop 
an  order  of  magnitude  between  300  and  350°C  and  then 
rises  again  at  450°C,  then  disappears  at  500 "C. 

This  defect  could  well  be  associated  with  the  "600°C 
defect"  seen  with  resistivity  measurements.  It  may 
be  that  we  have  the  same  energy  level  and  similar  cap 


ture  cross  sections  for  different  defects. 
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Identification  of  the  defect  structures  associated 
with  these  energy  levels  is  difficult.  In  Fig.  81  we 
show  neutron  damage  EFR  annealing  data  from  Ref.  22, 
Unfortunately,  the  times  and  temperatures  in  the  Purdue 
experiment  were  not  constant  and  we  have  observed  dif¬ 
ferences  in  the  annealing  from  data  shown  in  Fig.  81. 

Guldberg61  has  recently  summarized  most  of  the 

annealing  facts  known  about  neutron  irradiated  Si. 

Although  this  paper  shows  DLTS  data  for  n-type  material, 

his  annealing  does  not  start  until  100°C.  Very  few 

correlations  exist  betw'ecn  our  data  and  his  with  the 

exception  of  E n  which  Guldberg  claims  is  a  surface 

state.61  When  we  compare  our  annealing  data  with 
5  7 

Barnes,'  however,  the  agreement  is  good.  We  believe 
this  to  be  due  in  part  to  the  fact  that  we  and  Barnes 
have  investigated  CZ  while  Guldberg’s  work  is  in  FZ. 

Corbett1”  has  also  reviewed  know'n  annealing 
stages  in  Si.  His  review  is  similar  to  that  of 
Guldberg  ’  s  .  6 1  Kimcrling6'^  has  rcvicw'ed  DLTS  data  for 
electron  irradiation.  Some  similarities  occur  between 
these  data  and  ours.  From  all  of  these  sources  we 
can  make  some  educated  guesses  as  to  the  nature  of 
t lie  DLTS  defects  involved. 

Our  Ej  level  appears  to  be  well  established  as 
the  vacancy  -  oxygen  A  -  ccn  ter .  1 1  ’  6 It  should  be  noted, 
however,  that  this  level  does  not  anneal  as  soon  as 
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Isochronal  annealing  of  various  HPR  signals  in 
neutron  irradiated  Si  (data  from  Ref.  ?21  . 


2  1 1 


figure  81. 


reported  in  Refs.  59  and  64.  There  is  some  confusion 
about  the  F-^  level  at  different  annealing  temperatures, 
however,  the  first  annealing  stage  of  F-  at  100°C  is 
probably  related  to  the  breakup  of  the  F.-center  va¬ 
cancy  phosphorus.  From  250-400°C  the  level  F-(IIl) 
appears  to  be  a  divacancy  level  as  well  as  the  1 h  and 
H7  levels. All  of  these  anneal  in  about  the  same 
temperature  range. 

The  growth  of  F^  could  perhaps  be  related  to  the 
FPR  P-1  which  is  a  five  vacancy  cluster  which  increases 
as  the  divacancy  breaks  up.  The  rise  in  II,  should  then 
be  related  to  the  growth  of  the  FPR  signal  P-4  which 
is  known  to  contain  oxygen  and  is  thought  to  be  a  tri- 
vacancy  oxygen  center. ,f’“  The  rise  in  F5  in  the 
same  range  but  to  lower  maximum  concentration  is  most 
likely  to  be  associated  with  the  P-2  center  which  is 
thought  to  be  a  di vacancy-dioxygcn  center,^’1’'  The 
other  DI.TS  lines  at  high  annealing  temperatures  probabl 
involve  oxygen,  however,  little  more  can  be  said  about 
them.  The  low  temperature  II-  line  is  similar  to  the 
P-3  four  vacancy  center  while  the  behavior  of  II 5 
appears  to  be  related  to  the  P-6  resonance  which  is 
thought  to  be  a  d  i  -  in  t  e  rst  i  t  i  a  1  .  fl '  ’  “ 

We  have  summarized  our  best  guesses  as  to  the 
nature  of  the  HI,  IS  levels  in  Table  16.  these  guesses 
have  been  based  mainly  on  annealing  temperatures,  tala 


tive  concentrations  of  known  defects  in  neutron  irra¬ 


diated  Si,  and  in  the  case  oh  the  A  and  V-Y  centers, 
on  other  DLTS  papers.  It  must  be  realized  that  the 
situation  is  extremely  complicated  and  requires  much 
further  study. 

Table  16.  Tentative  identification  of  DLTS  levels 
with  known  defect  structures  in  Si. 
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V.  PROGRESS  AND  ACCOMPLISHMENTS 

The  objective  of  this  research  program  has  been  to 
extend  and  advance  the  technique  of  nuclear  transmutation 
doping  in  silicon  and  to  advance  the  state  of  the  art  in 
compensating  residual  boron  in  IR  extrinsic  detector 
material.  We  summarize  briefly,  here,  the  main  accomplishments 
of  this  program. 

•  The  theory  of  the  behavior  of  the  room  temperature 
electrical  parameters  in  NTD  silicon  vs.  fluence  has 
been  completed  and  accepted  for  publication.  A 
variety  of  intrinsic  parameters  have  been  found  which 
relate  the  extrema  of  these  curves  to  a  single  fluence 
parameter,  the  critical  fluence.  The  calculation  of  the 
critical  fluence  parameter  has  been  extended  to  all 
p-type  starting  material  and  is  now  an  exact  calculation, 
not  an  approximate  one. 

•  A  resistivity  homogeneity  model  has  been  completed 
and  accepted  for  publication  which  provides  a  calculation 
of  the  behavior  of  Ap/p  vs.  NTD  phosphorus  to  boron 

ratio.  This  theory  shows  that  the  resistivity  inhomogeneity 
is  no  worse  and  generally  better  than  the  starting 
material  inhomogeneity  for  phosphorus  to  boron  ratios 
greater  than  two. 

•  A  model  has  been  completed  which  allows  calculation 
of  the  maximum  resistivity  which  can  be  obtained  by  NTD 
compensation  of  p-type  silicon  before  inhomogeneity 
fluctuations  cause  mixed  typing  of  the  material.  This 
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model  should  be  useful  for  a  variety  of  other  envisioned 
detector  devices  which  depend  on  high  resistivity  for 
their  operation. 

A  new  calculation  of  the  damage  rate  due  to  fast 
neutron  displacements  has  been  made  which  takes  into 
account  the  neutron  flux  spectrum  and  secondary's  ioni¬ 
zation  for  the  first  time  and  is  relevant  to  row-2 
of  the  reactor,  the  position  where  EPR  and  optical 
samples  have  been  irradiated.  The  displacement  rate 
is  found  to  be  2  x  10^  displacements/cm^/sec  for  row-2. 

A  method  of  determining  In-concentration  from  room 
temperature  resistivity  measurements  is  presented  which 
should  be  accurate  to  better  than  101.  This  technique 
has  been  verified  by  neutron  activation  analysis. 

The  radioactive  decay  of  114mIn  has  been  investigated. 
It  is  found  that  the  case  of  heavy  In  doping  and  heavy 
irradiation  as  in  CZ-Si  leads  to  minor  activity  problems, 
however,  the  case  of  lighter  In  doping  and  lighter 
irradiation  as  in  FZ-Si  produces  no  activity  problems. 

The  vertical  flux  profile  in  RSAF  has  been  determined 
and  axial  uniformity  has  been  shown  to  be  adequate 
for  ingot  doping  over  restricted  lengths.  The  major 
limitation,  here,  will  be  axial  boron  uniformity  in 
the  starting  material. 

Ingot  and  wafer  programs  are  in  progress  with  all 
the  major  industrial  producers  of  IR  detector  arrays 
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using  NTD  compensated  Si.  General  experience  shows  that 
phosphorus  concentration  can  be  controlled  to  *3"  using 
industrially  determined  data  and  annealing.  We  believe 
this  to  be  an  extremely  significant  accomplishment 
since  it  moves  NTD  compensation  out  of  the  research 
laboratory  environment.  Furthermore,  minority  carrier 
lifetimes  of  550  usee  have  been  reported  doing  NTD 
compensation  under  these  circumstances. 

The  survey  of  isochronal  annealing  of  resistivity 
in  NTD  silicon  has  been  completed  for  high,  medium, 
and  low  fluences  in  C Z  and  FZ.  These  results  have 
been  published. 

Preliminary  experiments  have  been  completed  which 
suggest  that  oxygen  concentration  can  be  determined 
from  resistivity  changes  which  occur  during  annealing 
in  the  450-550°C  temperature  range  in  high  resistivity 
float  zone. 

A  number  of  defects  introduced  by  NTD  irradiation 
have  been  identified  by  FPR  and  optical  absorption  in¬ 
cluding  the  divacancy,  the  vacancy -oxygen  A  center,  the 
vacancy -phosphorus  ii-center,  the  four  vacancy  chain  P-3, 
the  five  vacancy  cluster  P-1,  one  of  the  HOB  bands  and 
the  phosphorus  concentration  related  free  carrier  reso¬ 
nance.  Most  of  these  defects  have  been  measured  in  a 
quantitative  way  to  some  extent  and  a  count  of  the  number 
of  vacancies  contained  in  these  defects  compared  to 
displacement  theory.  The  results  suggest  that  approx i- 


mutely  6°i  of  the  displacements  by  fast  neutrons  remain 
stable  following  room  temperature  NTH  doping.  The 
surprising  result  that  the  concentration  of  the  four 
vacancy  P-5  center  is  unaffected  by  interstitial  oxygen 
gettering  and  A-center  formation  suggests  that  this 
defect  is  a  primary  radiation  damage  product  which  is 
stabilized  during  the  displacement  cascade  and  not  formed 
from  a  general  vacancy  gas  by  coalescence. 

Raman  scattering  has  been  used  to  detect  highly 
disordered  regions  in  heavily  neutron  irradiated  Si. 

The  neutron  damage  induced  Raman  peaks  have  much  in 
common  with  those  Raman  peaks  observed  in  amorphous 
Si  and  tend  to  anneal  in  a  manor  characteristic  of 
amorphous  Si  regrowth.  The  frequency  shifts  of  these 
peaks  upon  annealing  have  been  interpreted  in  terms 
of  Brillouin  zone  unfolding.  As  an  added  benefit,  it 
was  discovered  that  boron  local  modes  could  bo  easily 
studied  in  neutron  irradiated  Si  because  the  radiation 
damage  compensated  the  free  carrier  interfering  effects 
normally  observed  in  such  material. 

DLTS  measurements  on  neutron  irradiated  silicon 
have  yielded  a  wealth  of  defect  energy  levels  in  the 
band  gap.  Upon  annealing  to  500°C,  fifteen  hole  traps 
and  seven  electron  traps  have  been  identified.  Although 
some  similarities  to  HPR  neutron  data  and  electron 
irradiation  data  occur,  the  situation  is  extremely 
complicated  and  certainly  deserves  much  further  study. 


It  is  also  clear  that  new  diode  fabrication  techniques 
will  have  to  be  devised  to  break  the  500°C  barrier  of 
thermally  induced  device  failure.  We  have  devised 
several  such  techniques  which  should  solve  this  prob¬ 
lem. 

We  feel  that  continued  progress  on  the  investigation 
of  specific  defect  production  rates,  annealing,  and  defect 
concentration  ratios  as  a  function  of  reactor  spectrum  will 
form  the  foundation  for  a  better  understanding  of  the  NTD 
process.  These  techniques  can  perhaps  be  used  to  determine 
the  overlap  of  defect  cascades  and  may  eventually  provide 
information  which  will  explain  the  observation  that  a  hard 
neutron  spectrum  produces  higher  minority  carrier  lifetimes. 


VI.  APPENDIX 


During  the  publication  of  this  report  additional  infor¬ 
mation  about  the  nature  of  the  displacement  cascade  has 
been  discovered  which  we  will  outline  briefly  in  this 
appendix. 

For  the  experiment  discussed  in  Section  IV-E,  where 
the  defect  concentrations  observed  by  optical  absorption 
and  EPR  are  listed  on  pp.  164-166,  it  is  instructive  to 
compare  these  concentrations  with  the  concentration  of  cas¬ 
cades.  For  this  experiment,  the  displacement  concentration 
19  -  3 

was  4  x  10  cm  while  the  mean  number  of  displacements 
per  neutron  was  181  from  Table  3.  The  mean  concentration 
of  cascades  is  therefore  given  by 

[cascade  concentration]  =  Np/v  =  2.2  x  10  '  cm  * 

Since  the  mean  recoil  energy  of  the  secondary  making  this 
cascade  is  ^  22  keV,  then  the  range  of  this  recoil  is  on 

O 

the  average  about  125  A.  Assuming  that  the  radius  of  this 
cascade  track  is  no  larger  than  the  range,  it  is  easy  to 
show  that  the  overlap  of  cascades  is  not  likely  for  this 
fluencc.  It  should  also  be  noted  that  the  concentration 
of  cascades  is  almost  identical  to  the  concentration  of 
P-3  cascade  remnants.  In  other  words,  there  are  about  one 
or  two  four-vacancy  P-3  centers  created  in  each  cascade. 

To  understand  this,  we  must  consider  the  energy  loss 
of  the  primary  knock-on  (PKA)  ion  along  the  displacement 
track.  As  the  ion  moves  along  the  track  displacements  are 


made  producing  simple  Irenkel  pairs.  As  the  recoil  secon¬ 
dary  energy  is  lost  hv  this  process,  the  ion's  mean  free 
path  between  collisions  drops  to  a  value  below  the  separa¬ 
tion  distance  between  atoms  along  the  track.  Once  this 
happens,  every  atom  along  the  track  must  be  displaced 
creating  a  Brinkman  displacement  spike.  An  estimate  of 
the  energy  below  which  a  spike  is  created,  fis  ,  is  given  by1 
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=  9  57  e  V 
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where  Ea  »  20  eV  ( Z2 ) 7/6  =  9448  eV  and  a  =  a  /21/5  = 

-  9  ° 

2.2  x  10  '  cm  where  a0  =  0.53  A  is  the  Bohr  radius.  This 
energy  should  be  compared  with  the  mean  energy  per  dis¬ 
placement  =  E/37  =  122  cY  and  with  the  energy  required  to 
make  a  four- vacancy  center  which  can  be  obtained  roughly 
f  rom 


E(P-3)  =  2Hdv  =  (50  cY)  v  =  200  eY  ,  (2-A) 

with  \T  =  4 .  Considering  the  approximations  of  the  above 
models,  and  the  possibility  of  some  defect  recombination, 
the  agreement  between  liqs.  (1-A)  and  (2-A)  is  reasonable. 

We  conclude,  therefore,  that  our  22  kcY  primary  ion  loses 
energy  along  the  track  by  creating  ionization  and  single 
displacements  at  a  rate  cf  about  122  cY  per  Irenkel  pair. 

Once  the  spike  energy  is  reached,  a  track  cf  displace 
ments  results  in  the  formation  of  a  few  P-3  four-vacancy 
centers.  Since  the  P-(>(IX)  interstitial  defect  complex, 
which  we  have  recently  observed,  occurs  in  equal  numbers 
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with  P-3  regardless  of  the  reactor  spectrum,  it  most  likely 
is  formed  from  the  atoms  ejected  from  the  displacement 
spike  also.  In  fact,  the  very  low  concentration  of  P-3 
and  P-6  centers  observed  in  high  Cd-ratio  samples  irradiated 
at  Argonne's  CP- 5  reactor  can  be  accounted  for  quant  i  t at i ve 1 
by  the  contamination  of  the  thermal  flux  by  a  small  fast 
neutron  component,  indicating  that  a  knock-on  Si,  with  an 
energy  less  than  or  of  the  order  of  a  typical  gamma  recoil 
knock-on  (C  -  474  eV) ,  has  a  negligible  probability  for 
producing  a  large  cluster  center  like  the  P-3.  This  lattice 
disorder  is  not,  however,  isolated  as  recognizable  isolated 
point  defects.  Only  t lie  fast  neutrons  accomplish  this  by 
way  of  the  spike  mechanism. 

It  should  also  be  noted  from  p.  164  that  there  arc  also 
about  two  divacancies  per  cascade.  It  is  not  known  at  this 
time  if  the  divacancy  is  also  a  result  of  the  displacement 
spike,  however,  the  energy  required  to  form  the  divacancy 
is  from  hq.  (2- A)  about  100  cV.  The  sum  of  the  divacancv 
and  four-vacancy  energies  is  an  appreciable  fraction  of  the 
total  spike  energy  (2  x  2 00  eV  +  2  x  100  eV) . 

Ihe  five-vacancy  center,  whose  concentration  never 
exceeds  10?,  of  the  four-vacancy  concentration,  does  not 
appear  to  be  a  cascade  remnant.  The  structure  of  this 
defect  (I  ig.  33)  is  non-planar  as  opposed  to  the  divacancy 
and  four  -  vacancy  centers  which  appeal'  as  tracks  in  |llo| 
directions.  be  believe  that  this  defect  is  formed  strict  1  v 
by  accident  whenever  four- vacanc  ics  and  <!  i  v.umiic  i  o~ 
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on  adjacent  sites  by  the  collapse  of  a  divacancv  and  four- 
vacancy  into  a  five-vacancy  and  free  vacancy.  The  isolated 
free  vacancy  can  then  migrate  to  interstitial  sites  for 
annihilation,  leaving  a  stable  f i ve- vacancy  center  to  be 
observed.  This  model  is  based  on  the  fact  that  the  con¬ 
centration  of  f ive- vacanc  ies  is  low  and  that  linear  live- 
vacancy  tracks  have  not  been  observed  by  f-PR. 
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